CERN-PH-TH/2004-139 
LMU 09/04 
|hj)-ph/0407244l 



MSSM Higgs Physics at Higher Orders 



S. Heinemeyer 1,2 * 



1 CERiV, Theory Division, Dept. of Physics, 1211 Geneva 23, Switzerland 

2 Institut fur theoretische Elementarteilchenphysik, LMU Miinchen, Theresienstr. 37, 

D-80333 Miinchen, Germany 



Abstract 

Various aspects of the Higgs boson phenomenology of the Minimal Supersymmetric Standard 
Model (MSSM) are reviewed. Emphasis is put on the effects of higher-order corrections. The 
masses and couplings are discussed in the MSSM with real and complex parameters. Higher- 
order corrections to Higgs boson production channels at a prospective e + e~ linear collider are 
investigated. Corrections to Higgs boson decays to SM fermions and their phenomenological 
implications for hadron and lepton colliders are explored. 
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Chapter 1 
Introduction 



The search for the lightest Higgs boson is a crucial test of Super symmetry (SUSY) [1] which 
can be performed with the present and the next generation of accelerators. The prediction 
of a relatively light Higgs boson is common to all supersymmetric models whose couplings 
remain in the perturbative regime up to a very high energy scale [2]. Within the Minimal 
Supersymmetric Standard Model (MSSM), where all parameters are assumed to be real 
(rMSSM), the mass of the lightest CP-even Higgs boson, m^, is bounded from above by 
about m h < 140 GeV 1 [5,7], including radiative corrections at the one-loop [8-11] and at 
the two-loop level [5, 7, 12-25]. In the case where also complex parameters are allowed 
(cMSSM) the evaluations of the Higgs boson sector are less advanced [26-33], but the upper 
bound of ~ 140 GeV still holds. This places the lightest MSSM Higgs boson in the reach of 
the currently operating Tevatron (depending on the luminosity performance), the LHC, and 
a prospective future e + e~ linear collider (LC). 

While at the Tevatron one can at most hope for the discovery of the Higgs boson and 
possibly a crude mass measurement [34], at the LHC already a high-precision determination 
of nih down to Srrf^ 9 = 100 — 200 MeV [35] as well as some coupling and total width 
measurements [36] seem to be feasible. At the LC eventually, rrih can be determined at 
the 50 MeV level [37]. The mass and width measurements are summarized in Tab. Il.ll for a 
Higgs boson with ~ 120 GeV. The various production cross sections [37,38], see Tab. 11.21 
as well as the branching ratios to SM fermions and gauge bosons [37,39,40], see Tab. 11.31 
can be determined with high precision down to a few per cent at the LC. Even the trilinear 
Higgs boson couplings seem to be in reach [41]. 

These expected accuracies make it mandatory to have a corresponding precision of the 
theoretical predictions in terms of the relevant SUSY parameters at hand. In this report we 
review several recent theoretical evaluations of higher-order corrections to the Higgs boson 
masses in the MSSM, to the dominant production cross sections at the LC, and to the decays 
to SM fermions, which are relevant for measurements at the Tevatron, the LHC, and the 
LC. 



1 This value has been obtained with the top-quark mass value of n\ xp = 178.0 GeV [3] that has become 
available while finalizing this report. Most results in this report, however, have been derived with the former 
value m° xp ' old = 174.3 GeV [4]. The use of the new value is always indicated. The main effect of the higher 
mj Xp value is an increase in the nih prediction by ~ 3 GeV [5,6]. 
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collider 


Sm h [GeV] 


sr h /r h 


Tevatron 


0(2) 




LHC 


0.2 


40% 


LC 


0.05 


6% 



Table 1.1: Expected precisions for the measurement of the lightest Higgs boson mass and 
width at the Tevatron, the LHC and the LC [34,35,37]. 



collider 


decay mode 


5 a/a 


LC 


e+e- -> Z* -> Zh 


1.5% 


LC 


e + e~ — > vvW + W~ — > z/z//i 


2% 



Table 1.2: Expected precisions for the measurement of Higgs production cross sections at 
the LC [37,38]. 



decay mode 


5BR/BR (v^ = 500 GeV) 


(ffiR/BR (Vs = 1 TeV) 


h^bb 


1.5% 


1.5% 


h — > T + T~ 


4.5% 


2% 


/l — > CC 


6% 




gg 


4% 


2.5% 


h -> V^Vy* 


3% 





Table 1.3: Expected precisions for the measurement of Higgs branching ratios at the LC for 
Higgs boson masses of ~ 120 GeV. For the analysis at y/s = 1 TeV not all channels have 
been investigated [37,39,40]. 
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Chapter 2 

The Higgs boson sector of the MSSM 
with real parameters 



In this section we will concentrate on the corrections in the Higgs boson sector in the MSSM 
with real parameters (rMSSM). The corresponding case with complex parameters (cMSSM) 
is treated in Sect. EJ An introduction to other ascpects of MSSM phenomenology can be 
found in Ref. [1]. 



2.1 The Higgs boson sector at tree- level 

Contrary to the Standard Model (SM), in the MSSM two Higgs doublets are required. The 
Higgs potential [42] 

V = ml\Hi\ 2 + m 2 2 \H2\ 2 - m 2 12 (6 ab H1H b 2 + h.c.) 

+\{9l + 9l) [m 2 -\n 2 \ 2 ] 2 + \gl\n\n 2 \ 2 , (2.1) 

contains m 1 ,m 2 ,m 12 as soft SUSY breaking parameters; g,g' are the SU{2) and U(l) gauge 
couplings, and ei 2 = — 1. 

The doublet fields Hi and H 2 are decomposed in the following way: 



H2 ~ U° ) - ) • (2 - 2) 

The potential (j2.1j) can be described with the help of two independent parameters (besides 
g and g')\ tan/5 = v%jv\ and M\ = — m^ 2 (tan/5 + cot @ ), where M4 is the mass of the 
CP-odd Higg boson A. 

The diagonalization of the bilinear part of the Higgs potential, i.e. of the Higgs mass 
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matrices, is performed via the orthogonal transformations 



h° 
G° 
A 
G ± 



cos a sin a 

— sin a cos a 
cos (3 sin (3 

— sin (3 cos (3 
cos (3 sin /? 

— sin (3 cos /3 



0? 



02 



The mixing angle a is determined through 



a = arctan 



-(Ml + M|)sin/3cos/3 



M| cos 2 /? + M 2 sin 2 /? -m 2 tree 

One gets the following Higgs spectrum: 

2 neutral bosons, CP = +1 
1 neutral boson, CP = — 1 
2 charged bosons 
3 unphysical Goldstone bosons 



7T 



- < a < 



h,H 
A 

H + ,H- 



M 



2, tree 
Higgs 



M\ sin 2 (3 + M 2 Z cos 2 /3 — (M| + Mf) sin/3 cos /? 
-(Mi + Mf) sin /? cos /3 Mi cos 2 /3 + M| sin 2 (3 



(2.3) 
(2.4) 
(2.5) 



(2.6) 



(2.7) 



At tree level the mass matrix of the neutral CP-even Higgs bosons is given in the <p\-(j)2- 
basis in terms of Mz, M4, and tan/3 by 



{21 



which by diagonalization according to eq. ()2.3|) yields the tree-level Higgs boson masses 



M 



2, tree a 



m 



if, tree 



m 



/i.trec 



Higgs ' \ Q 

The charged Higgs boson mass is given by 



m 2 H± M'\ ■ .Uf r . 
The masses of the gauge bosons are given in analogy to the SM: 



M 2 W 



■ 2/ 2 1 2\ 



M 2 Z 



1 



2 2 ); 



M^ = 0. 



(2.9) 



(2.10) 



(2.11) 
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2.2 The scalar quark sector 



The squark mass term of the MSSM Lagrangian is given by 

r 

*-"m. 



where 
Z 



M% + M\ cos 2ft (l( - Qfsl) + m) m f (A f - //{cot ft ; tan/5 }) 



(2.12) 



rrif(Af — //{cot /3; tan/?}) 



M| ; + Mf cos 2/3 + 



(2.13) 



and {cot ft ; tan /3 } corresponds to {u; d}-type squarks. The soft SUSY breaking term Mq, 
is given by: 



M 



D 



for u-type squarks 
for d-type squarks 



(2.14) 



In order to diagonalize the mass matrix and to determine the physical mass eigenstates 
the following rotation has to be performed: 



Si 

h 



cos 6j sin 9j 
— sin Qj cos Qj 



(2.15) 



The mixing angle Qj is given for tan ft > 1 by: 



COS#f 



sin Or 



(m 2 ~ -m 2 ) 2 
jr h 



\ m) (A f - //{cot/3 ; tan/3 }) 2 + (m 2 - R - m 2 -J 2 



(2.16) 



=F sgn 



Af — /t{cot ft ; tan /3 } 



m 2 (Af — //{cot ft ; tan ft }) s 



\ (A f - //{cot ; tan ft }) 2 + (m 2 ^ 



— m 



2 \2 



;2.i7) 



The negative sign in (|2.17|) corresponds to w-type squarks, the positive sign to d-type ones. 
m 2 ~ = M? ; + Mf cos 2ft QfS^ + m 2 denotes the lower right entry in the squark mass ma- 



trix (|2.13|) . The masses are given by the eigenvalues of the mass matrix: 



m ,- 



2 L 



Ml + M 2 Al 



+ -Mf cos 2ft il + ml 
2 ' 1 



(2.18) 



± 



± 



2L 

2 

c f 



M? - Ml, + Mf cos 2ft (l( - 2Q f s 2 jj * + Am 2 [A u - // cot ft ) 



sgn 



M? - Mj, + Mf cos 2/3 ( // - )] + Am 2 [A d - // tan ft ) ' 

M| - M|, + Mf cos 2/3 (// - 2Q/4 
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for w-type and <i-type squarks, respectively. For most of our discussions we make the choice 

M Q = M^=:m~ q = Msusy • (2.19) 

Since the non-diagonal entry of the mass matrix eq. (j2.13|) is proportional to the fermion 
mass, mixing becomes particularly important for scalar tops (/ = t), in the case of tan f3 ^> 1 
also for scalar bottoms (/ = b). 

Furthermore it is possible to express the squark mass matrix in terms of the physical 
masses m^,m^ and the mixing angle 9p 

( cos 2 6 f ~m 2 - + sin 2 fern 3 - sin 6 f ~ cos 6f(m\ - m 2 - ) \ 

Z = f fl f h o fl h ■ (2.20) 

\ sin 8? cos 8?(m 2 ~ — m 2 ~ ) sin 8fm 2 ~ + cos 2 8fm 2 ~ / 
\ J j ji fa' J n J fa / 

Af can be written as follows: 

sin Of cos 9f(m 2 f — m 2 f ) 
A f = £ ^_ + //{cot/3;tan/3}. (2.21) 

771/ 

Since the most relevant squarks for the MSSM Higgs boson sector are the t and b particles, 
here we explicitly list their mass matrices in the basis of the gauge eigenstates tjj,tji and 
b~L, b R : 

f M? L +m 2 + cos 2f3 (± - K)M 2 m t X t 
* ' \ m t X t M? r + m 2 + | cos 2/? s^Mf 

M? = j M^+m 2 + cos2/?(-| + |4)Mf m b X 6 



5 I m 6 X 6 M 2 ^ + m 2 — | cos 2/? s^Mf 




where 



m t X t = m t (A t — \i cot (3 ), m b X b = m h (A b — /itan/3 ). (2.24) 

Here A f denotes the trilinear Higgs-stop coupling, A& denotes the Higgs-sbottom coupling, 
and \i is the Higgs mixing parameter. SU(2) gauge invariance requires the relation 

M h =M bL . (2.25) 

2.3 Corrections in the Fey nman- diagrammatic approach 
2.3.1 Renormalization 

In order to calculate the higher-order corrections to the Higgs boson masses and effective 
mixing angle, the renormalized Higgs boson self-energies are needed. The parameters ap- 
pearing in the Higgs potential, see eq. (|2.1|) . are renormalized as follows: 

M 2 Z ^M 2 Z + 5M 2 Z} Th ^T h + 5T h , (2.26) 

■\/{V - .\/,Y • 5M^, t h ^T h + 6T H , 

M Hi gS s -> ^mggs + SM*^, tan/3 -> tan/3 (1 + 5tan/3 ). 
m^± — > m^± + 5m^± 
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M^ iggs denotes the tree-level Higgs boson mass matrix given in eq. (|2.8jl . T h and T H are the 
tree-level tadpoles, i.e. the terms linear in h and H in the Higgs potential. 

The field renormalization matrices of both Higgs multiplets can be set up symmetrically, 



h\^[l + \5Z hh \SZ, 
and for the charged Higgs boson 



2' 

Hj \ \5Z hH 1 + \5Z H h 




(2.27) 



E ± -> H^il + 6Z H - H+ ) . (2.28) 
For the mass counter term matrices we use the definitions 

HU *£) ' <2 ' 29) 

The renormalized self-energies, £(p 2 ), can now be expressed through the unrenormalized 
self-energies, S(p 2 ), the field renormalization constants and the mass counter terms. This 
reads for the CP-even part, 

£>hh(p 2 ) = Z hh (p 2 ) + 5Z hh (p 2 - m 2 htTee ) - 5m 2 h , (2.30a) 
t hH {p 2 ) = Z hH (p 2 ) + SZ hH (p 2 - |K )tree + m^ tree )) - 6m 2 hHl (2.30b) 
£hh(p 2 ) = ^hh(p 2 ) + SZ HH (p 2 - m^ tree ) - , (2.30c) 

and for the charged Higgs boson 

£ h -h+(p 2 ) = T, H - H +(p 2 ) + 5Z H ~ H +(p 2 - m 2 H ±) - 5m 2 H± . (2.31) 

Inserting the renormalization transformation into the Higgs mass terms leads to ex- 
pressions for their counter terms which consequently depend on the other counter terms 
introduced in ()2.2fij) . 

For the CP-even part of the Higgs sectors, these counter terms are: 

5m 2 h = 5M 2 A cos 2 (a - (3) + 5M 2 Z sin 2 (a + (3) (2.32a) 
+ 2M z l wCw ($Th cos(a - (3) sin 2 (a — 0) + 5T h sm(a — (3){l + cos 2 (a - (3))) 
+ 5tan/5 sin (3 cos (3 (M A sin 2 (a - (3) + M 2 Z sin 2 (a + (3)), 
Sm 2 hH = \{8M 2 A sm2{a- (3) - 5M| sin2(a + /?)) (2.32b) 

+ mk^( ST « sin3 (« 5T h cos 3 (a - (3)) 

- dtan/3 sin (3 cos (3 (M A cos 2(a - (3) + M 2 Z cos 2 (a + (3)), 

5m 2 H = 5M 2 A sin 2 (a - (3) + 5M 2 Z cos 2 (a + (3) (2.32c) 

- 2M ^ S c (STh cos(a — (3)(1 + sin 2 (a — (3)) + 8T h sin(a — (3) cos 2 (a — (3)) 

- 5taxif3 sin (3 cos f3 (M A sin 2(a - /3) + Mf sin 2(a + /?)) . 

For the charged Higgs boson it reads 

Sm 2 H± = 5M 2 A + . (2.33) 
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For the field renormalization we chose to give each Higgs doublet one renormalization 
constant, 

+ H 2 ^{l + \5Z H2 )H 2 . (2.34) 

This leads to the following expressions for the various field renormalization constants in 
eq. (EZZD: 

SZhh = sin 2 a SZ ni + cos 2 a SZ n2 , (2.35a) 
SZhH = sin a cos a (<5Z-h 2 — SZ^), (2.35b) 
SZ HH = cos 2 aSZ- Hl + sin 2 a SZ H2 , (2.35c) 
5Z H - H + = sin 2 /3 SZ Hl + cos 2 /3 8Z H2 . (2.35d) 

The counter term for tan (3 can be expressed in terms of the vaccuum expectation values as 

Stan/3 = ~ (5Z n , 2 - 5Z Hl ) + ^ - ^ , (2.36) 
2 v 2 vi 

where the Svi are the renormalization constants of the vf. 

v x -+ (1 + 5Z Hl ) ( Vl + M) , v 2 -> (1 + 5Z W2 ) (« 2 + <fo 2 ) . (2.37) 

The renormalization conditions are fixed by an appropriate renormalization scheme. For 
the mass counter terms on-shell conditions are used: 

5M| = ReS zz (Ml), bM 2 w = Re E W (M 2 / ), 5M 2 = Re S AA (M 2 ). (2.38) 

Here £ denotes the transverse part of the self-energy. Since the tadpole coefficients are 
chosen to vanish in all orders, their counter terms follow from T{h,H} + ST^,h} = 0: 

5T h = -T h , 5T H = -T H . (2.39) 

For the remaining renormalization constants for 5 tan (3, SZ^ and 5Zt-c 2 several choices are 
possible, see the discussion in Sect. 12.3.41 As will be shown there, the most convenient choice 
is a DR renormalization of <5tan/3, 5Z-^ 1 and SZ^, 

5tan/3 = 5t & nf3 m = -g^T [ Re S ^«trec) - ^ S'^^^J]^ , (2.40a) 
5Z Hl =5Z™ = - [Re E' HH |a=0 ] div , (2.40b) 
5Z H2 =5Z™ = - [ReS' M , Q=0 ] div . (2.40c) 

The corresponding renormalization scale, /ijjR, is set to fijy^ = m t in all numerical evalua- 
tions. 
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2.3.2 The concept of higher order corrections in the Feynman- 
diagrammatric approach 

In the Feynman diagrammatic (FD) approach the higher-order corrected CP-even Higgs 
boson masses in the rMSSM are derived by finding the poles of the (h, if )-propagator matrix. 
The inverse of this matrix is given by 



, , r l ■ I P 2 - m H,tree + ^HH(p 2 ) Z hH (p 2 , , 

(Affiggs) - -* | y ( 2 s 2 2 y ( 2) I • ^-41 J 



Determining the poles of the matrix AHi ggs in eq. (j2.41j) is equivalent to solving the equation 

P 2 - mltree + Z hh (p 2 )] [p 2 - m 2 H>tiee + ±hh( P 2 )] - [z hH ( P 2 )] * = . (2.42) 

The status of the available results for the self-energy contributions to eq. (j2.41j) can be 
summarized as follows. For the one-loop part, the complete result within the MSSM is 
known [8-11]. The by far dominant one-loop contribution is the 0(a t ) term due to top 
and stop loops (a t = h 2 /(4:7i), h t being the superpotential top coupling). Concerning the 
two-loop effects, their computation is quite advanced and has now reached a stage such that 
all the presumably dominant contributions are known. They include the strong corrections, 
usually indicated as C(a t a s ), and Yukawa corrections, 0(a 2 ), to the dominant one-loop 
0(a t ) term, as well as the strong corrections to the bottom/sbottom one-loop 0(a&) term 
(a& = hl/(4ir)), i.e. the 0(aba s ) contribution. The latter can be relevant for large values of 
tan/3 . Presently, the 0{a t a s ) [5,12,16,18,19], 0{a 2 t ) [12,20,21] and the 0{a b a s ) [22,23] con- 
tributions to the self-energies are known for vanishing external momenta. In the (s)bottom 
corrections the all-order resummation of the tan f3 -enhanced terms, 0(ab(ct s tan/3) n ), is 
also performed [43,44]. Recently the 0{a t cxh) and 0{al) corrections became available [24]. 
Most recently the C^a^av) corrections have been evaluated, which are, however, completely 
negligible [45]. Finally a "full" two-loop effective potential calculation (including even the 
momentum dependence for the leading pieces) has been published [25]. However, the latter 
results have been obtained using a certain renormalization in which all quantities, including 
SM gauge boson masses and couplings, are DR parameters. This makes them not usable in 
the approach and evaluations presented here. 

The charged Higgs boson mass is obtained by solving the equation 

p 2 -m 2 H± -t H - H+ (p 2 ) = . (2.43) 

The charged Higgs boson self-energy is known at the one-loop level [46,47]. We will not 
explore the corrections to the charged Higgs boson mass further in this report. For a detailed 
analysis, see Ref. [47]. 

2.3.3 The a e $ -approximation 

The dominant contributions for the Higgs boson self-energies can be obtained by setting 
p 2 = 0. Approximating the renormalized Higgs boson self-energies by 

±(p 2 ) -> E(0) = t (2.44) 
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yields the Higgs boson masses by re-diagonalizing the dressed mass matrix 



M 2 

1V1 Higgs 



m 



H 



'hH 



m 2 h ~ t hh 



Aa ( Mjj 



Ml 



(2.45) 



where Mh and Mh are the corresponding higher-order-corrected Higgs boson masses. The 
rotation matrix in the transformation ()2.45|) reads: 



D(Aa) 



cos Aa — sin Aa 
sin Aa cos Aa 



(2.46) 



The angle Aa is related to the renormalized self-energies and masses through the eigenvector 
equation 



m 



— E 



hH 



— Z->hH 

ml - t hh - Ml 



— sin Aa 

cos Aa 



which yields 



■'hH 



M l -™ 2 h + ^hh 
The second eigenvector equation leads to: 



tanAa . 







(2.47) 



'2.48) 



Using the relations 
and 



'hH 



M 2 H -m 2 +Z 



tanAa . 



hh 



D(a eS ) = D(a) D(Aa) 



^HH ^hH 
^hH ^hh 



D~ l (ot) 



D(a) 



(2.49) 



(2.50) 



(2.51) 



it is obvious that a c s = (a + Aa) is exactly the angle that diagonalizes the higher-order 
corrected Higgs boson mass matrix in the <pi, 02-basis: 



m^ ~ 



a 



m 



102 ^<t>i<h 
ml - XL 



"4'2 



J <(>2 



H 



'HH 



'hH 



— Z-'hH 
m 2 h - t hh 



The angle a e ff can be obtained from 

-{M 2 A + M|) sin/? cos (3 - % 



a e ff = arctan 



"off 



(j>2 



M 2 Z cos 2 (3 + Ml sin 2 (3 - E^ - m 2 h 



M% 
M 2 h 



M 2 H 
M 2 h 



7T 



7T 



— < a c ff < — 
2 2 



(2.52) 



(2.53) 
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2.3.4 Recently calculated higher-order corrections 

In order to discuss the impact of recent improvements in the MSSM Higgs sector we will make 
use of the program FeynHiggs2.1 [48,49], which is a Fortran code for the evaluation of the 
neutral CP-even Higgs sector of the MSSM including higher-order corrections to the renor- 
malized Higgs boson self-energies. The code comprises all existing higher-order corrections 
(except of the results of Ref. [25], which have been obtained in a pure DR renormalization 
scheme, see Sect. I2.3T2J) . This includes the well known full one-loop corrections [8-11], the 
two-loop leading, momentum-independent, 0{a t a s ) correction in the t/i sector [5,16,19], 
the two-loop leading logarithmic corrections at 0(af) [14,15] and all further corrections dis- 
cussed below. By two-loop momentum-independent corrections here and hereafter we mean 
the two-loop contributions to Higgs boson self-energies evaluated at zero external momenta. 
At the one-loop level, the momentum-independent contributions are the dominant part of the 
self-energy corrections, that, in principle, should be evaluated at external momenta squared 
equal to the poles of the h, if -propagator matrix, eq. (j2.41|) . 

With the implementation of the latest results obtained in the MSSM Higgs sector, Feyn- 
Higgs allows the presently most precise prediction of the masses of the CP-even Higgs bosons 
and the corresponding mixing angle. The latest version of FeynHiggs2.1 can be obtained 
from www . f eynhiggs . de. 

Hybrid renormalization scheme at the one-loop level 

FeynHiggs is based on the FD approach with on-shell renormalization conditions [5]. This 
means in particular that all the masses in the FD result are the physical ones, i.e. they 
correspond to physical observables. Since eq. (|2.42j) is solved iteratively, the result for 
and van contains a dependence on the field-renormalization constants of h and H, which 
is formally of higher-order. Accordingly, there is some freedom in choosing appropriate 
renormalization conditions for fixing the field-renormalization constants (this can also be 
interpreted as affecting the renormalization of tan/?). Different renormalization conditions 
have been considered in the literature, e.g. (E' denotes the derivative with respect to the 
external momentum squared): 

1. on-shell renormalization for E^, E^, E^, T^az, and Svx/vi = bv-ifvi [11] 

2. on-shell renormalization for E^, E^, ^az, and Svi = Svi^iv, i — 1,2 [10] 

3. on-shell renormalization for E^, E^ , DR renormalization (employing dimensional re- 
duction [50]) for 5Z Hl , 5Z n2 and tan/5 [49], see also the discussion in Sect. I2.3HH 

The original full one-loop evaluations of the Higgs boson self-energies [10,11] were based on 
type-1 renormalization conditions, thus requiring the derivative of the A boson self-energy. 
In Ref. [49] a hybrid DR/on-shell scheme, type 3, has been proposed. The choice of a DR 
definition for 5Zh, 5Zn and tan f3 requires to specify a renormalization scale Q 2 at which these 
parameters are defined, which is commonly chosen to be mj. Variation of this scale gives 
some indication of the size of unknown higher-order corrections, see Sect. 12.5.31 These new 
renormalization conditions lead to a more stable behavior around thresholds, e.g. Ma = 2 mt, 
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and avoid unphysically large contributions in certain regions of the MSSM parameter space 1 . 
This effect is demonstrated in Fig. 12.11 




Figure 2.1: The lightest Higgs boson mass is shown as a function of Ma in the m™ ax and the 
no- mixing scenario [53], see the Appendix, tan/3 has been set to tan/3 = 3 (lower curves) 
and tan/3 = 50 (upper curves). The hybrid DR/on-shell scheme (solid lines) is compared 
with the original on-shell (type-1) renormalization (dashed). The dotted lines indicate the 
effect of the "subleading" 0(af) correction. 

Two-loop 0(af) corrections 

Recently, the two- loop 0(a%) corrections in the limit of zero external momentum became 
available, first only for the lightest eigenvalue, m/,, and in the limit Ma S> Mz [20], then 
for all the entries of the Higgs propagator matrix for arbitrary values of Ma [21]. These 
corrections were obtained in the effective-potential approach, that allows to construct the 
Higgs boson self-energies, at zero external momenta, by taking the relevant derivatives of 
the field-dependent potential. In this procedure it is important, in order to make contact 
with the physical Ma, to compute the effective potential as a function of both CP-even 
and CP-odd fields, as emphasized in Ref. [19]. In the evaluation of the 0{af) corrections, 
the specification of a renormalization prescription for the Higgs mixing parameter \i is also 
required and it has been chosen as DR. In FeynHiggs2.1, which includes the two-loop 0(pif) 
corrections, the corresponding renormalization scale is fixed to be the same as for SZ^, SZ^ 
and tan (3. 

The availability of the complete result for the momentum-independent part of the 0(af) 
corrections allows to judge the quality of results that incorporate only the logarithmic con- 
tributions [14,15]. In Fig. 12.21 fsee also Fig. 12. ip we plot the two-loop corrected m h as a 
function of the stop mixing parameter, X t . For simplicity, the soft SUSY breaking para- 
meters in the diagonal entries of the stop mass matrix, Mi L , Mi R , are chosen to be equal, 
Mf L = Mi = Msusy- For the numerical analysis Msusy as well as Ma and fi, are chosen 
to be all equal to 1 TeV, while the gluino mass is rrig = 800 GeV, and tan j3 = 3. If not 

1 A more detailed discussion can be found in Refs. [49,51]; see also Ref. [52]. 
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Figure 2.2: Two-loop corrected as a function of X t in various steps of approximation. 
The relevant MSSM parameters are chosen as tan/? = 3 , m* olc = 174.3 GeV, M h = M~ tR = 
Ma = fJ> = 1 TeV and m- g = 800 GeV. The meaning of the different curves is explained in 
the text. 

otherwise stated, in all plots below we choose the trilinear couplings in the stop and sbottom 
sectors equal to each other, A& = A t , and set M 2 = Msusy, where M 2 is the SU(2) gaug- 
ino mass parameter (the U(l) gaugino mass parameter is obtained via the GUT relation, 
Mi = 5/3 4/4 M 2 ). 

The solid and dotted lines in Fig. 12. 21 are computed with and without the inclusion of the 
full corrections, while the dashed and dot-dashed ones are obtained including only the 

logarithmic contributions. In particular, the dashed curve corresponds to the result of [14], 
obtained with the one-loop renormalization group method. The dot-dashed curve, instead, 
corresponds to the leading and next-to-leading logarithmic terms. For a more detailed dis- 
cussion see Ref. [7]. From Fig. 12.21 it can also be seen that, for small X t , the full 0(af) 
result is very well reproduced by the logarithmic approximation, once the next-to-leading 
terms are correctly taken into account. On the other hand, when X t is large there are sig- 
nificant differences, amounting to several GeV, between the logarithmic approximation and 
the full result. Such differences are due to non-logarithmic terms that scale like powers of 
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X t /M S u SY . It should be noted that for more general choices of the MSSM parameters the 
renormalization group method becomes rather involved (see e.g. Ref. [54], where the case of a 
large splitting between M^ L and is discussed), and a suitable next-to- leading logarithmic 
approximation to the full result is much more difficult to devise. 



Two-loop sbottom corrections 

Due to the smallness of the bottom mass, the 0(ab) one-loop corrections to the Higgs boson 
self-energies can be numerically non-negligible only for tan (3 >> 1 and sizable values of the 
\x parameter. In fact, at the classical level hb/h t = {m b /m t ) tan/5, thus tan/3 ^> 1 is needed 
in order to have a b ~ a t in spite of m b <C m t . In contrast to the 0(a t ) corrections where 
both top and stop loops give sizable contributions, in the case of the 0(ab) corrections the 
numerically dominant contributions come from sbottom loops: those coming from bottom 
loops are always suppressed by the small value of the bottom mass. A sizable value of ji is 
then required to have sizable sbottom-Higgs scalar interactions in the large-tan (3 limit. 

The relation between the bottom-quark mass and the Yukawa coupling hb, which con- 
trols also the interaction between the Higgs fields and the sbottom squarks, reads at lowest 
order nib = hbVi/y/2. This relation is affected at one-loop order by large radiative correc- 
tions [43,44] (see also Ref. [55]), proportional to h b v 2 , in general giving rise to tan /3-enhanced 
contributions. These terms proportional to v 2 , often indicated as threshold corrections to 
the bottom mass, are generated either by gluino-sbottom one-loop diagrams, resulting in 
0(ab<y s ) corrections to the Higgs masses, or by chargino-stop loops, giving 0(ab<y t ) correc- 
tions. Because the tan /3-enhanced contributions can be numerically relevant, an accurate 
determination of hb from the experimental value of the bottom mass requires a resummation 
of such effects to all orders in the perturbative expansion, as described in Ref. [44]. 

The leading effects are included in the effective Lagrangian formalism developed in 
Ref. [44] . Numerically this is by far the dominant part of the contributions from the sbottom 
sector. The effective Lagrangian is given by 



£ = 9 rn b 



2M W 1 + Am b 



tan/5 Ai 6756 + V2 V tb tan f3 H + t L b 



R 



. sin a cos a , 

+ 7T - Am;,— hb L b R 

cos 3 sin p 



cos a . sin a: , 

+ Am b —- Hb L b R 



cos (3 sin (3 



+ h.c. . (2.54) 



Here nib denotes the running bottom quark mass including SM QCD corrections. In the 
numerical evaluations we choose nib = nTb(m t ) ~ 2.74 GeV. Furthermore, Anib is given at 
0(a s ) by 

Am b = —^rcigH tan (3 x I \m~ h m~ h nig) , (2.55) 

O 7T 

where I is given by 
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The large b — g loops are resummed to all orders of (a s tan f3) n via the inclusion of Amj [44]. 
The prefactor 1 / (1+ Am&) in eq. (J2.54)) arises from the resummation of the leading corrections 
to all orders. They are due to the shift in the relation of the bottom Yukawa coupling to 
the physical b quark mass, coming from the loop induced coupling of the </>i to the bottom 
quarks. The additional terms ~ Am;, in the hbb and Hbb couplings arise from the mixing 
and coupling of the "other" Higgs boson, H and h, respectively, to the b quarks. 

The b/b corrections to the renormalized Higgs boson self-energies in eq. (|2.41|) are given 
at the one-loop level by the Feynman diagrams shown in Figs. 12.31 12.41 

b bj 
H / \ H H / \ H 



b; 




bj d j 



H \ / H 



Figure 2.3: Generic Feynman diagrams for the b/b contributions to Higgs boson self-energies 
(H = h, H, A). 



H ' * u H 

i y D; 




Figure 2.4: Generic Feynman diagrams for the b/b contributions to Higgs tadpoles (H = 
h,H). 

The Ami, corrections included in the effective Lagrangian, eq. (|2.54|) . enter the one- loop 
diagrams in the following ways: 

• The factor (m& sin a / cos (3 ) in the hbb coupling changes to 

sin a m& /sin a . cos a \ , n „„. 

m b 77 -> ; r — 77 ~ Am fe — - . 2.57 

cos p 1 + Am b \ cos p sm p J 

The corresponding factor (— mj, cos a / cos (3 ) in the Hbb coupling is shifted to 

cos a rrih ( cos a . sin a \ , „„. 

-m b £ — - + Am b -J. (2.58) 

cos p 1 + Arrib \ cos p sin p J 

Finally, the Abb coupling receives the factor 

^ TTa^' (2 ' 59) 

17 



It should be noted that the b quark masses appearing in the bottom propagators do 
not receive a correction. 



• The b quark masses in the b mass matrix ()2.23|) are due to Yukawa couplings from the 
scalar bottoms to the Higgs doublets. Thus they receive the shift 

m b — ► . (2.60) 

• For the same reason the m b factors in the Higgs-66 couplings receive the same shift as 
given in eq. (|2.6U|) . 

Using the above changes in the one-loop corrections includes the leading b/b corrections 
beyond one-loop, resummed to all orders. 

Recently the complete two-loop momentum-independent, 0(a b a s ) corrections (which are 
not included in the 0(a b (a s tan/3) n ) resummation) have been computed [22,23]. The result 
obtained makes use of an appropriate choice of renormalization conditions on the relevant 
parameters that allows to disentangle the genuine two-loop effects from the large threshold 
corrections to the bottom mass, and also ensures the decoupling of heavy gluinos. 

To appreciate the importance of the various sbottom contributions, we plot in Fig. 12.51 
the light Higgs mass Mh as a function of tan/3. The SM running bottom mass computed 
at the top mass scale, m b (mt) = 2.74 GeV, is used in order to account for the universal 
large QCD corrections. The relevant MSSM parameters are chosen as Ma = 120 GeV, 
fi = -1 TeV, M h = M~ tR = M- bR = m~ g = 1 TeV, A t = A b = 2 TeV. The dot-dashed 
curve in Fig. 12.51 includes the full one-loop contribution as well as the two-loop 0(a t a s + a 2 ) 
corrections (the latter being approximately tan /^-independent when tan/5 is large). The 
dashed curve includes also the resummation of the tan /3-enhanced threshold effects in the 
relation between h b and m b . Finally, the solid curve includes in addition the complete 
0(a b a s ) two-loop corrections of Ref. [22]. In the last two curves, the steep dependence of 
nih on tan (3 when the latter is large is driven by the sbottom contributions. One sees that the 
tan /3-enhanced threshold effects account for the bulk of the sbottom contributions beyond 
one-loop. The genuine 0(a b a s ) two-loop corrections can still shift rrih by a few GeV for very 
large values of tan (3 and \i. 



2.3.5 Implications for tan/3 exclusion bounds 

The improved knowledge of the two-loop contributions to the Higgs boson self-energies results 
in a very precise prediction for the Higgs boson masses and mixing angle with interesting 
implications for MSSM parameter space analyses. In this section an important consequence 
of the various corrections is presented: the implications on the upper limit on rrih within the 
MSSM and on the corresponding limit on tan (3 arising from confronting the upper bound 
on rrih with the lower limit from Higgs searches, see also the discussion in Ref. [6]. 

The theoretical upper bound on the lightest Higgs boson mass as a function of tan/3 
can be combined with the results from direct searches at LEP to constrain tan/3. The 
diagonalization of the tree- level mass matrix, eq. ([2.8)1 . yields a value for m hjtree that is 
maximal when Ma 3> Mz, in which case m\ tree ~ M| cos 2 2(3 , which vanishes for tan (3 = 1. 
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Figure 2.5: The result for the lightest CP-even Higgs-boson mass in the MSSM, M h , as 
obtained with the program FeynHiggs is shown as a function of tan f3 for Ma = 120 GeV, 
pi = -1 TeV, M ih = M- tR = M 1r = m-g = 1 TeV, A t = A b = 2 TeV. The meaning of the 
different curves is explained in the text. 

Radiative corrections significantly increase the light Higgs boson mass compared to its tree- 
level value, but still m h is minimized for values of tan/3 around one. Thus, in principle, 
the region of low tan (3 can be probed experimentally via the search for the lightest MSSM 
Higgs boson [56]. If the remaining MSSM parameters are tuned in such a way to obtain the 
maximal value of Mh as a function of tan (3 (for reasonable values of MgusY and taking into 
account the experimental uncertainties of m t and the other SM input parameters as well as 
the theoretical uncertainties from unknown higher-order corrections), the experimental lower 
bound on M h can be used to obtain exclusion limits for tan (3. While in general a detailed 
investigation of a variety of different possible production and decay modes is necessary in 
order to determine whether a particular point of the MSSM parameter space can be excluded 
via the Higgs searches or not, the situation simplifies considerably in the region of small tan f3 
values. In this parameter region the lightest CP-even Higgs boson of the MSSM couples to 
the Z boson with SM-like strength, and its decay into a bb pair is not significantly suppressed. 
Thus, within good approximation, constraints on tan/3 can be obtained in this parameter 
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region by confronting the exclusion bound on the SM Higgs boson with the upper limit on 
Mh within the MSSM. We use this approach below in order to discuss the implications of 
the new Mh evaluation on tan /3 exclusion bounds. 

Concerning the upper bound on Mh within the MSSM, the one-loop corrections contribute 
positively to M%. The two-loop effects of 0(a t a s ) and 0(a^), on the other hand, enter 
with competing signs, the former reducing M% while the latter give a (smaller) positive 
contribution. The actual bound that can be derived depends sensitively on the precise value 
of the top-quark mass, because the dominant one- loop contribution to M^, as well as the 
two-loop 0(a t a s ) term, scale as m\. Furthermore, a large top mass amplifies the relative 
importance of the two-loop Oiptf.) correction, because of the additional ml factor. 

In order to discuss restrictions on the MSSM parameter space it has become customary in 
the recent years to refer to so-called benchmark scenarios of MSSM parameters [53,57]. The 
m max b e ri C fiHiark scenario [53] has been designed such that for fixed values of m t and Msusy 
the predicted value of the lightest CP-even Higgs boson mass is maximized for each value of 
Ma and tan (3. The value of the top-quark mass had been fixed to its experimental central 
value, m t = 174.3 GeV, while the SUSY parameters are taken as in the m™ ax scenario, see 
the Appendix. Furthermore, Ma has been fixed to Ma = 1 TeV. 

In Fig. 12.61 we plot Mh as a function of tan (3 in the m™ ax scenario. The dashed and 
dot-dashed curves correspond to the result obtained with the previous (used for the LEP 
evaluations so far [56]) and the more advanced version (where the improvements described in 
Sect. 12.3.41 are included, and which will be used for the final LEP evaluations [58]) of Feyn- 
Higgs, respectively. The two versions, FeynHiggsl.O and FeynHiggsl.3, differ by the recent 
improvements obtained in the MSSM Higgs sector which are described in Sect. 12.3.41 Also 
shown as a solid line is the result obtained with FeynHiggs2.l (which for the m™ ax scenario 
should yield quantitatively the same result as FeynHiggsl.3), but with mt set to the recently 
obtained new experimental value, m^ xp = 178.0 GeV [3]. For comparison, also the result ob- 
tained with a renormalization-group improved effective potential method is indicated. The 
dotted curve in Fig. 12.61 corresponds to the code subhpoledm [14,44,59] in the m™ ax scenario, 
for a t mixing parameter X t MS = v6 Msusy [14,59] (see also the Appendix). It deviates 
from the result of FeynHiggsl.O by typically not more than 1 GeV for tan/3 > 1. The LEP 
exclusion bound for the mass of a SM-like Higgs [60], Mh sm > 114.4 GeV, is shown in the 
figure as a vertical long-dashed line. As can be seen from the figure, the improvements on 
the theoretical prediction described in Sect. 12.3.41 in particular the inclusion of the complete 
momentum-independent 0(af) corrections into FeynHiggs, gives rise to a significant increase 
in the upper bound on Mh as a function of tan (3. Comparison of this prediction with the 
exclusion bound on a SM-like Higgs shows that the lower limit on tan/3 is considerably 
weakened. Also the new top mass value, see the solid line, shifting the resulting M^ values 
upwards, weakens the tan f3 exclusion bound. 

Concerning the interpretation of the results shown in Fig. 12.61 ft should be kept in 
mind that within the (pure) m™ ax benchmark scenario m< and Msusy are kept fixed, and 
no theoretical uncertainties from unknown higher-order corrections are taken into account. 
In order to arrive at a more general exclusion bound on tan (3 that is not restricted to a 
particular benchmark scenario, the impact of the parametric and higher-order uncertainties 
in the prediction for Mh has to be considered [6]. In order to demonstrate in particular the 
dependence of the tan (3 exclusion bound on the chosen value of the top pole mass, besides the 
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Figure 2.6: The result for the lightest CP-even Higgs boson mass, M^, as a function of tan/? 
in the m™ ax scenario. The dotted curve has been obtained with a renormalization-group 
improved effective potential method, using the code subhpoledm. The short-dashed curve 
corresponds to the result obtained with FeynHiggsl.O, while the dash-dotted curve shows the 
result of FeynHiggsl.3, where the improvements described in Sect. I2.3T41 are included. The 
full curve shows the result from FeynHiggs2.1, using the top quark mass of m t = 178.0 GeV. 
The dot-dashed curve, also obtained with FeynHiggs2.1, employs m t = 182.3 GeV and 
^susy = 2 TeV. The vertical long-dashed line corresponds to the LEP exclusion bound 
for the SM Higgs boson of 114.4 GeV. 



solid line obtained with m t = 178.0 GeV, the dot-dashed curve in Fig. 12.61 shows the result 
obtained with FeynHiggs2.l where the top-quark mass has been increased by one standard 
deviation, a mt = 4.3 GeV, to m t = 182.3 GeV, and Msusy has been changed from 1 TeV 
to 2 TeV. It can be seen that in this more general scenario no lower limit on tan/3 from the 
LEP Higgs searches can be obtained. 

Constraints from the Higgs searches at LEP do of course play an important role in regions 
of the MSSM parameter space where the parameters are such that Mh does not reach its 
maximum value. Also in this case, however, the remaining theoretical uncertainties from 
unknown higher-order corrections (see Sect. 12.51 below) have to be taken into account in 
order to obtain conservative exclusion limits. 
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2.4 Comparison with the RG approach 



The diagrammatic two-loop computation of the dominant contributions at 0(aa s ) to the 
neutral CP-even Higgs boson masses [16] had been obtained first in the on-shell scheme, 
was subsequently combined [5] with the complete diagrammatic one-loop on-shell result of 
Ref. [11]. Within the RG approach [14,15], on the other hand, the two-loop results was 
expressed in terms of the top-quark mass in the MS scheme. 

Besides a sizable shift in the upper bound of Mh of about 5 GeV, apparent deviations be- 
tween the explicit diagrammatic two-loop calculation and the results of the RG computation 
were observed in the dependence of Mh on the stop- mixing parameter X t . While the value of 
X t that maximizes the lightest CP-even Higgs mass is (X t ) max ~ ±\^6Ms ~ ±2.4M S (Ms 
denotes the soft SUSY-breaking parameter in the t and b sector, Ms = Mf L = Mi R = = 
M~ bR ) in the RG results, the corresponding on-shell two- loop diagrammatic computation 
found a maximal value for M h at (X t ) max ks 2M s . 2 Moreover, the RG result is symmetric 
under X t — > — X t and has a (local) minimum at X t = 0. In contrast, the two-loop diagram- 
matic computation yields M h values for positive and negative X t that differ significantly 
from each other and the local minimum in Mh is shifted slightly away from X t = [17]. 

In this section, it is shown that this apparent discrepancy is caused by the different re- 
normalization schemes employed in the two approaches, leading to differences in the leading- 
logarithmic contributions. It will furthermore be shown that in the analytic approximation 
for the leading m\ corrections, the dominant numerical contribution of the new genuine 
non- logarithmic two-loop contributions of the FD result [5, 16] can be absorbed into an ef- 
fective one-loop expression by choosing an appropriate scale for the running top-quark mass 
in different terms of the expression. 



2.4.1 Approximation formulas for 



Within the RG approach the following approximation for Mh has been obtained [13-15,59], 
taking into account terms up to 0(aa s ): 
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, (2.61) 



where we have introduced the notation Ms,X t to emphasize that the corresponding quan- 
tities are MS parameters, which are evaluated at the scale \i = M s : 



M s = M| lb (M 5 ), X t = X t Mb (M 5 ), 
and m t = rnf^ M (m t ) is the MS top mass 

M t 



m t = m t (m t ) 



1 + ^a s (m t ) 



(2.62) 



(2.63) 



2 A local maximum for Mh is also found for X t ~ ~2M$, although the corresponding value of Mh at 
X t ~ +2Ms is significantly larger [5, 17]. 



22 



M t denotes the top quark pole mass. Since in this section only the 0(aa s ) corrections are 
investigated, we do not include here the contributions of 0(a t ), see however Sect. 12.5141 

Concerning the FD result, we consider the dominant one-loop and two- loop terms and we 
focus on the case Ma 3> Mz, for which the result for M| can be expressed in a particularly 
compact form [17] 

Ml = m^ + ml a + ml aa % (2.64) 

and neglect the non-leading terms of 0{M\ /M|). Assuming that Ms 3> M t and neglecting 
the non-leading terms of 0(M t /Ms) and 0(M z /M^), one obtains the following simple result 
for the one-loop and two-loop contributions 

mh ~ 2 ^ Mt \ [U {mi) + MI\ L 12 Ml) J' {I - bb) 

(2.66) 

The corresponding formulae, in which terms up to 0(M t /M|) are kept, can be found in 
Ref. [59]. 

In eqs. fl2.65j) - ()2.66j) the parameters M t , Ms, X t are on-shell quantities. Using eq. ()2.63j) . 
the on-shell result for M%, eqs. (j2~Mjl -(|27HH |l . can easily be rewritten in terms of the running 
top-quark mass fn t . While this reparameterization does not change the form of the one-loop 
result, it induces an extra contribution at 0(aa s ). Keeping again only terms that are not 
suppressed by powers of m t /Ms, the resulting expressions read 



3G F y/2_ 



G F V2 



m h = _3 5 m t 

11 7T 2 7T 




+ t- 2 W s -3W s + ^Ml^ (2 ' 68) 

in accordance with the formulae given in Ref. [17]. 

We now compare the diagrammatic result expressed in terms of the parameters m t , Ms, 
X t , eqs. (I2.67|) - ()2.68J) . with the RG result of eq. ()2.61|) . which is given in terms of the MS 
parameters m t , Ms, X t eq. ()2.62)) . While the X t -independent logarithmic terms are the 
same in both the diagrammatic and RG results, the corresponding logarithmic terms at two 
loops that are proportional to powers of X t and X t , respectively, are different. Furthermore, 
eq. (|2.68j) does not contain a logarithmic term proportional to Xf, while the corresponding 
term proportional to X t appears in eq. ()2.61|) . To check whether these results are consistent 
with each other, one must relate the on-shell and MS definitions of the parameters Ms and 
X t . 

Finally, we note that the non-logarithmic terms contained in eq. ()2.68|) correspond to 
genuine two-loop contributions that are not present in the RG result of eq. (|2.61|) . They 
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can be interpreted as a two-loop finite threshold correction to the quartic Higgs self-coupling 
in the RG approach. In particular, note that eq. (|2.68|) contains a term that is linear in 
X t . This is the main source of the asymmetry in the two-loop corrected Higgs mass under 
X t — > —X t obtained by the diagrammatic method. The non-logarithmic terms in eq. f!2.68|) 
give rise to a numerically significant increase of the maximal value of Mh of about 5 GeV in 
this approximation. 



2.4.2 On-shell and MS definitions of M s and X 



t 

Since the parameters p = {m~^rri~ ,df,mt} of the t—t sector are renormahzed differently 
in different schemes, the parameters Ms and X t also have a different meaning in these 
schemes. The relation between these parameters in the MS and in the on-shell scheme have 
been derived in Refs. [59,61] and read in leading order in mt/Mg: 

M 2 S = M 2 s os --^M 2 Sl (2.69) 

3 7T 

X t = ^ OS =#rT + !- M - (2-70) 
m t (M s ) 3 vr 

As previously noted, it is not necessary to specify the definition of the parameters that 
appear in the 0(a s ) terms as long as higher orders are neglected. Thus, we use the generic 
symbol Mj in the 0(a s ) terms of eqs. ()2.69j) - (j2.7()j) . The corresponding results including 
terms up to O (mf/Mg) can be found in Ref. [59]. 

Finally, the evaluation of the ratio M t /fn t (Ms) is needed. In leading order in m t /Ms it 
is given by [59]: 



m t {M s ) = m t 



a s . ( mf\ a s X t 



1+ vr ln \M$) + 3ttM s 



[2.71) 



where m t = rn^ M {m t ) is given in terms of M t by eq. ()2.63|) . Note that the term in eq. (j2.71|) 
that is proportional to X t is a threshold correction due to the supersymmetry-breaking 
stop-mixing effect. Inserting the result of eq. (|2.71j) into eq. (|2.7(J|) yields 



AX t X? 3X t , fmf 
H — r~~ — r-pr — t-~t~ m ' 



M s M| M s V M i 



(2.72) 



It is interesting to note that X t ^ when X® s = 0. Moreover, it is clear from eq. ()2.72j) 
that the relation between X t defined in the on-shell and the MS schemes includes a lead- 
ing logarithmic effect, which has to be taken into account in a comparison of the leading 
logarithmic contributions in the RG and the two-loop diagrammatic results. 

A remark on the regularization scheme is in order here. In the effective field theory, 
the running top-quark mass at scales below Ms is the SM running coupling of eq. (|2.63j) . 
which is calculated in dimensional regularization. This is matched to the running top-quark 
mass as computed in the full supersymmetric theory. One could argue that the appropriate 
regularization scheme for the latter should be dimensional reduction (DRED) [50], which is 
usually applied in loop calculations in supersymmetry. 3 The result of such a change would 



3 In order to obtain the corresponding DRED result, one simply has to replace the term 4a s /37r in the 
denominator of eq. (|2.ti3|) by 5a s /37r. 
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be to modify slightly the two-loop non-logarithmic contribution to Mh that is proportional 
to powers of X t . Of course, the physical Higgs mass is independent of any scheme. One is 
free to re-express eqs. (|2.65|) - (j2.66|) (which depend on the on-shell parameters M t , M s , X t ) in 
terms of parameters defined in any other scheme. In this paper, we find MS-renormalization 
via DREG to be the most convenient scheme for the comparison of the diagrammatic and 
RG results for M^. 



2.4.3 Comparing the RG and the Feynman diagrammatic results 

In order to directly compare the two-loop diagrammatic and RG results, we must convert 
from on-shell to MS parameters. Inserting eqs. ()2.69|) - ()2.72|) into eqs. ()2.67|) - (J2.68J) . one finds 



3G F V2_ 




< a = ( =2 ) + =2(1- -=0 } } ■ (2.73) 

2,aa s Q G F V2a s 4 

m h = - 3 2 m t 

X t [2 7X? IXf 1X, 4 \1 ,„{m t \ 

+ =M =^ + — =^ + — =k)}+0 ( =£- J. 2.74 

M s V3 9M| 36M| 18M|yJ \M S J V ; 

Comparing eq. (|2.74j) with eq. f|2 . 6 lj) shows that the logarithmic contributions of the 
diagrammatic result expressed in terms of the MS parameters m t , Ms, X t agree with the 
logarithmic contributions obtained by the RG approach. The differences in the logarithmic 
terms observed in the comparison of eqs. (|2.67j) - (|2.68|) with eq. (|2.61|) have thus been traced 
to the different renormalization schemes applied in the respective calculations. The fact 
that the logarithmic contributions obtained within the two approaches agree after a proper 
rewriting of the parameters of the stop sector is an important consistency check of the 
calculations. In addition to the logarithmic contributions, eq. f!2.74|) also contains non- 
logarithmic contributions, which are numerically sizable. 

In Fig. 12.71 we compare the diagrammatic result for Mh in the leading m\ approximation 
to the results obtained by RG techniques. While the diagrammatic result expressed in terms 
of m t , Mg, X t agrees well with the RG result in the region of no mixing in the stop sector, 
sizable deviations occur for large mixing. In particular, the non-logarithmic contributions 
give rise to an asymmetry under the change of sign of the parameter X t , while the RG 
result is symmetric under X t — > —Xp In the approximation considered here, the maximal 
value for Mh in the diagrammatic result lies about 3 GeV higher than the maximal value of 
the RG result. The differences are slightly larger for smaller tan/? values. In addition, as 
previously noted, the maximal-mixing point (X 4 ) max (where the radiatively corrected value 
of M h is maximal) is equal to its one-loop value, (X t ) max ~ ±^/6M s , in the RG result 
of eq. (J2.61)) . while it is shifted in the two-loop diagrammatic result. However, Fig. 12.71 
illustrates that the shift in (X t ) max from its one- loop value, while significant in the two- loop 
on-shell diagrammatic result, is largely diminished when the latter is re-expressed in terms 
of MS parameters. 

The differences between the diagrammatic and RG results shown in Fig. 12.71 can be at- 
tributed to non- negligible non- logarithmic terms proportional to powers of X t . Clearly, the 
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Figure 2.7: Comparison of the diagrammatic two-loop 0(aa s ) result for Mh, to leading order 
in m 4 /M 5 (eqs. $T7ty - $T7ty ) with the RG result of eq. (J22IJ). Note that the latter omits 
the one-loop threshold corrections due to stop mixing in the evaluation of m t (Ms). Since 
this quantity enters in the definition of X t , see eq. (|2.70p . the meaning of X t plotted along 
the x-axis is slightly different for the diagrammatic curve, where X t MS = X t , and the RG 
curve, where X t MS = X t [1 + (a s /37r)(X t / Ms)). See text for further details. The graph is 
plotted for M S = M A = (mf + m 2 ) 1 ' 2 = 1 TeV. 

RG technique can be improved to incorporate these terms. In Ref. [15], it was shown that 
the leading two-loop contributions to M\ given by the RG result of eq. (j2.61|) could be 
absorbed into an effective one-loop expression. This was accomplished by considering sepa- 
rately the X t -independent leading double logarithmic term (the "no-mixing" contribution) 
and the leading single logarithmic term that is proportional to powers of X t (the "mixing" 
contribution) at 0(aa s ). Both terms can be reproduced by an effective one-loop expression, 
where m t in eq. ()2.73|) . which appears in the no- mixing and mixing contributions, is replaced 
by the running top-quark mass evaluated at the scales \i t and /i^, respectively: 



That is, at 0(aa s ), the leading double logarithmic term is precisely reproduced by the single- 
logarithmic term at 0(a), by replacing m t with rn t (fit), while the leading single logarithmic 
term at two loops proportional to powers of X t is precisely reproduced by the corresponding 
non-logarithmic terms proportional to X t by replacing m t with m t {Ms). 

Applying the same procedure to eq. ()2.74|) and rewriting it in terms of the running top- 



no mixing: fi t = (m t Ms) 



mixing: = Ms ■ 



(2.75) 



26 



Leading m, 4 approximation at O(oc oc s ) 



O 120 



I I I I I I I I I I I I I I I I I I 

diagrammatic 

"mixed scale" one-loop RG 


I 1 


1 1 

\ 
\ 
\ - 






\ - 


-/'/ 

V 




_ 


M S MS = 1 000 GeV, M A = 1 000 GeV, tanp 


= 30 

i 


1 


-3000 -2000 -1000 1000 

X t MS [GeV] 


2000 


3000 



Figure 2.8: Comparison of the diagrammatic two-loop 0(aa s ) result for Mh, to leading order 
in m t /M s (eqs. (j2~73jl - (j!T73Jl ) with the "mixed-scale" one-loop RG result of eq. (J2.76)) . Note 
that the latter now includes the threshold corrections due to stop mixing in the evaluation of 
m t (Ms) in contrast to the RG results depicted in Fig. 12.71 "Mixed-scale" indicates that in 
the no-mixing and mixing contributions to the one-loop Higgs mass, the running top quark 
mass is evaluated at different scales according to eq. (|2.75p . See text for further details. The 
graph is plotted for M s = M A = (m| + m t 2 ) 1/2 = 1 TeV. 

quark mass at the corresponding scales as specified in eq. (|2.75j) . we obtain 

- -3^^m 4 |iln^+^^-I^L + l^H (2 77) 

Indeed, the Xf-independent leading double logarithmic term and the leading single logarith- 
mic term that is proportional to powers of X t have disappeared from the two-loop expression 
eq. ()2.77|h having been absorbed into an effective one- loop result (I2.76J) . (denoted henceforth 
as the "mixed-scale" one- loop RG result). Among the terms that remain in eq. (|2.77|) . there 
is a subleading one-loop logarithm at two loops which is a remnant of the no-mixing con- 
tribution. But, note that the magnitude of the coefficient (1/6) has been reduced from 
the corresponding coefficients that appear in eqs. (|2.66|) - (j2.74j) (—2 and 2/3, respectively). 
In addition, the remaining leftover two-loop non- logarithmic terms are also numerically in- 
significant. We conclude that the "mixed-scale" one- loop RG result provides a very good 
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approximation to the leading m\ corrections to M\ at 0(a s a t ), in which the most significant 
two-loop terms have been absorbed into an effective one-loop expression. 

To illustrate this result, we compare in Fig. 12.81 the diagrammatic two-loop result ex- 
pressed in terms of MS parameters (eqs. ()2.73j) - (j2.74j0 with the "mixed-scale" one- loop RG 
result of eq. ()2.76|) as a function of X t . The difference between the solid and dashed lines of 
Fig. I2.8l is precisely equal to the leftover two- loop term given by eq. (j2.77J) . which is seen to be 
numerically small. Hence, within the simplifying framework under consideration (i.e., only 
leading t-i sector-contributions are taken into account assuming a simplified stop squared- 
mass matrix ()2.22|) . with Mg, Ma 3> fn t and m g = nig), we see that the "mixed-scale" 
one-loop result for M h provides a very good approximation to a more complete two-loop 
result for all values of X t . 4 

2.5 Missing higher-order corrections and parametric 
uncertainties 

The prediction for Mh in the MSSM is affected by two kinds of uncertainties, parametric 
uncertainties from the experimental errors of the input parameters and uncertainties from 
unknown higher-order corrections. 

2.5.1 Parametric uncertainties 

Currently the parametric uncertainties dominate over those from unknown higher-order cor- 
rections, as the present experimental error of the top-quark mass of about ±4 GeV [3] 
induces an uncertainty of AMh ~ ±4 GeV [6]. However, at the next generation of colliders 
m t will be measured with a much higher precision, reaching the level of about 0.1 GeV at 
an e + e~ LC [37]. Thus, the m r induced parametric error will be drastically reduced, see also 



Besides nit, the other SM input parameters whose experimental errors can be relevant 
for the prediction of Mh are Mw, ct s , and nib- The W boson mass Mw mainly enters via 
the reparameterization of the electromagnetic coupling a(0) in terms of the Fermi constant 



where the quantity Ar summarizes the radiative corrections. 

The present experimental error of the W boson mass of 34 MeV leads to a parametric 
theoretical uncertainty of Mh below 0.1 GeV. In view of the prospective improvements 
in the experimental accuracy of Mw the parametric uncertainty induced by Mw will be 
substantially smaller than the one induced by nit even for Sm t = 0.1 GeV. 

The current experimental error of the strong coupling constant, 5a s (Mz) = ±0.002 [64], 
induces a parametric theoretical uncertainty of Mh of about 0.3 GeV. Since a future im- 
provement of the error of a s (M z ) by about a factor of two can be envisaged [65,66], the 

4 Strictly speaking, the analytic approximations discussed here break down when rfitXt ~ M s . Thus, 
one does not expect an accurate result for the corresponding formulae when X t is too large [13-15,17]. In 
practice, one should not trust the accuracy of the analytic formulae once X t > (^t)max- 



Refs. [62,63]. 



Gp 




(2.78) 
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parametric uncertainty induced by m t will dominate over the one induced by a s [Mz) down 
to the level of 5m t = 0.1-0.2 GeV. The effect of the experimental uncertainties in mj are 
negligible, once the proper resummation of the leading effects [43,44] is taken into account. 

Also the (so far unknown) SUSY parameters have a large impact on the value of Mh- 
Most prominent here are the parameters of the t sector. However, the induced uncertainty 
will depend strongly on the future experimental uncertainty, and can thus currently not very 
well estimated. Therefore these uncertainties are not further investigated here. 

2.5.2 Estimating missing two-loop corrections 

Given our present knowledge of the two-loop contributions to the Higgs boson self-energies, 
see Sect. 12.3.41 the theoretical accuracy reached in the prediction for the CP-even Higgs boson 
masses is quite advanced. However, obtaining a complete two-loop result for the Higgs boson 
masses and mixing angle requires additional contributions that are not yet available 5 . In this 
and the following sections we discuss the possible effect of the missing two-loop corrections, 
and we estimate the size of the higher-order (i.e. three-loop) contributions. 

It is customary to separate the corrections to the Higgs boson self-energies into two parts: 

i) the momentum- independent part, namely the contributions to the self-energies evaluated 
at zero external momenta, which can also be computed in the effective potential approach; 

ii) the momentum-dependent corrections, i.e. the effects induced by the dependence on the 
external momenta of the self-energies that are required to determine the poles of the (h, H)- 
propagator matrix. 

All the presently available two-loop contributions are computed at zero external momen- 
tum, and moreover they are obtained in the so-called gaugeless limit, namely by switching 
off the electroweak gauge interactions (with the already discussed exception in Ref. [25], see 
Sect. I2.3.2J) . The two approximations are in fact related, since the leading Yukawa corrections 
are obtained by neglecting both the momentum dependence and the gauge interactions. In 
order to systematically improve the result beyond the approximation of the leading Yukawa 
terms, both effects from the gauge interactions and the momentum dependence should be 
taken into account. 

To try to estimate, although in a very rough way, the importance of the various contri- 
butions we look at their relative size in the one- loop part. There, in the effective potential 
part, the effect of the 0(a t ) corrections typically amounts to an increase in of 40-60 
GeV, depending on the choice of the MSSM parameters, whereas the corrections due to the 
electroweak (D-term) Higgs-squark interactions usually decrease M h by less than 5 GeV [9] . 
Instead, the purely electroweak gauge corrections to Mh, namely those coming from Higgs, 
gauge boson and chargino or neutralino loops [11], are typically quite small at one- loop and 
can reach at most 5 GeV in specific regions of the parameter space (namely for large values of 
\i and M 2 ). Concerning the effects induced by the dependence on the external momentum, as 
a general rule we expect them to be more relevant in the determination of the heaviest eigen- 
value Mh of the Higgs boson mass matrix, and when Ma is larger than Mz, see Sect. 13.3.21 
Indeed, only in this case the self-energies are evaluated at external momenta comparable to 
or larger than the masses circulating in the dominant loops. In addition, if Ma is much larger 

5 Concerning the "full" two- loop effective potential calculation [25], see Sect. 12.3^21 
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than Mz, the relative importance of these corrections decreases, since the tree- level value of 
inn grows with Ma- In fact, the effect of the one-loop momentum-dependent corrections on 
Mh amounts generally to less than 2 GeV. 

Assuming that the relative size of the two-loop contributions follows a pattern similar 
to the one- loop part, we estimate that the two-loop diagrams involving D-term interactions 
should induce a variation in of at most 1-2 GeV, while we expect those with pure gauge 
electroweak interactions to contribute to not very significantly, probably of the order of 
1 GeV or less. Given the smallness of the one-loop contribution it seems quite unlikely that 
the effect of the momentum-dependent part of the 0{a t a s ) corrections to M^, which should 
be the largest among this type of two-loop contributions, could be larger than 1 GeV. As 
already said, the situation can, in principle, be different for the heavier Higgs boson mass. 
The momentum-dependent corrections turn out to be more relevant in processes where the 
H boson appears as an external particle, see Refs. [47,67] and Sect. I3.3T21 



2.5.3 Effects of the variation of the renormalization scale 

Another way of estimating the uncertainties of the kind discussed above is to investigate 
the renormalization scale dependence introduced via the DR definition of tan /3, /i, and 
the Higgs field-renormalization constants [51], see Sect. 12.3141 The variation of the scale 
parameter between 0.5m t and 2rrit is shown in Fig. 12.91 It gives rise to a shift in Mh of about 
±1.5 GeV. This intrinsic error is in accordance with the estimates in Sect. 12.5.21 




Figure 2.9: The renormalization scale dependence of introduced via the DR definition 
of tan/3 and the Higgs field renormalization constants is shown as a function of Ma (left 
plot) and tan/3 (right). The lower curves correspond to tan/3 = 2 (left) and Ma = 100 GeV 
(right). For the upper curves we have set tan/3 = 20 (left) and Ma = 500 GeV (right), 
/ijjp; has been varied from writ/ 2 to 2mt- The other parameters are chosen according to the 
m max sceriar j 0; see the Appendix. The dotted line corresponds to the full on-shell scheme, 
see Sect. HI 
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2.5.4 Estimate of the uncertainties from unknown three-loop cor- 
rections 



Even in the case that a complete two-loop computation of the MSSM Higgs masses is 
achieved, non-negligible uncertainties will remain, due to the effect of higher-order correc- 
tions. Although a three-loop computation of the Higgs masses is not available so far, it is 
possible to give at least a rough estimate for the size of these unknown contributions. 

A first estimate can be obtained by varying the renormalization scheme in which the 
parameters entering the two-loop part of the corrections are expressed. In fact, the resulting 
difference in the numerical results amounts formally to a three-loop effect. Since the 0(a t a s + 
al) corrections are particularly sensitive to the value of the top mass, we compare the 
predictions for Mh obtained using in the two-loop corrections either the top pole mass, 
777,P ole = 174.3 GeV, or the SM running top mass m t , expressed in the MS renormalization 
scheme, i.e. 



pole 
4_ 

1 + 4a s (m t )/3iT - a t {m t )/2-7r 



m t ee m t (m t )fi ~ 1 , A . /"« 7-^7^- (2-79) 



Inserting appropriate values for the SM running couplings a s and a t we find m t = 
168.6 GeV. In Fig. I2.1UI we show the effect of changing the renormalization scheme for m t 
in the two-loop part of the corrections. The relevant MSSM parameters are chosen as in 
Fig.ES i-e. tan/? = 3 , M SUS y = M A = fi = 1 TeV and m- g = 800 GeV. The dot-dashed and 
dotted curves show the 0(a t a s ) predictions for M h obtained using m^ ole or fn t , respectively, 
in the two-loop corrections. The solid and dashed curves, instead, show the corresponding 
0(a.ta s + a>t) predictions for Mh- The difference in the two latter curves induced by the shift 
in mt, which should give an indication of the size of the unknown three- loop corrections, is 
of the order of 1-1.5 GeV. However, as can be seen from the figure, the effect of the shift in 
m t partially cancels between the 0(a t a s ) and 0(otf) corrections, and there is no guarantee 
that such a compensating effect will appear again in the three-loop corrections. 

An alternative way of estimating the typical size of the leading three-loop corrections 
makes use of the renormalization group approach. If all the supersymmetric particles (in- 
cluding the CP-odd Higgs boson A) have the same mass Msusy, and (3 = ir/2, the effective 
theory at scales below Msusy is just the SM, with the role of the Higgs doublet played by 
the doublet that gives mass to the up-type quarks. In this simplified case, it is easy to apply 
the techniques of Refs. [14, 15] in order to obtain the leading logarithmic corrections to Mh 
up to three loops (see also Ref. [68]). Considering, for further simplification, the case of zero 
stop mixing, we find 



3 „ \ t 1 23 2 5 33 2 1 ^ 



(2.80) 

where t = log (M| USY /mj), m t is defined in eq. f|2.79|) . a t and a s have to be interpreted 
as SM running quantities computed at the scale Q = fn t , and the ellipses stand for higher 
loop contributions. It can be checked that, for Msusy = 1 TeV, the effect of the three-loop 
leading logarithmic terms amounts to an increase in Mh of the order of 1-1.5 GeV. If Msusy 
is pushed to larger values, the relative importance of the higher-order logarithmic corrections 
obviously increases. In that case, it becomes necessary to resum the logarithmic corrections 
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Figure 2.10: M h as a function of X t , using either m]? olc or m t in the two-loop corrections. 
The relevant MSSM parameters are chosen as tan (3 = 3 , M SUSY = M A — /j, — 1 TeV and 
rrig = 800 GeV. For the different lines see text. 

to all orders, by solving the appropriate renormalization group equations numerically. Since 
it is unlikely that a complete three-loop diagrammatic computation of the MSSM Higgs 
boson masses will be available in the near future, it will probably be necessary to combine 
different approaches (e.g. diagrammatic, effective potential, and renormalization group), in 
order to improve the accuracy of the theoretical predictions up to the level required to 
compare with the experimental results expected at the next generation of colliders. 

To summarize this discussion, the uncertainty in the prediction for the lightest CP-even 
Higgs boson arising from not yet calculated three-loop and even higher-order corrections can 
conservatively be estimated to be 1-2 GeV. From the various missing two-loop corrections 
an uncertainty of less than 1-2 GeV is expected. However, it is extremely unlikely that all 
these effects would coherently sum up, with no partial compensation among them. Therefore 
we believe that a realistic estimate of the uncertainty from unknown higher-order corrections 
in the theoretical prediction for the lightest Higgs boson mass should not exceed 3 GeV. 
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2.6 Outlook 



The current precision in the theory prediction of M^, with an uncertainty of about 3 GeV 
from missing higher-order corrections and at least 4 GeV from parametric uncertainties can 
certainly not match the anticipated experimental precision at the LHC (Am h ks 200 MeV) 
or at the LC (Am h « 50 MeV). 

The intrinsic uncertainties from unknown higher-order corrections could be reduced to the 
level of 0.5 GeV or below, if a full two-loop calculation (including a proper renormalization) 
and leading/subleading three-loop (and possibly leading four-loop) calculations are available. 
On the other hand, the top quark mass will be measured extremely precisely at the LC [69], 
bringing down the parametric error induced by m t to the level of 0.1 GeV. Also the other 
parametric uncertainties from SM parameters are expected to be reduced to this level. 

If both the intrinsic and the parametric error will reduce as anticipated, the mass of the 
lightest Higgs boson M h can serve as an extremely accurate electroweak precision observable. 
It will be possible to use it to constrain unknown SUSY parameters, especially from the scalar 
top sector, see e.g. Refs. [70-72]. 
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Chapter 3 

The Higgs boson sector of the cMSSM 



The radiative corrections to the Higgs boson masses in the CP conserving MSSM (rMSSM) 
are meanwhile quite advanced, as has been described in detail in Sect. 121 Comparisons of 
different methods, see e.g. Sect. 12.41 showing agreement where expected, have lead to deeper 
insight into the radiative corrections in the MSSM Higgs sector and thus to the confidence 
that the higher-order contribution, although being large, are under control. 

In the case of the MSSM with complex parameters (cMSSM) the higher-order corrections 
have yet been restricted, after the first more general investigations [26], to evaluations in 
the EP approach [27] and to the RG improved one-loop EP method [28,29]. These results 
have been restricted to the corrections coming from the (s)fermion sector and some leading 
logarithmic corrections from the gaugino sector. A more complete calculation has been at- 
tempted in Ref. [30]. More recently the leading one-loop corrections have also been evaluated 
in the FD method, using the on-shell renormalization scheme [31]. Some preliminary results 
including the full one- loop evaluation have been presented in Ref. [32]. The full one- loop 
result including a detailed phenomenological analysis can be found in Ref. [33,47]. 

In this chapter we describe the full one-loop calculation in the FD approach using the on- 
shell renormalization scheme (MS renormalization for tan/3 and the field renormalizations). 
Following Ref. [31], we present all analytical details of the calculation. The results are then 
briefly analyzed in view of the corrections coming from the non-(s)fermion sectors and the 
effects of the non-vanishing external momentum that can only be incorporated completely 
in the FD approach. For a more detailed analysis, see Ref. [33]. All results are incorporated 
into the public Fortran code FeynHiggs2.1 [73]. 

3.1 Tree- level relations and on-shell renormalization 
scheme 

In this section we review the tree-level structure of of the MSSM, however with an emphasis 
on the CP-violating parameters and their effects. 
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3.1.1 The scalar quark sector in the cMSSM 

The mass matrix of two squarks of the same flavor, q^ and (Jr, is given by 

M~-( Ml + m * m « X «\ (31) 



m q X q M R + m q 



with 



Ml = M? + Mf cos 2/3 (If-Q 9 4) 

M^ = M|, + Mf cos 2/3 Q 9 4 (3.2) 
X g = A g — //{cot /3 , tan/3}, 

where {cot /3 , tan/3} applies for up- and down-type squarks, respectively. In the Higgs and 
scalar quark sector of the cMSSM N q + 1 phases are present, one for each A q and one for /i, 
i.e. N q + 1 new parameters appear. As an abbreviation it will be used 

(f) q = arg (X q ) . (3.3) 

As an independent parameter one can trade arg(A 9 ) = (p* for (fi q . 
The squark mass eigenstates are obtained by the rotation 

IH(t) <M> 

with 

V q = ( % S } ) , U«*U* = 11 , (3.5) 



~S q C q 



The mass eigenvalues are given by 

1 



2 2 i 

m~, = m„ H — 

91,2 <? 2 



Ml + M 2 R TJ(Ml-M 2 R y + 4ml\X q \* , (3.6) 



and are independent of the phase of X q } The elements of the mixing matrix U can be 
calculated as 



I Mj + ml -ml x* 

C q — - = , 5g — ■ 

\ m 2 - — m 2 - x M 2 + m 2 — m 2 ~ x m 2 ~ — m 2 ~ 

Y 91 <?2 \l L 1 92 y <?1 92 

The parameter = cos 6^ is real, whereas = e -1 * 9 sin ^ can be complex with the phase 

<f> 8 . = -<p q = arg (X* q ) . (3.8) 



In the limit of having all parameters real, this definition differs slightly from the one given in Sect. 
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3.1.2 The chargino/neutralino sector in the cMSSM 

The physical masses of the charginos can be determined from the matrix 



X 



M 2 sin /3M W 

^/2cospM w /I 



(3.9) 



which contains the soft breaking term M 2 and the Higgs mass term //, both of which may take 
on complex values in the cMSSM. The rotation to the physical chargino mass eigenstates 
is done by transforming the original Wino and Higgsino fields with the help of two unitary 
2x2 matrices U and V, 



XI 

xi 



V 



Xi 

X2 



u 



w 



(3.10) 



These definitions lead to the diagonal mass matrix 

= u*xv+ 



m~± 



m^± 

X-2 



(3.11) 



From this relation, it becomes clear that the chargino masses m ~± and m~± can be determined 
as the (real and positive) singular values of X. The singular value decomposition of X also 
yields results for U and V. 

A similar procedure is used for the determination of the neutralino masses and mixing 
matrix, which can both be calculated from the mass matrix 



/ M 1 


—M z s w cos (3 
\ M z s w sin (3 





M 2 
M z c w cos P 
M z c w sin (3 



-M z s w cos (3 
M z c w cos f3 




M z s w sin (3 \ 
M z c w sin (3 

/ 



(3.12) 



This symmetric matrix contains the additional complex soft-breaking parameter Mi. Its di- 
agonalization is achieved by a transformation starting from the original bino/wino/Higgsino 
basis, 



/X°i] 






xl 


= N 


w° 




xl 




HI 


\x°J 







X-2 



( m x\ 




V o o 



\ 



m^J 



^3 



N* YN+. 



(3.13) 



The unitary 4x4 matrix N and the physical neutralino masses again result from a numerical 
singular value decomposition of Y. The symmetry of Y permits the non-trivial condition 
of using only one matrix N for its diagonalization, in contrast to the chargino case shown 
above. 
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3.1.3 The cMSSM Higgs potential 

The Higgs potential Vh contains the real soft breaking terms m\ and mf, the potentially 
complex soft breaking parameter m\ 2 and the U(l) and SU(2) coupling constants g\ and g^. 



V H = miH^Hu + m z 2 H* 2l H 2i + e tJ (m( 2 H u H 2j + m\ 2 H* U H* 2] ) 



+ Ugf + gl)(H* u H u - H* 2t H 2l f + \gl\H* u H, 



2,\ 



(3.14) 



The indices {i,j} = {1,2} refer to the respective Higgs doublet component, and e±2 = +1. 
The Higgs doublets are decomposed in the following way, 



n 2 



H22 



Vl + 



T2 y 



■A 



V2 + 7|(02 + %; 



(3.15) 



As compared to the real case, see eq. ()2.2|) . besides the vacuum expectation values V\ and 
v 2 , eq. ()3.15|) introduces a real, as yet undetermined, possible new phase £ between the two 
Higgs doublets. Using this decomposition, Vh can be rearranged in powers of the fields, 



V, 



H 



Tfafa - T 02 2 - T xi xi - T X2 X2 

A 



+ \ (01,02,Xl,X2) M n 



W2 

Xi 

W 



+ 



+ 



(3.16) 



where the coefficients of the linear terms are the tadpoles and those of the quadratic terms 
are the mass matrices M n and M c . The tadpole coefficients read 



Tfa = -V2{m\vi + cos£'\m 2 l2 \v 2 + \{g\ + gl){v\ 

T <t>2 = -V2(mlv 2 + cos^VLki - \(j9i + 9i)( v l 
T Xl = — v2 sin£'|mj 2 |wi = — T, 



vl)v 2 ) 



X2' 



(3.17a) 
(3.17b) 
(3.17c) 



H ._ 



£ + arg(mj 



With £ .— t, -r cx,l &V il 12 j 

The real, symmetric 4x4- matrix M n and the hermitian 2x2-matrix M c contain the 
following elements, 



M 



<i>x 



M 



<)>x 



(3.18a) 
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M 



™\ + i(9i + #DO? - v i) cos^|m? 2 | - \{g\ + g 2 2 )v x v 2 
v cos£'|m? 2 l - \{9x + 9i)viv 2 m\ + \{g\ + g$)(3v% - v\) 
— sin | | 



4>x 



sin£'|mf 2 | 

'm\ + \{gl + gl){vl-vl) 
cosf'|m? 2 | 



cosf'|m? 2 | 
m 2 2 + Ugj + g 2 2 )(v% 



J 2 




m 



121 



\g\v\v 2 



1.2 - I^I^l^ 



(3.18b) 
(3.18c) 
(3.18d) 
(3.18e) 



The non- vanishing elements of M<^ x lead to CP- violating mixing terms in the Higgs potential 
between the CP-even fields 0i and 02 and the CP-odd fields xi an d Xi if £' 7^ 0- The 
physically relevant mass eigenstates in lowest order follow from a unitary transformation of 
the original fields, 



U n ( ) 

such that the resulting mass matrices 



(h\ 
H 

A 

\GJ 



02 

Xi 

W 



U 



c(0) 



(3.19) 



Mfs = u n(0 )M n U+ 0) and Mf ^ = U c(0) M c U 



c(0) 



(3.20) 



of the new fields will be diagonal. The new fields correspond to the three neutral Higgs 
bosons h, H and A, the charged pair H ± and the Goldstone bosons G and G . 

The lowest-order mixing matrices can be determined from the eigenvectors of M n and M c , 
calculated under the additional condition that the tadpole coefficients (I3.17|) must vanish in 
order that v\ and v 2 are indeed stationary points of the Higgs potential. This automatically 
requires £' = 0, which in turn leads to a vanishing matrix M0 X and a real, symmetric matrix 
M c . Therefore, no CP- violation occurs in the Higgs potential at the lowest order, and the 
corresponding mixing matrices can be parametrized by real mixing angles as 



U n (o) 



/ — sin a cos a \ 

cos a sin a 

— sin (3 n cos j3 n 

V cos/5 n sin/5 n/ / 



u 



c(0) 



• sin (3 C cos (3 C 
cos p c sin j3 c 



(3.21) 



The mixing angles a, (3 n and /3 C can be determined from the requirement that this transfor- 
mation will result in diagonal mass matrices for the physical fields. 



38 



3.1.4 Higgs mass terms and tadpoles 

The terms in Vu that are linear or quadratic in the physical fields are denoted as follows, 
V H = const. -T h -h-T H - H -T A - A-T G -G 



+ l(h,H,A,G) 



( m l 


m hH 


m hA 






(h\ 


™Ih 




m HA 


™hg 




H 


m hA 


m HA 


m\ 


m AG 




A 


\™>hG 


m HG 


m AG 


™>G ) 







+ 



(3.22) 



+ (H~,G-)- 



m 

2 

mi 



in 



H-G^ 



"G-H+ m G± 

In order to perform a simple renormalization procedure, the parameters in V# have to be 
expressed in term of physical parameters. In total, Vh contains eight independent real 
parameters: v±, t>2, gf, g\ y mf, m 2 , \m\ 2 \ and arg(m 2 2 ) (or £'). The coupling constants g\ 
and (?2 can be replaced by the electric coupling constant e and the weak mixing angle # w 
fs w = sin^w, c w 



H + 
G + 



sin# w , c w = cos# w ), 

e = gi = g 2 s w , 
while the Z boson mass Mz and tan/3 substitute for v± and v-i- 

tan P 

The W boson mass is then given by 



V ^Ml = \{gl + gl){vl + vl). 



M w = 7;9l{vl + vl). 



(3.23) 
(3.24) 

(3.25) 



The tadpole coefficients in the physical basis follow from the original ones (j3.17|) by applying 
the transformation (j3.21|) . 

,2„. „„„ „ 2„ 



T, 



ii - V2[— m\v\ cos a — m^sina — cos^'|m^ 2 |(t>i sin a +t> 2 cosa 
- \{g\ + 92)( v i -«I)(«icosa -v 2 sma)], 
\/2[-\-m\v i sin a — m\v 2 cos a — cos£'|m 2 2 |(t>i cos a; — t> 2 sina! 



(3.26a) 
(3.26b) 



Th = V 'Z[+m\v i sin a — m 2 -y 2 cos a 

+ \{g\ + 9l)(vl - vDiyisuia + v 2 cos a. )], 
Ta = Vzsin^'Im^K^i cos/3 n + t> 2 sin/3 n ), (3.26c) 
T G = -tan{/3-p n )T A . (3.26d) 

Due to the linear dependence of Tq on Ta, eq. (|3.26|) provides only three replacements for 
the original parameters. Typically, the remaining parameter is replaced by either the mass 
of the neutral A-boson, m A , or the mass of the charged pair, m H ±. Their expressions in 
terms of the original parameters are given by 



2 2 • 2 a 1 2 2 a ■ no H I 2 1 

m A = m 1 sm p n + m 2 cos p n — sm2p n cos £ |m 12 | 
-cos 2/3 + </|)(u? 



7 2;> 



m^-± = m 2 sin 2 /3 C + ml cos 2 /3 C — sin 2f3 c cos £'|m 2 2 1 

- cos2f3Mgl + 5f 2 )(w 2 - vj) + \g\{v\ cos/3 c + v 2 sin/3 c ) s 



(3.27a) 
(3.27b) 



39 



If (3 n = (3 C (what will be shown below) the relation between m\ and rn? H ± becomes 

m 2 H± = m\ + M 2 v . (3.28) 

Using (|3.26b -c) and ()3.27a|) or ()3.27b|) . all of m 2 , mf, \m\ 2 \ and £' can be substituted by T h , 
Tjj, Ta and [vn A or m^±). In summary, this leads to the following replacements: 



v 2 



y/2 cos (3 s w c w M z 
e 

V2 sin (3 s w c w M z 



g x -> e/c v 
92 — > e/ s v 



+ 
+ 



- 4m^ + 2(2m 2 A + M|) cos 2/3 + M§ cos(4/3 - 2/3 n ) + Mf cos(2/3 n 
eT h cos/? n 



(3.29) 

(3.30) 

(3.31) 
(3.32) 

(3.33) 



2c w s w M z 
eT H cos /3 n 

2c w s w M z 



cos (3 cos /3 n sin a + sin (3 (cos a cos (3 n + 2 sin a sin /3 n ) ) 
(cos (3 n sin a sin /3 — cos a (cos /9 cos (3 n + 2 sin /3 sin /3 n )) / (cos 2 (/3 — /3 n )) 



Am\ + 2(2m 2 A + M 2 Z ) cos 2/3 + M| cos(4/3 - 2/? n ) + Mf cos(2/3 n ) 



(3.34) 



eT fe sin f3 n 
2c w s w M z 
eTn sin f3 n 
4c w s w M z 



K 2 I - + f 2 ) 

1 



cos /3 sin /3 n sin a + cos a (sin (3 sin /3 n + 2 cos (3 cos /3 n )) 

(2 sin a sin/3sin/3 n + cos/5 (3sin(a — /3 n ) + sin(a + /3 n )) / (cos 2 (/3 — /3 n )) 

(3.35) 

cos(/5 + a) + cos(/3 — a) cos(2/? n )) 



- 2 • fl eT h 
-m A sin 2a — — — 

2 ^ 4cwS w M z 



eT, 



4cwS w M z 



sin 



((3 + a) + sin(/3 - a) cos(2/3 n )) / (cos 2 (/3 - &)) 



eT A 



2s w c w M z (cos /3 cos f3 D + sin /? sin (3 D ) 



(3.36) 
(3.37) 



The resulting physical mass terms are given either in terms of m A or m#±, depending on 
which parameter leads to more compact expressions. 
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The charged Higgs sector contains, apart from m 2 H ±, the mass terms 

m 2 H - G + = —m 2 H± tan(/3 — f3 c ) (3.38a) 
T H sin(a - (3 C )/ cos(/5 - (3 C ) 



2M z Sv,c w 
e 



2M zSvJ cJh cos(a - (3 C )/ cos(/3 - (3 C ) 
2tf ; SwCw ffi/cos(/3-A), 



m G-H+ = ( m H-G+)*> (3.38b) 
m 2 G ± = m 2 H± tan 2 (/9 — /3 C ) 

- mf^-Tv cos (« + ^ - W/ cos 2 (/3 - ft) (3.38c) 

+ sin (« + - 2ft)/ cos 2 (/5 - ft), (3.38d) 

where the star denotes a complex conjugation. 

The neutral mass matrix is more easily parametrized by tha, as can be seen from the 
2x2 sub-matrix of the A and G boson: 

m 2 AG = —m 2 A ta,n(P — ft) (3.39a) 

" 2M z e s vC J h cos(a - ft)/ cos(/3 - ft), 

= m \ tan 2 (/3 - ft) (3.39b) 

+ 2uf^- T ^ sin (« + P~ 2/5n)/ cos 2 (/3 - /3 n ). (3.39c) 
The CP- violating mixing terms connecting the h-/H- and the A-/G-sector are 

T A sin(a- ft)/cos(/3- ft), (3.40a) 

T A cos(a;-/3 n )/cos(/3-/3 n ), (3.40b) 

(3.40c) 

T A sin{a- (3 Q )/cos{(3- (3 n ). (3.40d) 





e 


m hA = 


2M z s w c w 




e 


m hG = 


2M z s w c w 


™ha = 


-mla, 


m HG = 


e 

2M z s„c w 
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Finally, the mass terms of the CP-even h and H bosons are: 



rrii 



rn 



hH 



m 



H 



T h \ sin(a - /9 n )(cos(2a - (3 — (3 n ) + 3 cos((3 - AO)/ cos 2 (/3 - f3 n ] 



M 2 z sm 2 {a + f3) 

+ m 2 A cos 2 (a — 13)/ cos 2 (/3 — /3 n ) 

+ 2m^c v Th cos ( a ~ ^) sin2 ( a _ Ai)/ cos 2 (/3 - /? n ) 

2M z s w c„ J " ft 2 

: -M| sin(a + /3) cos(a + (3) 
+ m 2 ^ sin(a — (3) cos(a — (3)/ cos 2 (/5 — /3 n ) 
+ o,,- 6 . ^ Tg sin(a - /3) sin 2 (a - (3 n )/ cos 2 (/3 - /?„) 
T/ l cos(a — /3) cos 2 (a — /3 n )/ cos 2 (/5 — /3 n ), 
M|cos 2 (a + (3) 
+ sin 2 (a — (3) / cos 2 (/? — /3 n ) 
+ 



(3.41a) 



(3.41b) 



e 

2M z s„c w 
r2 „ 2 



(3.41c) 



2M z s w Cw 
e 

2M z s„c 



— Tjy| cos(o; — /3 n )(cos(2a — (3 — (3 n ) — 3cos(/3 — AO)/ cos 2 (/3 
T ft sin(a — (3) cos 2 (a — f3 n )/ cos 2 (/5 — (3 n ). 



-AO 



3.1.5 Masses and mixing angles in lowest order 

The masses and mixing angles in lowest order follow from eqs. ()3.38|) - (j3.41j) and the additional 
requirement that both the tadpole coefficients T^^^a} and all non-diagonal entries of the 
mass matrices must vanish. From ()3.38a|) and ()3.39a|) . it immediately follows that 

Pc = Pn = A (3.42) 

which in turn determines the tree-level value of a = o;(o) (up to a sign) from (I3.41bl) as 
(which is equivalent to eq. ()2.6|) in the rMSSM) 

tan2a (0) = tan 2(3^ -f , -- < a (0) < . (3.43) 

171 Z 

The Higgs masses mh(o) and rriH(o) are the eigenvalues of their 2x2 mass matrix with en- 
tries f)3.41l a set to zero), 

{m 2 m ,m 2 m } = \(rn\ + M 2 z T yj {m\ + M 2 f - Am 2 A M 2 cos 2 2/?) . (3.44a) 

Finally, combining eqs. ()3.27|) and ()3.42|) relates the remaining masses m A and m H ± with 
each other, 

ml ± =m 2 A + c 2 w M 2 z = m\ + M 2 W . (3.45) 

Specifying one Higgs boson mass as an input parameter therefore unambiguously determines 
the other ones. Since the CP-odd A boson will -due to the CP-violating mixing in the neutral 
Higgs sector- no longer be an eigenstate in higher orders, the charged Higgs mass will 
be used as input parameter in the cMSSM. 
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3.1.6 Renormalization of the Higgs potential 



In this and the following sections we will focus one the one-loop corrections to the cMSSM 
Higgs sector. However, two-loop corrections (taken from the rMSSM, see Refs. [5, 7, 16]) will 
be included in our numerical evaluation, see Sect. 13.31 

In order to calculate the first-order corrections to the Higgs boson masses and effective 
mixing angles, the counter terms for the mass and tadpole coefficients in the Higgs potential 
are needed. Apart from these, counter terms are needed for several other parameters which 
appear in the Higgs potential: 



M 2 Z 
M 2 W 
M n 
IVL 



Mf + 8 M 2 Z , 
+ 8M^ } 
M n + 8M Q , 
M r + 8M r , 



T h ^T h + 8T h , 
T H ->T H + 8T H , 
T A —>T A + 8T A , 
tan (3 — > tan /3 (1 + <5tan/3 



(3.46) 



These definitions explain why the expressions (j3.38H3.4T]) for the Higgs masses must differ- 
entiate between the mixing angles f3 n and (3 C (which, like a, are not renormalized) and the 
parameter j3 (which is). This distinction is necessary to arrive at the following expressions 
for the counter terms. 

The field renormalization matrices of both Higgs multiplets can be set up symmetrically, 



fh\ 

H 

A 
\GJ 



+ \8Z h 

\8Z hA 
\ \8Z hG 



\8ZhH 
1 + \8Z H 

\8Z H a 
\8Z HG 



\8Z h A 

\8Zha 

1 + \8Z AA 



\8Z 



AG 



\8ZhG 

\8Z HG 
\8Z AG 
1 + UZggJ 



fh\ 

H 

A 
\GJ 



(3.47a) 



and 















[g-i 





1 + \8Z H+ \8Z H - G+ 

Uz G - H+ l + Uz G+ 



\8Z* H+ 



\8Z H - 



G+ 



1 + \SZ* r . 



G+. 

















KG- 



(3.47b) 
(3.47c) 



For the mass counter term matrices we use the definitions 
/ 8m\ 



8M r 



5m hH 



8m\ A 



6m hH 
5m hA 



\8m 



2 

hG 



8m 2 H 
8m 2 HA 
8m 2 HG 



8m 2 HA 
8m A 
8m\ G 



8m 2 HG 
8m 2 AG 
8m% , 



8m 2 H± 

Sm G-H+ 



8m 2 H . G , 
8m 2 G± 



(3.48) 



The renormalized self-energies, S(p 2 ), can now be expressed through the unrenormalized 
self-energies, S(p 2 ), the field renormalization constants, and the mass counter terms. This 
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reads for the CP-even part, 



T, hh (p ) = T, hh (p ) + 5Z hh (p - m m ) - 5m h , (3.49a) 

£hn(p 2 ) = ^hH(p 2 ) + SZ hH (p 2 - \{m 2 m + m 2 H(0) )) - 5m 2 hH , (3.49b) 

£hh(p 2 ) = ^hh(p 2 ) + 5Z HH (p 2 - m 2 H{Q) ) - 5m 2 H , (3.49c) 
the CP-odd part, 

±aa(p 2 ) = ^aa{p 2 ) + SZ AA (p 2 - m 2 m ) - 6m 2 A , (3.49d) 

£ AG (p 2 ) = £ AG (p 2 ) + 5Z AG (p 2 - |mi (0) ) - 5m 2 AG , (3.49e) 

Sgg(p 2 ) = S GG (p 2 ) + 6Z GG P 2 - 5m 2 G , (3.49f) 
the CP-violating self energies, 

t hA {p 2 ) = ^ hA {p 2 ) + 5Z hA {p 2 - \{m 2 m + m 2 m )) - 5m 2 hA , (3.49g) 

± hG (p 2 ) = £, G (p 2 ) + 8Z hG (p 2 - ±m 2 (o)) - 5m 2 G , (3.49h) 

^ha(p 2 ) = ^ha(p 2 ) + SZ HA (p 2 - ~(m 2 H{0) + m 2 m )) - 5m 2 HA , (3.49i) 

±hg(p 2 ) = Z HG (p 2 ) + 5Z HG (p 2 - \m 2 m ) - 5m 2 HG} (3.49j) 
and finally for the charged sector: 

?, H - H +(p 2 ) = Y, H - H +(p 2 ) + SZ H - H +(p 2 - m 2 H±{0) ) - 5m 2 H ±, (3.49k) 

£ H - G +(P 2 ) = Z H - G+ (P 2 ) + 5Z H - G+ (p 2 - \m 2 H±{Q) ) - 5m 2 H . G+ , (3.491) 

£ g -h+(p 2 ) = £h- g+ (p 2 ), (3.49m) 

Sg-g+ (p 2 ) = X G - G + (p 2 ) + 5Z G - G+ p 2 - 5m 2 G± . (3.49n) 

It follows from the definition of the field renormalization matrices (j3.47b|) and (j3.47cj) that 
SZ H - H + = 2Re8Z H + and 5Z G - G + = 2 Re SZ G + are real quantities. 

Inserting the renormalization transformation into the Higgs mass terms leads to ex- 
pressions for their counter terms which consequently depend on the other counter terms 
introduced in ()3.46|) . Since the counter terms themselves are of first order, the zero-order 
equalities T^ H A ^ = and (3 n = (3 C = (3 can afterwards be used to simplify these expressions. 

For the CP-even part of the Higgs sectors, these counter terms are: 

5m 2 h = 5m 2 A cos 2 (a - (3) + 8M 2 Z sin 2 (a + (3) (3.50a) 
+ 2A/ z s wCw {ST H cos(a - (3) sin 2 (a - (3) + 5T h sin(a - f3)(l + cos 2 (a - (3))) 
+ <5tan/3 sin/5 cos (3 (m A sin 2 (a - (3) + M§ sin 2 (a + (3)), 
Sm 2 hH = l(5m 2 A sm2{a - /3) - 5M 2 z sm2{a + (3)) (3.50b) 
+ 2M^( ST h sin 3 (" -13)- 5T h cos 3 (a - (3)) 

- 5tan.p sin (3 cos (3 (m 2 A cos 2(a — (3) + M| cos 2(a + /3)), 

5m 2 H = 5m A sin 2 (a - (3) + 5M 2 Z cos 2 (a + (3) (3.50c) 

- 2M z s w c w (STh cos(a - (3){1 + sin 2 (a - (3)) + 8T h sm(a - (3) cos 2 (a - (3)) 

- 5tan/3 sin (3 cos f3 (m A sin 2(a - f3) + M§ sin 2(a + (3)), 
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while for the CP-odd part they follow as 

e 

2Mzs v c 



5m 2 AG = tjtf- 5 {—5Th sin(a — (3) — 5Th cos(a — (3)) — <5tan/3 m 2 A sin (3 cos (3 , (3.50d) 



Sm G = 2Mzs v cS ~ 5Th cos ( a ~ ^) + ^ sin ( a _ ^))> (3.50e) 
for the CP-violating mixing as 

$ m fiA = +2Mi^ 5T Asin(a -/?), (3.50f) 

= + 2M z C SwCw ^ cos ( a ~ (3.50g) 

^ m< HA = - Sm hG> (3.50h) 

= + 2m/. wCw ^ sin ( a - z 3 )' ( 3 - 50i ) 

and for the charged Higgs bosons as 

Sm 2 H - G+ = 2Mzs v cS ~^ Th sin ( a ~ @)~ ST h cos ( a - P) - i $ t a), (3.50j) 

— Stan/3 sin /3 cos (3 , 
5m 2 G - H+ = (5m 2 H - G +)*, (3.50k) 

^ m c± = 2M z l wCw (~^ Tg cos ( a ~ £0 + ^ sin ( a ~ Z 5 ))- (3.501) 

Note that neither <5m^ nor 5m 2 H± are listed here, since one of these masses can be a free 
input parameter whose definition depends on the renormalization. However, from 1)3.45)1 the 
relation 

5m 2 H± = bm\ + SM^ (3.51) 

can be derived between them, which, being generally valid, is used to replace 5m 2 A in other 
expressions. 

For the field renormalization we chose to give each Higgs doublet one renormalization 
constant, 

Hx -> (1 + l6Z Hl )H u H 2 ^(l + \5Z H2 )n 2 . (3.52) 
This leads to the following expressions for the various field renormalization constants: 



5Z 



H~G^ 



8Z hh 


= sm 2 a5Z Hl + cos 2 a 5Z n2 , 


(3.53a) 


SZ A A 


= sin 2 (3 bZ^ + cos 2 /? $Z?{ 2 , 


(3.53b) 


sz hH 


= sin a cos a; {SZ n2 — SZ Hl ), 


(3.53c) 


5Z AG 


= sin (3 cos j3 (SZ^ 2 — 5Z^J, 


(3.53d) 


8Z HH 


= cos 2 a SZ Hl + sm 2 a SZ n2 , 


(3.53e) 


SZ G g 


= cos 2 /? SZui + sin 2 /3 8Z-n 2 , 


(3.53f) 


5Z h -h+ 


= sin 2 /? SZ^ + cos 2 /? 5Z-n 2 , 


(3.53g) 


$Z G - H + 


= sin (3 cos (3 (8Z n2 — SZ ni ), 


(3.53h) 


5Z G - G + 


= cos 2 (3 SZn-i + sin 2 /? 5Zu 2 ■ 


(3.53i) 
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For the field renormalization constants of the CP-violating self-energies it follows, 

5Z hA = 5Z hG = 5Z HA = 5Z HG = 0, (3.54) 

which corresponds to the fact that the CP-violating self-energies do not possess divergences 
depending on the external momentum. 



3.1.7 Hybrid on-shell/MS renormalization 

Up to now, the counter terms for Mz, Mw, rriH±, tan P and T^h,a} as wen as the field renor- 
malization constants are undetermined. For the mass counter terms, an on-shell definition 
is appropriate, 

<5M 2 =ReE zz (M 2 ), SM^ = Re E WW {M^), Sm 2 H± = Re E H - H+ {m 2 H± ). (3.55) 

Here E denotes the transverse part of the self-energy. Since the tadpole coefficients are 
chosen to vanish in all orders, their counter terms follow from T^h,a}{\) + ST^,h,a} = 0: 

5T h = -T ht ft, 5T H = -T H (i), 8T A = -T A{1 y (3.56) 

Concerning the field and tan/3 renormalization, we adopt the DR scheme (see also 
Ref. [49]), 

5tan(3 = 5tan(3 m = [Re E' hh (m 2 m ) - Re E^(m^ (0) )] div , (3.57a) 

5Z Hl =5Z^ = - [Re E' HH ]a=0 ] div , (3.57b) 

5Z m =5Z™ = - [Re E' hh |a=0 ] div , (3.57c) 

i.e. only the divergent parts of the renormalization constants in eqs. ()3.57)1 are taken into 
account. As renormalization scale we have chosen /x^p; = m t . 

3.2 Higher-order corrections 

3.2.1 Calculation of the renormalized self-energies 

In order to obtain the higher-order corrections in the cMSSM Higgs sector the renormalized 
self-energies eqs. ()3.49a|) - (|3.49n|) have to be evaluated. A renormalized self-energy can be 
decomposed as 

E(p 2 ) = E (1 V) + £ (2 V) + ■ ■ ■ , (3-58) 

where E^ denotes the contribution at the i-loop order. In this section we present in detail the 
full one-loop contribution to S(p 2 ), i.e. E^Qo 2 ) in the cMSSM. However, for the numerical 
evaluation in Sect. 13. 3| also corrections beyond the one-loop level (from the rMSSM) are 
taken into account. 

The generic Feynman diagrams for the one-loop contribution to the Higgs and gauge 
boson self-energies are shown in Figs. 13.11 13.21 The one-loop tadpole diagrams entering via 
the renormalization are generically depicted in Fig. 13.31 
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The diagrams and corresponding amplitudes have been obtained with the program Feyn- 
Arts [74,75] and further evaluated with FormCalc [76]. As a regularization scheme differen- 
tial regularization [77] has been used, which has been shown to be equivalent to dimensional 
reduction [50] at the one-loop level [76]. Thus, the employed regularization preserves SUSY. 



3.2.2 Masses and mixing at higher orders 

The masses of particles in a multiplet are determined by the poles of their propagator matrix. 
In higher orders, self-energy terms appear in this matrix. Its inverse is given in the case of 
the three physical neutral Higgs bosons by 

T hHA (p 2 )=p 2 l-M n (p 2 ), (3.59) 
(m 2 m -t hh (p 2 ) -t hH {p 2 ) -t hA (p 2 ) \ 



M n (p 2 



-j: hH (p 2 ) m 2 H{0) - ^ HH {j> 2 ) -Z HA (p 2 ) 
\ -t hA (p 2 ) -t HA (p 2 ) m 2 m -t AA (p 2 )) 



The mixing terms between the Goldstone boson G and physical Higgs bosons will not be 
considered in the actual calculations of this article, since the size of these mixing self-energies 
is absolutely negligible compared to the size of the self-energies containing only physical Higgs 
bosons. 

The loop corrected pole masses correspond to the roots of det ThHA(p 2 )- A full calculation 
therefore involves solving 

det(p 2 l-M n (p 2 )) = (3.60) 

for all its roots. The CP-even Higgs bosons h and H and the CP-odd boson A mix to form 
new mass eigenstates Hi, H 2 and H 3 with 

n^H-L < m 'H 2 — m H z ■ (3.61) 

A simpler approximation for calculating the Higgs masses consists of setting p 2 = 
in (J3.60j) . This u p 2 = 0-approximation" identifies the masses with the eigenvalues of M n (0) 
instead of the true pole masses and is mainly useful for comparisons with effective-potential 
calculations and the determination of effective mixing in higher orders. 

Another simple approximation consists of choosing the p 2 values as follows, 

£ h h(p 2 ) -> t hh (m 2 m ) 

^hh(p 2 ) -> t hh (m 2 H{0) ) (3.62) 
^huip 2 ) -> t hH ((m 2 m +m 2 H{0) )/2) . 

This "p 2 =on-shell" approximation removes all dependencies from the field renormalization 
constants. It results in Higgs boson masses much closer to the true pole masses than the 
"p 2 = 0" approximation, see Sect. 13.3.21 
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Figure 3.1: Generic Feynman diagrams for h, H, A, G self-energies (/ = {e, fi, r, d, s, b 
c, t} ). Corresponding diagrams for the Z boson self-energy are obtained by replacing 
external Higgs boson by a Z boson. 
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Figure 3.2: Generic Feynman diagrams for H ± , G ± self-energies (I = {e, fi, r}, d = {d, s, 
b}, u = {u, c, t} ). Corresponding diagrams for the W boson self-energy are obtained by 
replacing the external Higgs boson by a W boson. 
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Figure 3.3: Generic Feynman diagrams for h, H, A, tadpoles (/ = {e, //, r, d, s, 6, u, c, £}). 



In each of these approximations, the mass eigenstates Hi,H 2 , H 3 are connected to h, H 
and A by an orthogonal transformation matrix U n which diagonalizes M n : 






U Q - \ H\, U n 



m 2 Hl 

m 2 H2 | = U n M n U+. (3.63) 


The elements of U n are used in the following to quantify the extent of CP-violation. For 
example, uf 3 can be understood as the CP-odd part in Hi, while u[ x + u\ 2 make up the 
CP-even part. The unitarity of U n ensures that both parts add up to I. 



3.2.3 The Higgs boson couplings 

The leading corrections in the neutral MSSM Higgs boson sector are taken into account 
by the Higgs boson self-energies at vanishing external momentum (or in the "p 2 =on-shell" 
approximation). The matrix U n then also provides the leading corrections to the neutral 
Higgs boson couplings to SM gauge bosons and fermions, see e.g. Ref. [28]. 

Taking complex phases into account, all three neutral Higgs bosons contain a CP-even 
part, thus all three Higgs bosons can couple to two gauge bosons, VV = ZZ, W + W~ . The 
couplings normalized to the SM values are given by 

QHiW = un sin(/3 - a) + u i2 cos(/3 - a) . (3.64) 
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The coupling of two Higgs bosons to a Z boson, normalized to the SM value, is given by 

9HiHjZ = Ui3 (uji cos(/3 - a) - u j2 sin(/5 - a) ) 

— uj 3 {un cos(/3 — a) — u i2 sin(/3 — a) ) . (3.65) 

The Bose symmetry that forbids any anti-symmetric derivative coupling of a vector particle 
to two identical real scalar fields is respected, gUiUiV — 0. 

Concerning the decay to light SM fermions, the decay width of the if, can be obtained 
from the SM decay width of the Higgs boson by multiplying it with 

(3.66) 
with 

9%uu = ( u *i cos a + u i2 sin a )/s , g v H . uu = u i3 c p / s p (3.67) 

9% ld d = (~ u ii sin « + Mucosa )/cp, g v Hidd = u i3 s p /c p (3.68) 

for up- and down-type quarks, respectively. For more details, see Ref. [73]. 

3.3 Phenomenological implications 

In this section we briefly describe some of the phenomenological implications of a complete 
one-loop evaluation of the cMSSM Higgs sector. More details can be found in Ref. [33]. 

The higher-order corrected Higgs boson sector has been evaluated with the help of the 
Fortran code FeynHiggs2.1 [48,49,73]. The code includes the full one-loop calculation, see 
Sect. 13.2.11 Furthermore, the two-loop corrections are taken over from the rMSSM [5,7,25] 
and for the b/b sector [43,44] from the cMSSM. The code can be obtained from the FeynHiggs 
home page: www.feynhiggs.de . 

For a more detailed phenomenological analysis, constraints on CP-violating parameters 
from experimental bounds, e.g. on electric dipole moments (EDMs), have to be taken into 
account [64]. However, in our analysis below we only take non-zero phases for A t = 
and M 2 ,M 1; which are not severely restricted from EDM bounds. However, our analysis is 
confined to y?^ = 0, since this is the most restricted phase, see e.g. Ref. [78] and references 
therein. 

The numerical analysis given below has been performed in the the following scenario (if 
not indicated otherwise): 

Msusy = 500 GeV, \A t \ = \A b \ = \A T \ = 1000 GeV, VAb = VAt = 0, 
fj, = 1000 GeV, M 2 = 500 GeV, M l = 250 GeV 

M H ± = 150 GeV, tan /3 = 4, 15, fj, m = m t , (3.69) 

where the parameter under investigation has been varied. 

Larger effects of the CP- violating phases are observed for smaller values of m H ± . Values 
of m H ± and tan/3 as given in eq. ()3.69|) could already be challenged by the Higgs search 
performed at LEP [56,60] (depending on the other parameters). It has been shown, however, 



(&//) +04//) 
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that within the cMSSM the obtained limits on cannot be taken over directly to the 
complex case [79,80]. Therefore, we do not apply the bounds of Refs. [56,60], but one 
should be aware that effects for experimentally not excluded parameters might be slightly 
weaker than shown below. 

3.3.1 Dependence on the gaugino phases 

First we analyze the dependence on the gaugino phases <fM x and <fM 2 - 111 Fig. 13.41 the 
dependence of the lightest cMSSM Higgs boson mass on ipM 2 is shown. Am^ := (all 
sectors) - (/// sector) is evaluated at the pure one-loop level for three different values of 
\M^\, IM2I = 500, 1000,2000 GeV from the upper to the lowest line. The other parameters 
are chosen as in eq. ([3.69)1 . The result including the full momentum dependence is given by 
the solid lines, while the p 2 = approximation is shown as dashed lines. In the left plot we 
have chosen tan/3 = 4, in the right one tan/5 = 15. For the low tan/3 value the effects from 
the gaugino (and Higgs) sector are about 2 GeV if the external momentum is not neglected, 
and about 3 GeV if the external momentum is neglected as e.g. in the effective potential 
approach. The effect coming from varying the gaugino phase (fM 2 itself is of (9(1 GeV). 
Both types of effects become smaller for larger tan (3 values. However, being in the ballpark 
of 1-2 GeV the effects from the gaugino sector are non-negligible and have to be taken into 
account in a precision analysis. 

We now turn to the effects from varying ipMi as shown in Fig. 13.51 The parameters are 
as in Fig. E3 but with M 2 = 500 GeV, and \M X \ = 250, 500, 1000 GeV from the most upper 
to the most lower line. The size of the effects from the gaugino sector is of course the same 
as in Fig. 13.41 However, the dependence on <pM\ is niuch smaller, being at (9(200 MeV). 
Aiming to match the anticipated LC accuracy of Sm ^ = 50 MeV even these relatively 
small corrections have to be taken into account. 

3.3.2 Threshold effects for heavy Higgs bosons 

In this subsection the threshold effects on the masses of the heavy neutral Higgs bosons are 
analyzed. In Fig. 13.61 the mass difference Am^ = mH 3 — l m>H 2 is shown as a function of 
m H ± for tan/3 = 4 (left) and tan/3 = 15 (right) for two different values of (fA t , (pA t = tt/2 
(black) and (p^ t = (gray). In general it can be observed that the mass differences are 
larger in the case of non- vanishing complex phases 2 . The dashed curves are evaluated with 
the p 2 = approximation. The mass difference monotonuously decreases with increasing 
m#±. The full calculation shown in the solid curves, however, exhibits strong threshold 
effects coming from the scalar top quarks in the Higgs boson self-energies (e.g. the second 
and sixth diagram in Fig. 13. lj) . Their effects can be larger than ~ 10 GeV and thus can be 
more relevant than the effect induced by the complex phases for A t and A^. Thus, neglecting 
the external momentum can lead to large uncertainties in the calculation of the heavy Higgs 
boson masses. On the other hand, it turns out that the p 2 = on-shell approximation, see 
eq. (|3.62j) . shown as dotted lines, gives a rather good approximation to the full result. The 
remaining deviations stay below the level of 1 GeV. For LC precisions for the heavier Higgs 

2 However, in Ref. [33] it is shown that all mass differences that appear for complex parameters can also 
be realized (for other parameter combinations) in the rMSSM. 
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Figure 3.4: Am^ := (all sectors) - (/// sector) is shown as a function of ipM 2 f° r P 2 7^ 
(solid) and p 2 = (dashed) and for tan/3 = 4 (left) and tan/3 = 15 (right). \M2\ is chosen 
as 500, 1000, 2000 GeV. Only one-loop corrections are taken into account. 
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Figure 3.5: Am^ := (all sectors) - (/// sector) is shown as a function of ifiMx for p 2 ^ 
(solid) and p 2 = (dashed) and for tan/3 = 4 (left) and tan/3 = 15 (right). |M 2 | is chosen 
as 250, 500, 1000 GeV. Only one-loop corrections are taken into account. 
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boson of (9(1 GeV) these corrections should be taken into account. It should be kept in 
mind that the inclusion of the finite widths of the scalar top quarks will change the results 
close to the threshold peaks. These effects will have to be taken into account in order to 
obtain a reliable result at the thresholds. 




all sectors, all sectors, p 2 = 0, all sectors, p 2 = on-shell approx. 

Figure 3.6: Am 32 = m H , A — m H2 is shown as a function of m H ± for tan/? = 4 (left) and 
tan/3 = 15 (right) for two different values of ifA t , fA t = tt/2 (black) and ip At = (gray). 
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Chapter 4 

Higgs boson production at the LC 



Controlling the Higgs boson production properties at the per-cent level is mandatory to 
perform high-precision measurements at the LC. The Higgs boson production cross sections 
themselves will be measured at 0(1%) [37,38]. In order to be able to exploit these mea- 
surements the experimental accuracy has to be matched with a corresponding theoretical 
precision. At the LC, the possible channels for neutral Higgs-boson production are the 
production via Z-boson exchange (Higgs-strahlung), 

e + e~ -> Z^Z{h,H}, 

e + e~ -> Z-^A{h,H}, (4.1) 

and the jyW-fusion channel, 

e + e~ -> v e W + v e W~ -> v e v e {h, H, A} . (4.2) 
Charged Higgs bosons can be produced in pairs, 

e+e~ -> {7, Z} -> H + H- (4.3) 

or singly, 

e + e~ -> {7, Z} -> , (4.4) 

e + e~ -> z/ e H/ + e"{7,Z} -> v e e~ H + (4.5) 

(and also with the charged conjugated processes). We do not discuss here the possibility 
of Higgs-particle production by bremsstrahlung off heavy quarks (e.g. e + e~ — > bb{h, H, A}, 
which can be significant for large tan/3 [81]). 

The Higgs-strahlung processes, eq. (|4.1j) . are possible at the tree-level. The full MSSM 
one-loop corrections can be found in Refs. [82-84]. The leading two-loop corrections have 
been incorporated in Ref. [85]. While the WW fusion channel, eq. (|4.2jh for the CP-even 
Higgs bosons are possible already at the tree-level [86,87], the CP-odd Higgs boson can 
only be produced at the loop level. Within the SM recently the full one-loop corrections 
to the WW fusion channel have been obtained by two groups [88] (see also Ref. [89] for 
a partial analytical calculation). In the MSSM so far only the corrections from all [67] or 
third generation [90] fermion and sfermion loops are available. The CP-odd channel has been 
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evaluated including 3- and 4-point function one-loop corrections [91]. The pair production of 
the charged Higgs bosons, eq. (|4.3|) . is possible at the tree-level, the full one-loop corrections 
in the MSSM can be found in Ref. [92]. The charged Higgs boson production in association 
with a W boson, eq. ()4.4|) . is mediated via loop diagrams, where the full one-loop corrections 
can be found in Refs. [93,94]. Finally the single Higgs production in the gauge boson fusion 
channel, eq. ()4.5|) . possible only at the loop level, has been evaluated including the full 
one-loop SM fermion and scalar fermion corrections in Ref. [95]. 

In the following subsections we will review the leading two-loop corrections to the Higgs- 
strahlung process, the (s)fermion one-loop corrections to the WW fusion channel and the 
(s)fermion one-loop contribution to the single charged Higgs production in gauge boson 
fusion. The main phenomenological consequences are also discussed. 

4.1 Corrections to the Higgs-strahlung channel 

The most promising channels for the production of the CP-even neutral MSSM Higgs bosons 
in the first phase of a LC are the Higgs-strahlung processes [86], 

e + e-^ZH i} (4.6) 

(Hi 2 = h, H) and the associated production of a scalar and a pseudoscalar Higgs boson, 

e + e-^AHi. (4.7) 

We review the computation of the MSSM predictions for the cross sections of both channels 
in the Feynman-diagrammatic (FD) approach using the on-shell renormalization scheme. 
We take into account the complete set of one-loop contributions, thereby keeping the full de- 
pendence on all kinematical variables. The one-loop contributions consist of the corrections 
to the Higgs- and gauge-boson propagators, where the former contain the dominant elec- 
troweak one-loop corrections of 0(a t ), and of the contributions to the 3-point and 4-point 
vertex functions [82-84] 1 . We combine the complete one-loop result with the dominant two- 
loop QCD corrections of 0(a t a s ) [5, 16] and further sub-dominant corrections. In this way 
the currently most accurate results for the cross sections are obtained. 

Furthermore we show analytically that the Higgs-boson propagator corrections with ne- 
glected momentum dependence can be absorbed into the tree-level coupling using the effec- 
tive mixing angle from the neutral CP-even Higgs-boson sector. We compare our results for 
the cross sections with the approximation in which only the corrections to the effective mix- 
ing angle, evaluated within the renormalization-group-improved one-loop effective potential 
approach, are taken into account. For most parts of the MSSM parameter space we find 
agreement of the two approaches of better than 10% for the highest LEP energies, while for 
\/s = 500 GeV the difference can reach 25%. 

1 Only photonic corrections to the Ze + e~ vertex are omitted. These virtual IR-divergent photonic cor- 
rections constitute, together with real-photon bremsstrahlung, the initial-state QED corrections, which are 
conventionally treated separately and are the same as for the SM Higgs-boson production. 
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4.1.1 Cross sections for Higgs production in e + e _ collisions 

Classification of radiative corrections 

The set of diagrams taken into account for Higgs-strahlung e + e~ — ► hZ is schematically 
shown in Fig. 14.11 where a) is the tree-level diagram. The shaded blobs summarize the 
loops with all possible virtual particles, except photons in the Ze + e~ vertex corrections, see 
above. More details can be found in Refs. [82-84]. An analogous set has been evaluated for 
the second process e + e~ — > hA. 

For completeness, in Fig. 14.11 also contributions are shown that are proportional to the 
electron mass or vanish completely after contraction with the polarization vector of the 
Z boson (e.g. the A, G-Z mixing contributions and the longitudinal parts of the Z and j—Z 
self-energies). The different types of corrections can be summarized as follows: 

(i) Corrections to the e, Z, 7 and 7— Z self-energies on the internal and external lines and 
to the (initial state) Ze + e~ and r ye + e~ vertices, b) - g). 

(ii) Corrections to the scalar and pseudoscalar propagators, h). 

(iii) Corrections to the ZZHi (ZAHj) vertex, i). 

(iv) Box-diagram contributions and t-channel-exchange diagrams, j) - 1). 
The corrections (i)-(iv) have a different relative impact: 

- Electroweak corrections of type (i) are typically of the order of a few percent (like in the 
Standard Model) and do not exhibit a strong dependence on any SUSY parameters. 

- The main source of differences between the tree-level and higher-order results are the 
corrections to the Higgs-boson self-energies (ii). They are responsible for changes in 
the physical masses M/, and M# and the effective mixing angle a e ff (via contributions 
to the renormalization constants, Z ext , for the external Higgs particles in the 5*-matrix 
elements, see Ref. [85]) predicted for given values of tan/5 and Ma- At the one- loop 
level these propagator corrections constitute the only source for the large correction of 
0(a t ). At the two-loop level they exclusively give rise to contributions of 0(a t a s ) and 
of 0{af). In this sense the propagator corrections define a closed subset of diagrams, 
being responsible for a numerically large contribution. 

- Corrections to the final-state vertices (iii) are typically larger than those of type (i), 
but smaller than the Higgs-boson propagator corrections. At LEP2 energies they can 
reach at most 7-10% [82] for very low or very large values of tan (3, when the Yukawa 
couplings of the top or bottom quarks become strong. 

- Finally, the box-diagram contributions (iv) depend strongly on the center-of-mass en- 
ergy. They are of the order of 2-3% at LEP2 energies and may reach 20% for yfs = 500 
GeV [84]. 
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It should be noted that initial-state QED corrections as well as finite-width effects (allowing 
for off-shell decays of the Higgs and the Z boson) are not included in our calculation. How- 
ever, by incorporating our result into existing codes, e.g. HZHA [96], QED corrections and 
finite-width effects can be taken into account. 

In Ref. [85] it has been shown analytically that the leading Higgs boson propagator 
corrections can be incorporated by changing the tree-level a in the Higgs-gauge boson cou- 
plings to a e ff, see also Sect. 15.1.21 For instance, for Zh production, the Born coupling 
Vzzh ~ sin(a — (3) is changed by the leading propagator corrections to sin(a e fj — (3). Anal- 
ogous results hold for all Higgs vertices, including the AH-i vertices. While the a e s approxi- 
mation, i.e. using an improved Born result for the cross sections where the tree-level angle a 
is replaced by a e fr, incorporates the dominant one- loop and two- loop contributions, it is ob- 
vious from the discussion above that this approximation neglects many effects included in a 
full FD calculation. These are, in particular, the process-specific vertex and box corrections. 



Cross sections 

In this subsection we briefly summarize the analytical formulae for the cross sections for 
the on-shell production of the Higgs bosons e + e~ — > ZHi, e + e~ — ► AHi including the 
corrections (i)-(iii). Box diagrams (iv) give another, more complicated, set of formfactors 
that make the expressions quite lengthy and are, hence, omitted here; more details can 
be found in Ref. [83]. However, we include the box-diagram contributions, as described 
in Ref. [84], in our numerical programs [85] and in the figures shown in this section. 

The presented formalism for cross sections is general enough to accommodate corrections 
of any order to 2- and 3-point vertex functions. Beyond the one-loop level, however, currently 
only two-loop corrections to the scalar propagators have been included [5,16]. Therefore, 
in the cross section calculations we include all possible types of one-loop corrections and 
the available two-loop corrections to scalar self-energies. This is well justified because, as 
discussed above, propagator corrections constitute a closed subset of the leading 0(a t a s ) 
and 0(a 2 ) contributions. Therefore, these two-loop corrections are of particular relevance 
and interest. 

The cross sections (in the center of mass system (CMS)) for both processes ()4.6|) and 
(J4.7|) have the form: 



da z{A)Hi _ x(s,M* z(A) ,Ml 



dtt 64ttV \D z (s 

where A is the standard phase space factor, 



,2 



(Ai + A 2 cos 2 flcMs) , (4.8) 



A (s,m 2 ,m 2 ) = \J s 2 + m\ + m\ — 2sm\ — 2sml — 2m 2 m 2 ,, (4.9) 
and At, A2 are defined by 

Ai + A 2 cos 2 #cms = \J2 ( MM *) > ( 4 - 10 ) 



pol 
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with M.zs an d M-ps as given below. In the following, 6*cms denotes the scattering angle 
#cms =< M e ~ ) Hi) in the CMS. The momenta of the incoming electron and positron are 
denoted as k\ and k 2 , respectively. The momentum of the outgoing h/H is labeled with q, 
whereas the outgoing Z/A momentum is denoted as p, see Fig. 14.1b . The matrix elements 
for the Higgs-strahlung process and the associated Higgs production read 



M 



zs 



ev(k 2 )Y 



VSb (cv - Sa7 5 ) + V^s 



+ box corrections, 



Dz{s) 



zzs 



£> 7 (a) 



(4.11) 



M 



PS 



ev(k 2 )^^ 



Dz(s) 



+ box corrections. 



T 7(o)w S 7 z( g ) 
Vzps D 7 (s) 



u(ki) 



(4.12) 



For the corresponding expressions for the box contributions see Ref. [83]. 

In the above expressions u(ki) and v(k 2 ) are spinors of the incoming electron-positron 
pair, tfxip) is the polarization vector of the outgoing Z. &v,ca are the renormalized vector 
and axial couplings of the Z boson to an electron-positron pair, at the one-loop level ca = 



-l/4s w c w + c^ 1} , Cy = (-1 + 4s^) /As w c w + c v 



(1) 



cL = 1 - si 



/M 



2 



and S^ z (s) denote the renormalized photon and transverse Z boson self-energies, 
and D~(s) are the inverse Z and photon propagators defined as 



E 7 (s) 
Dz{s) 



D Z (s) 



D 7 (s 



Ml + ±l(s) 
) = s + ±Js) 



(4.13) 



Finally, V denotes the effective neutral Higgs-gauge-boson vertices with the one-loop form 
factors. The explicit expression for those vertices and for the matrix elements for Higgs- 
strahlung and associated Higgs production can be found in the Appendix of Ref. [85] and in 
Ref. [82]. 



4.1.2 Numerical results 

In the following we present numerical examples for the dependence of the neutral Higgs- 
boson couplings and cross sections on tan/3, M^, and the mixing in the scalar top sector. 
The results are obtained in the m™ ax and the no- mixing scenario [53], see the Appendix. 
Only for the scalar soft SUSY-breaking parameters we have set the slepton mass parameters 
to Mf = 300 GeV. The squark mass parameter is denoted as Mg. 

Below we will also perform comparisons with results obtained in the framework of the 
RG improved one-loop EPA, where the input parameters are understood as MS quantities. 
To ensure consistency, in the latter case we have transformed the on-shell SUSY input 
parameters into the corresponding MS values as discussed in Ref. [59]. The results shown 
below for the higher-order corrected Higgs-boson masses and the mixing angle within the 
RG improved one-loop EPA have been obtained with the Fortran program subhpole (based 
on Refs. [14,59]). 
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Differences in the Higgs-production cross sections between our FD result (containing the 
complete one-loop result and the dominant two-loop corrections) and the RG a c s approxi- 
mation have a two-fold origin: the different predictions for the values of Mh and a e Q , and the 
additional contributions contained in the FD result (i.e. the one-loop 3- and 4-point vertex 
functions. 




90 100 110 120 130 140 150 90 1 00 1 1 1 20 1 30 1 40 1 50 



M„ (GeV) M h (GeV) 

Figure 4.2: a Zh and a ah as a function of M h at ^/s = 500 GeV, for two values of M4, in 
the m™ ax scenario. The solid (dot-dashed) line represents the two-loop FD result including 
(excluding) box contributions, the dotted line shows the RG a c s approximation and the 
dashed line shows the one- loop FD result. 

In Fig. I4.2I the cross sections for the Higgs-strahlung process and the associated produc- 
tion are shown for a typical Linear Collider energy, yfs = 500 GeV [37], in the m™ ax scenario 
(in the no-mixing scenario similar results have been obtained). We also show the result for 
the two-loop FD calculation where the box contributions have not been included, in order 
to point out their relative importance for high-energy collisions. For -^s = 500 GeV the 
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sqrt(s) (GeV) sqrt(s) (GeV) 



Figure 4.3: ozh and a Ah as a function of ^fs for tan/? = 5, in the no-mixing scenario. 
The solid (dot-dashed) line represents the two-loop FD result including (excluding) box 
contributions, the dotted line shows the RG a c s approximation and the dashed line shows 
the one-loop FD result. 



differences between the FD result and the RG a c g approximation can be large, an effect 
that is more pronounced for the higher value of Ma- For Ma = 200 GeV, typically they 
are of the order of 10-15% for o Zh and even up to 25% for a Ah (for M h > 90 GeV). The 
difference between the two-loop and one-loop FD result can be sizable. The two-loop result 
for azh is in general larger than the one-loop value, again increasing with Ma, where for 
Ma = 200 GeV the difference can amount up to 15%. a Ah, on the other hand, is decreased 
at the two-loop level for Ma = 200 GeV and the difference may be sizable. 

The box contributions become more important for higher \fs and change the total cross 
section by 5-10%. This result remains unchanged even if sleptons are significantly heavier 
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Figure 4.4: azh and a Ah as a function of Ma at yfs = 500 GeV, shown for tan/5 = 5 
and tan/5 = 50, in the no-mixing scenario. The solid (dot-dashed) line represents the two- 
loop FD result including (excluding) box contributions, the dotted line shows the RG a cS 
approximation and the dashed line shows the one-loop FD result. 



than Mi = 300 GeV used in our numerical analysis, as the dominant contributions to box 
diagrams are given by W and Higgs boson exchanges [83,84], which do not depend on M\. 
Also, one should recall that box contributions lead to an angular distribution of the final- 
state particles different from the effective Born approximation and thus give much larger 
corrections to the differential rather than to the total cross section, at least for some range 
of the scattering angle. Therefore, the box diagrams have a significant effect at Linear 
Collider energies and thus have to be included. The same conclusions can be drawn for 
\fs = 500 GeV in the no-mixing scenario, which we do not show here. The differences 
between the FD result and the RG a e s approximation are only slightly smaller than in the 
m™ ax case. 
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In Fig. 14.31 the results for ozh an d oak are shown as a function of y/s in the no-mixing 
scenario for Ma = 100, 200 GeV. Besides the obvious kinematical drop-off of the cross 
sections, one can observe that the relative differences between the FD two-loop result and 
the RG a e s approximation grow with y/s. The differences remain almost constant or even 
increase slowly in absolute terms, whereas the full cross sections decrease, a Ah becomes 
very small for large Ma, as can be seen in more detail in Fig. 14.41 There we show the 
dependence of o Zh and a Ah on Ma in the no-mixing scenario for tan (3 = 5 and tan (3 = 50. 
For (Tzhi the A boson decouples quickly; the dependence on Ma becomes very weak for 
Ma ^ 250 GeV, when azh is already practically constant (compare e.g. Ref. [5]). In the 
same limit, a ah g° es quickly to zero due to suppression of the effective ZhA coupling, which 
is ~ cos(a e ff — (3); also the kinematical suppression plays a role, but this becomes significant 
only for sufficiently large Ma, Ma > 350 GeV. For large tan/? the decoupling of Ma is even 
more rapid. The differences between the FD two-loop result and the RG a e g approximation 
for the Higgs-strahlung cross section tend also to a constant, but they increase with Ma for 
the associated production. The latter can be explained by the growing relative importance of 
3- and 4-point vertex function contributions compared to the strongly suppressed Born-like 
diagrams. As can be seen from Fig. 14.41 for tan (3 = 50 and Ma > 300 GeV the FD two-loop 
result is almost an order of magnitude larger than the result of the RG a c g approximation, 
and starts to saturate. This can be attributed to the fact that the (non-decoupling) vertex 
and box contributions begin to dominate the cross section value. However, such a situation 
occurs only for very small a Ah values, a Ah ~ 10~ 3 fb, below the expected experimental LC 
sensitivities. 



4.2 Corrections to the WW fusion channel 

While the discovery of one light Higgs boson might well be compatible with the predictions 
both of the SM and the MSSM, the discovery of one or more other heavy Higgs bosons would 
be a clear and unambiguous signal for physics beyond the SM. 

In the decoupling limit, i.e. for Ma ^ 200 GeV, the heavy MSSM Higgs bosons are nearly 
degenerate in mass, Ma ~ M H ~ M H ±. The couplings of the neutral Higgs bosons to SM 
gauge bosons are proportional to 

VVh ~ VHA ~ sin(/3 - a) , (4.14) 
VVH ~ VhA ~ cos(/3 - a) , (V = Z, W±) . (4.15) 

In the decoupling limit one finds (3 — a — > 7r/2, i.e. sin(/5 — a) — > 1, cos(/5 — a) — ► 0. 

At the LC, the possible channels for neutral Higgs-boson production are the production 
via Z-boson exchange, eq. (|4.1j) . and the WW fusion channel, eq. (|4.2j) . As a consequence of 
the coupling structure, in the decoupling limit the heavy Higgs boson can only be produced 
in (H,A) pairs. This limits the LC reach to Mh ^ \/s/2. Higher-order corrections to the 
WW — > H channel from loops of fermions and sfermions, however, involve potentially large 
contributions from the top and bottom Yukawa couplings and could thus significantly affect 
the decoupling behavior. 

In this section we review the one-loop corrections of fermions and sfermions to the process 
e + e~ — > vv{h,H}, i.e. to the production of a neutral CP-even Higgs boson in association 
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with a neutrino pair, both via the VTVT-fusion and the Higgs-strahlung mechanism. In the 
latter case the Z boson is connected to a neutrino pair, e + e~ — > Z{h,H} — > i>ivi{h, H}, 
with I = e, //, r (where the latter two neutrinos result in an indistinguishable final state in 
the detector). 

While the well-known universal Higgs-boson propagator corrections turned out not to 
significantly modify the decoupling behavior of the heavy CP-even Higgs boson, an analysis 
of the process-specific contributions to the WWH vertex has been obtained recently. Taking 
into account all loop and counter-term contributions to the process e + e~ — ► vv {h, H} from 
fermions and sfermions and including also the effects of beam polarization in our analysis, we 
review in this section the LC reach for the heavy CP-even Higgs boson. We have obtained 
results for values of the M SSM parameters according to the four benchmark scenarios defined 
in Ref. [53]. While within these benchmark scenarios we find that the loop corrections do not 
significantly enhance the LC reach for heavy CP-even Higgs boson production and in some 
cases even slightly reduce the accessible parameter space, we have also investigated MSSM 
parameter regions where the loop effects do in fact lead to a significant improvement of the 
LC reach. In "favorable" MSSM parameter regions an e + e~ LC running at y/s = 1 TeV can 
be capable of producing a heavy CP-even Higgs boson with a mass up to M H < 700 GeV. 

Concerning the production of the light CP-even Higgs boson, an accurate prediction of 
the production cross section for precision analyses will be necessary. Aiming for analyses at 
the percent level [38] also requires a prediction of the production cross section in this range 
of precision. Besides the already known universal Higgs propagator corrections, in particular 
loops from fermions and sfermions (especially from the third family) are expected to give 
relevant contributions. We analyze our results for the parameters of the four benchmark sce- 
narios defined in Ref. [53] and study the results as a function of different SUSY parameters. 
We discuss the relative importance of the fermion- and the sfermion-loop contributions and 
furthermore evaluate the fermion-loop correction within the SM for comparison purposes 
(the complete 0(a) corrections within the SM can be found in Ref. [88]). 



4.2.1 Renormalization in the Higgs sector 

The renormalized Higgs-boson self-energies, see eqs. ()3.49aJ) - (j3.49c|) . are as usual given by 
^hh{q 2 ) = T, HH {q 2 ) + 5Z H {q 2 - m 2 H ) - 5m 2 H , 

£hn(q 2 ) = Z hH (q 2 ) + -5Z Hh (q 2 -m 2 H ) + -5Z hH (q 2 -m 2 h ) - 5m 2 hH , 

± hh (q 2 ) = Z hh (q 2 )+6Z h (q 2 -m 2 h )-5m 2 h . (4.16) 

The mass counter terms arise from the renormalization of the Higgs potential, see Ref. [51]. 
They are evaluated in the on-shell renormalization scheme. The field-renormalization con- 
stants, see also eq. ()3.57|) can be obtained in the DR scheme, leading to 

1 div 



5Z H = - [ReT,' HH (m 2 H )] 
5Z h = - [Re E; 



div 



sin a cos a 

oZ hH = {dZ h -dZ H ), 

cos 2a 

8Z Hh = 5Z hH , (4.17) 
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i.e. only the divergent parts of the renormalization constants in eqs. (j4.17j) are taken into 
account. As renormalization scale we have chosen /zj^ = m t . The finite LSZ factors required 
in this renormalization scheme are described in the next subsection. 

4.2.2 The process e + e~ — > v e v e {h, H} 
The tree-level process 

The tree-level process [86,87] consists of the two diagrams shown in Fig. 14.51 Besides the 
W-W-fusion contribution (left diagram), we also take into account the Higgs-strahlung con- 
tribution (right diagram), where a virtual Z boson is connected to two electron neutrinos. 




Figure 4.5: The tree-level diagrams for the process e + e — ► uu{h,H}, consisting of the 
lUIU-fusion contribution (left) and the Higgs-strahlung contribution (right). 



An analytical expression for the tree-level cross section for a SM Higgs boson can be 
found e.g. in Ref. [97]. For relatively low energies and moderate values of the SM Higgs- 
boson mass (y/s < 400 GeV, M HsM < 200 GeV) the resonant production via the Higgs- 
strahlung contribution dominates over the lypU-fusion contribution. At higher energies, 
however, the VFVy-fusion contribution becomes dominant. The cross section, containing 
both contributions, in the high-energy limit takes the simple form [86] 



a e T e 



7T° 



1 + — ^ ] log 



Ml 



2 1 



Ml 

-nSM 



, (4.18) 



where the t-channel contribution from the PUiy-fusion diagram gives rise to the logarithmic 
increase. 

The coefficients for the couplings WWh and WWH are denoted by and Ijy at the 
tree level, respectively (and analogously for the ZZh and ZZH couplings): 



,(o) 

H 



\eM w 

■ sm{p — a) 



ieM w 



cos(/3 — a) . 



(4.19) 
(4.20) 



The SM coupling is obtained by dropping the SUSY factors sin(/3 — a) or cos((3 — a) . 
In the decoupling limit, Ma ^ 200 GeV, (3 — a — >• ir/2, so that sin(/3 — a) — > 1 and 
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cos(/3 — a) — > 0, i.e. the heavy neutral CP-even Higgs boson decouples from the W and 
Z bosons. 

We parametrize the Born matrix element by the Fermi constant, Gp, i.e. we use the 
relation 



2 s w M w 



1/2 



1 + Ar 

where Ar incorporates higher-order corrections, see the next subsection. 



(4.21) 



Higher-order corrections 

In the description of our calculation below we will mainly concentrate on the VTVT-fusion 
contribution. The Higgs-strahlung contribution, which we describe in less detail, is taken 
into account in exactly the same way (with the only exception that a finite Z width has to 
be taken into account). 

We evaluate the one-loop 0(oi) contributions from loops involving all fermions and 
sfermions. Especially the corrections involving third-generation fermions and sfermions, 
i.e. t,b,T,v T , and their corresponding superpartners, ti,t 2 , bi,b 2 , Ti,T2, v t , are expected 
to be sizable, since they contain potentially large Yukawa couplings, y t , y^ y T , where the 
down-type couplings can be enhanced in the MSSM for large values of tan f3. This class of 
diagrams in particular contains contributions enhanced by m1/M w . 

The contributions involve corrections to the WW{h, H} vertex and the corresponding 
counter-term diagram, shown in Fig. I4.6[ corrections to the H^-boson propagators and the 
corresponding counter terms, shown in Fig. 14. 7\ and the counter-term contributions to the 
ev e W vertex as shown in Fig. 14.81 Furthermore, Higgs propagator corrections enter via 
the wave-function normalization of the external Higgs boson, see below. There are also 
H^-boson propagator corrections inducing a transition from the W ± to either G ± or H ± . 
These corrections affect only the longitudinal part of the W boson, however, and are thus 
oc m e /M w and have been neglected. 

While the renormalization in the counter terms depicted in Figs. 14.71 and 14.81 is as in the 
SM (see e.g. Ref. [98]), the WW{h, H} vertices are renormalized as follows, 



WWh : r 



(o),ct 

WWh 



VV VV ll . WWH — 1 H 



(0) 



2 M w 



(4.22) 



. n n cos(/3 - a) 
+ sin p cos p — 



sin(/3 — a) 



i r (0) 

5 tan/? + -5Z h + -^_5Z Hh 



1 

2' 



l + 5Z e 



1SM W 5s v 

2l^ + SZw + T w 



2 r (o) 

1 h 



<0) 



(4.23) 



- sm/3 cos (3 — f- L 5tan/3 + -5Z H + --^5Z 

cos(p — a) 2 2 YvJ 



hH 



Analogous expressions are obtained for T z ^ ($ = h,H). In the above expressions 5Z e 
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Figure 4.6: Corrections to the WW{h, H} vertex and the corresponding counter-term 

diagram. The label / denotes all (s)fermions, except in the presence of a /', in which 
case the former denotes only the isospin-up and the latter the isospin-down members of the 
(s)fermion doublets. 



incorporates the charge renormalization and the Ar contribution arising from eq. (j4.21|) . 

5Z e = 5Z e - ^Ar, 5Z e = JlP(O) - — » ( 4 - 24 ) 



2 ' 2 cv M 



where E T denotes the transverse part of a self-energy. <5M^ is the PF-mass counter term, 5Zw 
is the corresponding field-renormalization constant, and Ss w denotes the renormalization 
constant for the weak mixing angle. The field-renormalization constants, 5Zh, SZ^, and 
SZffh — 8Z hH are given in eq. (|4.17j) . The counter term for tan/5 (with tan/? — > tan (3(1 + 
5 tan (3)) is derived in the DR renormalization scheme [51] (using DRED). The parameter 
tan/3 in our result thus corresponds to the DR parameter, taken at the scale /j,^ = m t . We 
list here all contributing counter terms except for the Higgs field renormalization which has 
already been given in eq. (|4.17j) : 

8M 2 W = ReS^(M^), (4.25) 
6 Ml = ReE|(M|), (4.26) 



6s w \cl(bM\ 5M^ 



2 si \ M 2 Z M 2 W 



(4.27) 



aan/3 = 5 tan /3 DR = - ^ ^ — [Re E' hh (m 2 h ) - Re T! HH {m 2 H )} ^ . (4.28) 
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Figure 4.7: Corrections to the VF-boson propagator and the corresponding counter-term 
diagrams. The label / denotes all (s)fermions, except in the presence of a /', in which 
case the former denotes only the isospin-up and the latter the isospin-down members of the 
(s)fermion doublets. 




Figure 4.8: Counter-term contributions entering via the eu e W vertex. 



For 5Z e we find, taking into account only contributions from fermion and sfermion loops, 

^-2rUt"^M as Jr (429) 

The gauge-boson field-renormalization constants, 5Z W , SZ Z , SZ lZ , drop out in the result for 
the complete S- matrix element. 

In order to ensure the correct on-shell properties of the outgoing Higgs boson, which are 
necessary for the correct renormalization of the S- matrix element, furthermore finite wave- 
function normalizations have to be incorporated, see below. For the various checks that have 
been performed to ensure the reliability of our result, see Ref. [67]. 
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The individual contributions from fermions and sfermions constitute two subsets of the 
full result which are individually UV-finite. This is in contrast to the evaluation of renormal- 
ized Higgs boson self-energies (see e.g. Refs. [10,11]), where a UV-fmite result is obtained 
only after adding the fermion- and sfermion-loop contributions. 



The Higgs-boson propagator corrections and the effective Born approximation 

For the correct normalization of the S'-matrix element, finite Higgs-boson propagator cor- 
rections have to be included such that the residues of the outgoing Higgs bosons are set to 
unity and no mixing between h and H occurs on the mass shell of the two particles. The cor- 
rections affecting the Higgs-boson propagators and the Higgs-boson masses are numerically 
very important. Therefore we go beyond the one-loop fermion/sfermion contribution used 
for the evaluation of the genuine one-loop diagrams and include Higgs-boson corrections also 
from other sectors of the model [10,11] as well as the dominant two-loop contributions [5,7] 
as incorporated in the program FeynHiggs [48]. 

For the WW{h, H} vertex, these contributions can be included as follows, yielding the 
correct normalization of the S matrix at one- loop order 2 : 



WWh : 
WWH : 

This gives rise to the following terms: 

wwh ■ r WF 

W Vv II . 1 W Wh 



r (0) , 1 7 r (0) 

L h + 2 Hh h 



(0) , l y p(0) 

2 ZhH 1 h 



z H [ r% + ^z 



(4.30) 
(4.31) 



WWH : 



1 WWH 



r (o) 

1 h 



H 



IT 



2 r (o) 
1 h 



Z h Z 



HI, 



IT 



(0) 



2 r (o) 

1 H 



z„z 



H ^hH 



(4.32) 
(4.33) 



Analogous expressions are obtained for r^J $ ($ = h, H). In the above expressions, the finite 
Higgs-mixing contributions enter, 



ZhH 



Re± hH (M 2 h 



MZ-m 2 H + ReZ HH (M 2 h 



Ret hH {M%) 



M 2 H -m 2 h + ReZ hh (M 2 H ) 



(4.34) 
(4.35) 



involving the renormalized self-energies S(g 2 ), see eq. f!4.16|) . which contain corrections up to 
the two-loop level. The wave-function normalization factors Zh, Zh are related to the finite 



2 Note that our notation is slightly different from Refs. [85,99] 
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residue of the Higgs-boson propagators: 



>i2 



1 

1 + ReS' (a 2 ) ( ( Ref: /f (<?2) ) 2 



1 

1 + ReS' (a?) ( ( Rea ^ fa2) ) 2 
i + KeL ™ i9j U-^+ReW) 



? 2 =M| 



(4.36) 



(4.37) 



' denotes the derivative with respect to the momentum squared. 

If in eqs. ()4.34j) - ()4.37|) the renormalized self-energies were evaluated at q 2 = 0, the above 
wave-function correction would reduce to the a e fi approximation [85,99]. In this approxima- 
tion, however, the outgoing Higgs boson does not have the correct on-shell properties. 

In order to analyze the effect of those corrections that go beyond the universal Higgs 
propagator corrections, we include the Higgs propagator corrections according to eqs. (|4.32J1 - 
()4.37|) into our Born matrix element, see the next subsection. Concerning our numerical 
analysis, see Sect. 14.272*1 we either use this Born cross section (thus the difference between 
our tree-level and the one-loop cross sections indicates the effect of the new genuine loop 
corrections), or we use the a e g approximation (so that the difference between the tree- level 
and the one-loop cross section directly shows the effect of our new calculation compared to 
the previously used results). 



The higher-order production cross section 

The amplitude for the process e + e~ — > v e v e {h, H} is denoted as 

M<£] e , ($ = h,H; * = 0,1), (4.38) 

where i = denotes the lowest-order contribution and i — 1 the one-loop correction. 

The tree-level amplitude involves the VTiy-fusion channel (left diagram of Fig. 14.5") and 
the Higgs-strahlung process (right diagram of Fig. 14.5") where the virtual Z boson is con- 
nected to two electron neutrinos. As explained above, we include the Higgs propagator 
corrections into our lowest-order matrix element. We use 

Mt\ = M% cc e +M^ F e , (4.39) 

where M.% e ^ is the contribution of the two tree-level diagrams, parametrized with a = 
octree, eq. ()2.6|) . and M.$f^ denotes the wave-function normalization contributions given in 
eqs. ()4.32|) - (|4.37|) (and analogously for the ZZ{h,H} vertices). 

At one-loop order (i = 1), the diagrams shown in Figs. I4.6TH4.8I contribute (and corre- 
sponding diagrams for the Higgs-strahlung process), involving fermion and sfermion loops. 
The counter-term contributions given in eqs. (|4.22|) . (J4.23)) enter via the WW{h, H} vertices 
(and analogously for the ZZ{h,H} vertices), while the other counter-term contributions 
have the same form as in the SM. 



71 



In order to evaluate the cross section that is actually observed in the detector, e + e~ — > 
(h, H + missing energy), we furthermore take into account the amplitude of the Higgs- 
strahlung process where the Z boson is connected to UfUf (/ = //, r), 

M% = M% c ] + M^ f , = h,H; i = 0, 1; / = r) . (4.40) 

Of course there is no interference between the M-^ f for different flavors. 

For all flavors, on the other hand, the Z-boson propagator connected to the two outgoing 
neutrinos can become resonant when integrating over the full phase-space, and therefore 
a width has to be included in that propagator. We have incorporated this by using the 
running width in the Z-boson propagators, T z{s) = (s/M|) T^ p , where T c ^ v = 2.4952 GeV, 
and dropping the imaginary parts of the light-fermion contributions to the Z-boson self- 
energies. 

The cross-section formulas for h production thus become 

< « E K/l 2 ' (4-41) 
« E (M£/I a + 2 Re i( M h!fT M§] ) . (4.42) 

The formulas for if production are analogous, except that we have also included the square 
of the one-loop amplitude. This is because the decoupling behavior of the WWH coupling 
can make the tree-level cross section very small so that the square of the one-loop amplitude 
becomes of comparable size: 

4« E \M { Z\\ (4-43) 

f=e,n,T 

4oc E (\MZ + M% f | 2 ) . (4.44) 

/=e,/u,T 

In this way, at 0(a 2 ) only contributions ~ (M-HfT •^■hj are neglected, which are expected 
to be very small, if M-Hf ls suppressed. 

Parameters for the numerical evaluation 

For the numerical evaluation we followed the procedure outlined in Ref. [100]: The Feyn- 
man diagrams for the contributions mentioned above were generated using the FeynArts [74] 
package. The only necessary addition was the implementation of the counter terms for the 
VV{h,H} (V = W,Z) vertices, eqs. (jO^ - fjOSjl . into the existing MSSM model file [75]. 
The resulting amplitudes were algebraically simplified using FormCalc [76] and then auto- 
matically converted to a Fortran program. The LoopTools package [76, 101] was used to 
evaluate the one-loop scalar and tensor integrals. The numerical results presented in the 
following subsections were obtained with this Fortran program. 3 

3 The code is available at www.hep-processes.de . 



72 



While we have obtained results both for the total cross sections and differential distri- 
butions, in the numerical examples below we will focus on total cross sections only. For the 
various cross checks of our results with (existing) literature, see Ref. [67] and the second and 
third reference in [88]. 

For our numerical evaluation we have chosen for simplicity a common soft SUSY-breaking 
parameter in the diagonal entries of the sfermion mass matrices, Msusy ; and the same 
trilinear couplings for all generations. Our analytical result, however, holds for general 
values of the parameters in the sfermion sector. The SM parameter are chosen as in the 
Appendix. 

The further SUSY parameters entering our result via the Higgs boson propagator cor- 
rections are the SU(2) gaugino mass parameter, M 2 , (the U(l) gaugino mass parameter is 
obtained via the GUT relation, M x = (5/3) (4/4) M ?)i and the g luino 

mass, trig. 

For our numerical analyses we assume all soft SUSY-breaking parameters to be real. Our 
analytical result, however, holds also for complex parameters entering the loop corrections 
to e + e~ — > uu{h,H}. 

4.2.3 SM Higgs-boson production 

For comparison purposes, we start our analysis with the fermion-loop corrections to the 
process e + e~ — > uuHsm in the SM. For the full one-loop calculation see Ref. [88]. 

Fig. 14.91 shows the tree-level and one-loop-corrected production cross section for a SM 
Higgs-boson mass of M HsM = 115 GeV. The absolute values are shown in the upper plot. 
The sharp rise in the cross section for ^fs > 200 GeV is due to the threshold for on-shell 
production of the Z boson in the Higgs-strahlung contribution, see the right diagram of 
Fig. 14.51 Above the threshold the 1/s behavior of the Higgs-strahlung contribution competes 
with the logarithmically rising t-channel contribution from WW fusion. 

The lower plot shows the relative correction coming from all fermions, as well as the 
correction from the third-generation fermions only. The correction from all fermions ranges 
from about +5% at low ^/s to —1.2% at high y^. Restricting to the contribution of third- 
generation fermions only, we obtain corrections in the range from +1.3% to —1.8%. These 
corrections (both from the third family only as well as from the first two generations) are 
at the level of the expected sensitivity for the H^H^-fusion channel at the LC 4 . For a LC 
running in its high-energy mode with ^fs ~ 800 GeV, in particular, a measurement of the 
total cross section with an accuracy of better than 2% seems to be feasible [38] . 

In Fig. 14.101 the SM production cross section is shown as a function of M HsM for s/s = 
800, 1000 GeV. A SM Higgs boson possesses a relatively large production cross section, 
O(10 fb), depending on the available energy, even for M HsM > 500 GeV. Thus it should 
easily be detectable at a high-luminosity LC. 

4.2.4 Light CP-even Higgs-boson production 

Since the mass of the lightest CP-even Higgs boson in the MSSM is bounded from above 
by Mh ^ 140 GeV [5,7], its detection at the LC is guaranteed [102]. In order to exploit the 

4 The full O(a) corrections in the SM can be even larger due to the QED/bremsstrahlung contribu- 
tions [88]. 



73 



350 fb 



300 fb - 



150 fb - 



100 fb - 



50 fb 



fb 




6% 



m - 



0% I 



-m 



all fcl lui atia 
3 rd gen. anlj- 



L --. 



300 



400 



600 



800 



1000 



Figure 4.9: 
and crir 

-HSM -HgM 



The tree-level and the one-loop cross sections for e + e~ — > uuH^m i n the SM, 
are shown as a function of \fs for M HsM = 115 GeV. The upper plot shows 
the absolute values, the lower plot shows the relative corrections for all fermions and for the 
third-generation fermions only. 



precision measurements possible at the LC, a precise prediction at the percent level of its 
production cross section (and its decay rates) is necessary. 

In the following we analyze the h production cross section. To begin with, we focus on 
the four benchmark scenarios given in the Appendix [53] (proposed for MSSM Higgs-boson 
searches at hadron colliders and beyond). Ma and tan/3 are kept as free parameters. 

Figure (|4.11j) shows in the four benchmark scenarios the h production cross section, cr^, 
eq. (|4.42|) . as well as the relative size of the loop corrections, 



Rl 



(4.45) 



(which is nearly always negative), including the contributions from all generations as well as 
from the third generation only. The results are shown in the M^-tan f3 plane for 100 GeV < 
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Ma < 400 GeV and 2 < tan/3 < 50. For larger Ma values the behavior of the lightest 
CP-even MSSM Higgs boson is very SM-like, i.e. the results hardly vary with Ma any more. 

For very low values of Ma, Ma < 150 GeV, the cross section is relatively small. This is 
due to the fact that the WWh coupling at tree level, being ~ sin(/3 — a) , can become very 
small. In this region of parameter space, however, the heavy CP-even Higgs boson is still 
very light and couples to the gauge bosons with approximately SM strength, the tree-level 
coupling of WWH being ~ cos(/3 — a) ~ 1. 

For the interpretation of the two middle and the right column of Fig. 14. Ill it is important 
to keep in mind that we have absorbed the universal Higgs propagator corrections, which are 
numerically very important, into our tree- level cross section. Thus, the relative corrections, 
R\, shown in Fig. 14. Ill display the effects of the other genuine one-loop corrections only. 
We first compare the corrections in the two cases where the Higgs propagator corrections 
are implemented according to eqs. ()4.40j) . ()4.41j) . which ensures the correct on-shell proper- 
ties of the outgoing Higgs boson (second column), and where an a e ff approximation is used 
(third column), which is often done in the literature. In the a c fr approximation, the leading 
contribution of the process-independent corrections entering via the Higgs-boson propaga- 
tors is included by replacing the tree-level coupling of {h, H}WW = {sin,cos}(/5 — a) by 
{sin, cos}(/3 — a e fr). The difference between the full on-shell prescription and the a e s approxi- 
mation turns out to be sizable. It amounts to several percent even for relatively large values 
of Ma- As a consequence, including the Higgs propagator corrections in an a e fj approxima- 
tion will not be sufficient in view of prospective precision measurements of the e + e~ -^vvh 
cross section. 

The results in the second column of Fig. 14.111 show that the size of the corrections from 
fermion and sfermion loops is somewhat different in the four scenarios. While corrections 
of more than 5% only occur for Ma ^ 130 GeV, we obtain corrections of 2-5% in the 
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Figure 4.11: Results for the cross section for e + e~ — > vvh at one-loop order, cr^ (left 
column), the relative corrections including all three fermion and sfermion generations (two 
middle columns), and the relative corrections including only the third generation (right 
column) are shown in the M^-tan/? plane for four benchmark scenarios at y/s = 1 TeV. 
The universal Higgs propagator corrections have been absorbed into the tree-level cross 
section and thus do not appear in the relative corrections. In the second and fourth column 
the Higgs propagator corrections are implemented according to eqs. (|4.4()jl . 1)4.41)1 . while in 
the third column the a e ff approximation is used. In the results for o\ (left column) the 
black region corresponds to <j\ < 50 fb, the dark-shaded region to 50 fb < o\ < 100 fb, the 
light-shaded region to 100 fb < o\ < 150 fb, and the white region to o\ > 150 fb. For the 
relative corrections (second to fourth column) the black region corresponds to R\ > 5%, the 
dark-shaded region to 2% < R\ < 5%, the light-shaded region to 1% < R\ < 2%, and the 
white region corresponds to R\ < 1% (see text). 
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whole parameter space of the no-mixing scenario. Corrections of 1-2% can be found in large 
parts of the parameter space of the m™ ax and the small-a e ff scenario. The situation in the 
four benchmark scenarios, which have been chosen to represent different aspects of MSSM 
phenomenology, shows that the corrections investigated here are typically of the order of 
about 1-5%. A measurement of the e + e~ — * vv h cross section at the percent level will thus 
be sensitive to this kind of corrections. 

In the right column of Fig. 14.111 we show derived including the contributions from 
the third family of fermions and sfermions only. Thus the differences between the second 
and the fourth column reflect the relevance of the loop corrections coming from the first two 
families. While in the m™ ax and the no-mixing scenario differences can mostly be found for 
small Ma, Ma ^ 150 GeV, in the other two scenarios the effect of the first two families can 
be relevant also for larger Ma- Within the gluophobic-Higgs scenario, the first two families 
play a role for small tan/3 and large Ma- In the small-a e fj scenario differences can be found 
for larger tan f3 over the whole Ma range. In the latter two scenarios, the corrections coming 
from the first and second family lead to a partial compensation of the corrections from the 
third family. The light fermion generations can give rise to a contribution of ~ 1%, which is 
non-negligible for cross section measurements at the percent level. 

The Higgs propagator corrections, which we have absorbed into our tree-level cross sec- 
tion, mainly affect the numerical value of Mh, which enters the final-state kinematics, while 
the numerical effect of the corrections to the WWh coupling is less important. The compar- 
ison between the prediction for e + e~ — > vv h in the MSSM and the corresponding process in 
the SM for the same value of the Higgs boson mass (which is not shown here) yields devia- 
tions of more than 5% for Ma ^ 200 GeV, which to a large extent are due to the suppressed 
WWh coupling in the MSSM case. Deviations of more than 1% are found in all scenarios 
up to rather large values of Ma- 

Fig. 14 .121 shows our results for e + e~ — > vv h in the four benchmark scenarios as a function 
of y/s for Ma = 500 GeV and tan f3 = 3, 40. Note that the difference in the cross sections for 
the four benchmark scenarios for given Ma and tan {3 is entirely due to SUSY loop corrections 
(which, as explained above, affect in particular the value of Mh)- 

The numerically important effects of the Higgs propagator corrections become apparent 
in particular from Fig. 14.131 where the tree-level and the one-loop cross sections are shown 
as a function of X t , i.e. the mixing in the scalar top sector. The plots are given for the four 
combinations of Ma = 150,500 GeV and tan/3 = 3,40, and the other parameters (besides 
X t ) are chosen as in the m™ ax scenario (see the Appendix). The variation of the tree- level 
cross sections indicates the effect of the Higgs propagator corrections affecting both the 
value of Mh and the Higgs coupling to gauge bosons. These corrections can change the cross 
section by up to ~ 25%, while the other loop corrections typically stay below 2.5%. 

Finally, in Fig. 14.141 we analyze the relative importance of the purely sfermionic loop 
corrections (corresponding to the Feynman diagrams with sfermion loops in Figs. !4.fiH4.8| 
as before, the Higgs propagator corrections absorbed into the tree-level result contain both 
fermion- and sfermion- loop contributions). These corrections constitute, as explained earlier, 
a UV-finite and gauge-invariant subset of the loop contributions. The relative size of the 
sfermion corrections as compared to the purely fermionic one-loop corrections is shown in 
Fig. 14.141 The upper row shows the relative size as a function of X t in the m™ ax scenario for 
all combinations of M A = 150, 500 GeV and tan/? = 3,40. While for X t rs 0, i.e. for small 
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Figure 4.12: The tree-level and the one-loop cross sections for e + e — > vv h, a° and a\, 
are shown as a function of y/s in the four benchmark scenarios for Ma = 500 GeV and 
tan/3 = 3,40. 



splitting in the scalar top sector, the sfermionic corrections are small, their contribution 
becomes more important for increasing \X t \. They have the opposite sign of the purely 
fermionic corrections and thus partially compensate their effects. For |X 4 /Msusy| ~ 2 the 
sfermionic corrections are about half as large as the fermion corrections. For very large \X t \ 
(which also lowers Mh substantially) they can become even bigger than the fermionic ones. 

In the lower part of Fig. 14.141 (middle and lower row) we analyze the relative size of the 
purely sfermionic corrections in the m™ ax (middle) and the no-mixing scenario (lower row) as 
a function of Msusy- I n the no-mixing scenario, for increasing MgusY the relative size of the 
sfermion corrections becomes smaller, as can be expected in the decoupling limit [103,104]. 
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Figure 4.13: The tree-level and the one-loop cross sections for e + e~ — > uuh, and a\, 
are shown as a function of X t . In the upper (lower) plot Ma has been set to Ma = 150 
(500) GeV, tan/? is fixed at tan/5 = 3,40. The other parameters are as given in eq. (J2J). 



In the m™ 5 ™ scenario, however, the situation is different. Here X t pa A t is fixed to A t pa X t = 
2Msusy- In the htit2 vertex, being ~ A t cosa + /isma , the coupling is proportional to 
the SUSY mass scale. This results in a term ~ A t /MsvsY in the one-loop corrected WWh 
vertex, which for large Msusy goes to a constant and can be of the order of the purely 
fermionic correction. The cross section then behaves as ~ X 4 2 /M| USY as can be seen in the 
upper row of Fig. 14.141 

4.2.5 Heavy CP-even Higgs-boson production 

We now investigate the effects of loop corrections on the cross section for heavy CP-even 
Higgs-boson production in the MSSM. As explained above, these corrections are of particular 
interest in the decoupling region, i.e. for large values of Ma- If Ma ^ the heavy Higgs 

bosons can be pair-produced at the LC via e + e~ — > Z* — > HA. Beyond this kinematical 
limit, H production is in principle possible via the lyiy-fusion and the Higgs-strahlung 
channels. This production mechanism is heavily suppressed at tree-level, however, owing to 
the decoupling property of the H coupling to gauge bosons. If loop-induced contributions 
turn out to be sizable in the mass range Mjj > y/s/2, an enhanced reach of the LC for 
H production could result. 
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Figure 4.14: The relative size of the purely sfermionic loop corrections to o\ as compared 
to the purely fermionic one-loop corrections is shown as a function of X t (upper) and MgusY 
(middle and lower row) for all combinations of Ma = 150, 500 GeV and tan (3 = 3, 40. The 
other parameters are chosen as in the m™ ax (upper and middle) or as in the no-mixing 
scenario (lower row). 
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Figure 4.15: The tree-level cross section for e + e~ — ► vv H , o"^(a e fr), using an a cff ap- 
proximation (left), and the one-loop cross sections, cr^, without (middle), and with beam 
polarization (right column), are shown in the M^-tan/3 plane for four benchmark scenarios 
for ^fs = 1 TeV. The different shadings correspond to: white: a < 0.01 fb, light shaded: 
0.01 fb < a < 0.02 fb, dark shaded: 0.02 fb < a < 0.05 fb, black: a > 0.05 fb. 
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In Fig. 14.151 we first compare the tree-level cross section evaluated in the a c s approxi- 
mation, a H (a e g) (left column), and the one- loop cross section according to eq. (|4.44|) . a x H 
(middle column), in the four benchmark scenarios. For phenomenological analyses of MSSM 
Higgs-boson production in this channel the cross section has so far mostly been evaluated 
using an a e g Born approximation (see also Sect. 14.2.41) . 

We concentrate on the case of \fs = 1 TeV. Since we are interested in Mh ~ Ma > 
•y/s/2 = 500 GeV, we focus on the region 400 GeV < Ma < 900 GeV, and scan over 
the whole tan/3 region. For a LC like TESLA, the anticipated integrated luminosity is of 
(9(2ab _1 ). For this luminosity, a production cross section of about an = 0.01 fb constitutes 
a lower limit for the observation of the heavy CP-even Higgs boson. (For the scenarios 
discussed below, the dominant decay channel of H is the decay into tt or bb, depending on 
the value of tan (3, and also sizable branching ratios into SUSY particles are possible in some 
regions of parameter space; a more detailed simulation of this process should of course take 
into account the impact of the decay characteristics on the lower limit of observability, while 
in this work we use the approximation of a universal limit.) The area with an < 0.01 fb 
is shown in white in Fig. 14.151 In the four benchmark scenarios shown in Fig. 14.151 the 
inclusion of the loop corrections that go beyond the a e g Born approximation turns out to 
have only a moderate effect on the area in the M^— tan (3 plane in which H production could 
be observable 5 . While for the m™ ax and the no-mixing scenario the area is slightly decreased 
to smaller tan (3 values and somewhat enlarged to higher Ma values, the area is slightly 
decreased in tan (3 in the gluophobic-Higgs scenario (and stays approximately the same in 
Ma), while the area is slightly enlarged both in Ma and tan/5 in the small-a; e ff scenario. 
For the four benchmark scenarios, an observation with Mh > 500 GeV is only possible for 
low tan/5, tan/5 < 5, where the LC reach in Mh can be extended by up to 100 GeV. It 
should be noted at this point that while the area of observability is modified only slightly 
in the plots as a consequence of including the loop corrections, the relative changes between 
the tree-level and the one-loop values of the cross sections can be very large, owing to the 
suppressed WWH coupling in the tree-level cross section. 

The prospects for observing a heavy Higgs boson beyond the kinematical limit of the 
HA pair production channel become more favorable, however, if polarized beams are used. 
The cross section becomes enhanced for left-handedly polarized electrons and right-handedly 
polarized positrons. While a 100% polarization results in a cross section that is enhanced 
roughly by a factor of 4, more realistic values of 80% polarization for electrons and 60% 
polarization for positrons [105] would yield roughly an enhancement by a factor of 3. The 
right column of Fig. 14.151 shows the four benchmark scenarios with 100% polarization of 
both beams. The area in the M^-tan/3 plane in which observation of the H boson might 
become possible is strongly increased in this case. In the m™ ax and the no-mixing scenario, 
H observation could be possible for small tan/3 up to Ma ^5 700 GeV. In the gluophobic- 
Higgs and the small-a e ff scenario this effect is somewhat smaller. In the latter scenario 
the discovery of a heavy Higgs boson in the parameter region Ma > 500 GeV will become 

5 The neglected higher-order corrections (e.g. from bosonic loops and initial state radiation) can in principle 
be sizable, see e.g. Fig. 3 in the second reference in [88]. However, they constitute corrections to the Born 
matrix element, i.e. they contain the WWH coupling and thus show the decoupling for the Ma values 
considered here. Therefore the neglected corrections would not substantially change the contours shown in 
this section. 
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possible also for large tan (3, tan (3 > 35. 

While without the inclusion of beam polarization we do not find a significant enhancement 
of the LC discovery reach as a consequence of the loop corrections in the four benchmark sce- 
narios analyzed in Fig. 14.151 this behavior changes in other regions of the MSSM parameter 
space. As a particular example, we investigate the M^-tan (3 plane in the "cxf^ h " ( "enhanced 
cross section" ) scenario, which is defined by 



Msusy = 350 GeV, /i = 1000 GeV, 



(4.46) 



with the other parameters as in the m™ ax scenario, eq. (J2J). The ajf^ scenario is thus char- 
acterized by a relatively small value of Mgusy and a relatively large value of \x. Large tan j3 
values, tan/3 > 30, can result in low and experimentally ruled-out b masses and are therefore 
omitted. 



'tan 




poLaxi'Ded 



iCt, lii, &ju itil aClH iilH un 1U *fifl UH £iu Iiu 6 LIU \n nu 

.^r^, ?,r A ?,r A 

Figure 4.16: The cross sections for e + e~ — > vv H are shown in the M^-tan/3 plane for the 
cr^ h scenario, eq. (|4.46|) . The tree- level cross section (left) including the finite wave-function 
corrections is compared to the a e g approximation (middle) and the one-loop corrected cross 
section (right column). The upper (lower) row shows the production cross section for unpo- 
larized (100% polarized) electron and positron beams. The color coding is as in Fig. 14.151 



In the upper row of Fig. 14. 161 we compare ajj(a e e), and a x H for the parameters of the 
a]f h scenario, eq. (|4.46|h in the unpolarized case. The figure shows that both the inclusion of 
the finite Higgs propagator wave-function corrections, see eqs. (I4.32|) - (j4.33|) . as compared to 
the a e fj approximation (left vs. middle column) and of the genuine one-loop corrections (right 
vs. left column) is very important in this scenario. According to the a e g Born approximation, 
the parameter area in which observation of H is possible would not be significantly larger 
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Figure 4.17: The tree-level cross section for e + e~ — > vv H, er#(a: e ff), using an a c g- approxi- 
mation (upper), and the one- loop cross section, a l H (lower row), are shown in the /x-Msusy 
plane for the parameters of eq. (J2J) and Ma = 600 GeV and tan/3 = 4. The left column 
shows the unpolarized case, while in the right column the effects of beam polarization are 
included. The color coding is as in Fig. 14.151 
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Figure 4.18: The tree-level cross section for e + e~ — > vv H, o#(a: e ff), using an a e fr approxi- 
mation (upper), and the one-loop cross section, a\j (lower row), is shown in the /x-Msusy 
plane for the parameters of eq. (0) and Ma = 700 GeV and tan/? = 4. The left column 
shows the unpolarized case, while in the right column the effects of beam polarization are 
included. The color coding is as in Fig. 14.151 
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than in the benchmark scenarios discussed in Fig. 14.151 For Ma ^ 500 GeV, observation of 
the H boson is only possible for either rather small, tan/3 < 5, or rather large, tan/3 > 28, 
values of tan/3. Inclusion of the finite Higgs propagator wave-function corrections, which 
ensure the correct on-shell properties of the outgoing Higgs boson, changes the situation 
considerably. While for small and large values of tan/3 the additional corrections suppress 
the e + e~ —>■ vvH cross section, restricting the H observability to values of Ma below 
500 GeV, observation of the heavy CP-even Higgs boson of the MSSM becomes possible 
for Ma < 550 GeV for a significant range of intermediate values of tan/3. Including also 
the genuine one-loop corrections (right column) leads to a further drastic enhancement of 
the parameter space in which the H boson could be observed. The observation of the H 
boson will be possible in this scenario for all values of tan/3 if Ma ^ 600 GeV, i.e. the 
discovery reach of the LC in this case is enhanced by about 100 GeV compared to the HA 
pair production channel. 

The prospects in this scenario for observation of the heavy CP-even Higgs boson of the 
MSSM become even more favorable if polarized beams are used. The lower row of Fig. 14.161 
shows the situation with 100% polarization of both beams for the ajf 1 scenario, eq. ()4.46|) . 
While in this case the tree-level result (both for the case including the finite wave-function 
corrections and for the a c s approximation) gives rise to observable rates for Ma ^ 600 GeV 
for a certain range of tan/3 values only, the further genuine loop corrections enhance the 
cross section significantly. In this situation the observation of the heavy CP-even Higgs 
boson might be possible for values of Ma up to about 700-750 GeV for all tan/3 values, 
corresponding to an enhancement of the LC reach by more than 200 GeV. Cross-section 
values in excess of 0.05 fb are obtained in this example for all tan (3 values for Ma ^5 600 GeV. 

In order to investigate whether this result is a consequence of a very special choice of 
SUSY parameters or a more general feature, in Fig. 14.171 we choose the parameters of the 
m max scenar j f or a fixed combination of Ma and tan /3, Ma = 600 GeV, tan f3 = 4, but 
scan over /i and Msusy- The choice M# ~ Ma = 600 GeV implies that the HA production 
channel is clearly beyond the reach of a 1 TeV LC. The upper row of Fig. 14.171 shows that 
according to the tree-level cross section (using the a e ff approximation) an observable rate for 
a heavy CP-even Higgs boson with Ma = 600 GeV cannot be found for any of the scanned 
values of /i and Msusy- Inclusion of the further loop corrections changes this situation sig- 
nificantly and gives rise to observable rates in this example for nearly all Msusy 500 GeV 
if [i > 500 GeV. The visible "structure" at M S usy ~ 500 GeV is the result of several 
competing effects that affect the finite Higgs wave-function corrections. 

The same analysis, but with 100% polarization of both beams, is shown in the right 
column of Fig. 14.171 The tree-level a e g approximation results in an observable rate in nearly 
the whole plane apart from the area with fi > 1000 GeV and Msusy ^ 350 GeV. Adding 
the loop corrections again improves the situation. No unobservable holes remain in the \i- 
Msusy plane, i.e. the heavy CP-even Higgs boson with M# « 600 GeV should be visible 
at a 1 TeV LC with (idealized) beam polarization in this scenario. The production cross 
section is larger than 0.02 fb for all /i and Msusy ~ 500 GeV and mostly even larger than 
0.05 fb for fx > 500 GeV. 

The same analysis, but for an even larger Higgs boson mass scale, Ma = 700 GeV is 
shown in Fig. 14.181 For the a c s approximation the production cross section is lower than 
0.01 fb in the whole /i-M S usY-plane. However, for the loop corrected cross section, especially 



85 



including polarization, the heavy CP-even Higgs boson could be observable for all \i values 
for Msusy < 500 GeV. 

In Fig. 14.191 we compare for all five scenarios analyzed in this section the results for the 
cross section when different parts of the generic one- loop corrections are taken into account. 
In all cases the full result for the Higgs propagator corrections being absorbed in the lowest- 
order cross section are employed. The left column shows the result containing the corrections 
from all fermion and sfermion loops, ajj, which is repeated from previous plots for comparison 
purposes. The middle column shows the cross-section prediction based on taking into account 
only the corrections from the third-generation fermions and sfermions, which have turned 
out to be the leading corrections for the light Higgs-boson production. The result shown in 
the right column have been obtained by including only the purely sfermionic contributions 
from all generations, i.e. the fermion-loop corrections are omitted in this case. As expected 
from Fig. 14.151 m the four scenarios defined in the Appendix, where the loop corrections 
turned out to modify the parameter regions in which H observation becomes possible only 
slightly, omitting the contributions of the first two generations of fermions and sfermions 
and of the fermion loops of all three generations does not lead to a qualitative change in the 
H discovery reach. In the a'jf 1 scenario, on the other hand, the genuine one-loop corrections 
had a considerable impact on the area in the M^-tan/5 plane in which H observation is 
possible, see Fig. 14.161 The result displayed in the middle column for this scenario shows 
that the bulk of the corrections comes from the third generation of fermions and sfermions, 
i.e. omission of the first two generations does not lead to significant effects in the M^-tan (3 
plane. The result in the right column for this scenario shows furthermore that the omission 
of all fermion-loop corrections leads only to very moderate changes of the parameter regions 
where H observation is possible. As a consequence, the by far dominant corrections in this 
scenario can be identified as the ones from the sfermions of the third generation. This is 
contrary to the h production, where we found that the fermionic corrections are mostly larger 
than the sfermionic ones. 

4.3 Single charged Higgs boson production 

In this section the dominant contributions to the process e + e~ — > ev e H ± [95] are briefly 
described. Since there is no {7, Z}W ± H T coupling, the single charged-Higgs production 
starts at the one-loop level. As for e + e~ — > uu{h,H}, we take into account the leading 
corrections arising from the full set of SM fermion and sfermion loops. Bremsstrahlung 
diagrams are suppressed by an additional 0(a) as compared to the leading contribution and 
have thus been omitted. 

4.3.1 Calculation 

The loop corrections that enter the process e + e~ — > ev e H ± at the one-loop level are generi- 
cally depicted in Figs. 14.201 14.211 The contributions involve all corrections from fermion and 
sfermion loops (which give contributions only to self-energies and vertices). Contributions 
oc m e /M w were neglected. 

Furthermore, counterterm contributions are needed for the W ± H T self-energy correc- 
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Figure 4.19: a^- is shown in the M^-tan/3 plane for the four benchmark scenarios defined 
in eqs. (0)-© and the scenario for y/s = 1 TeV in the unpolarized case. The one- 
loop result containing the corrections from all fermion and sfermion loops (left) is compared 
with the result from only third-generation fermions and sfermions (middle) and the purely 
sfermionic corrections from all families (right column). The color coding is as in Fig. 14.151 
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tions, see Ref. [94]. In order to generate the counterterms it is sufficient to introduce the 
field renormalization for the H ± — W ± mixing, SZhw- This yields the Feynman rules: 



rcT^w^P)] 



mLsz hw 



rcrbuWflF] = -ieM w g^5Z HW , 
T CT [Z,W±H^} = teM w ^g^5Z HW 



In the on-shell scheme SZhw is given by 



5Z 



HW 



M w 



ReZ HW (M 2 H± ] 



(4.47) 
(4.48) 
(4.49) 

(4.50) 



4.3.2 Results 



The results for H + and H~ production are the same if CP is not violated (which we assume 
in this section). In Fig. 14.221 we show the typical size of the production cross section for 
e + e~ — > e + v e H~ for unpolarized external particles. The parameters are chosen according to 
the four benchmark scenarios described in the Appendix, with Ma = 250 GeV and tan (3 = 2 
and 10 (with Mn± « 262 GeV). Concerning the discovery of the charged Higgs boson, the 
cross section has to be doubled, due to the production of both, H + and H . 

In Fig. I4.22l the cross section for e + e _ — > e + u e H~ is shown as a function 
of y/s. The rise of the cross section at \/s ~ M H ± + M w is due to the 
W propagator in the type of diagram on the right becoming resonant. 
The resonance was treated with a fixed W width. No cut on the electron 
emission angle had to be put, since the result is finite in the limit of forward electron 
scattering. 
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Figure 4.22: The e + e~ — > e + v e H~ production cross section as a function of \fs. 



The variation within the four benchmark scenarios is small. For tan (3 = 10 the charged- 
Higgs production cross section stays at a negligible level. Even for tan (3 = 2 it stays below 
0.01 fb for 2 Mh± w \fs < 500 GeV. Using polarized e + and e~ beams, the cross section 
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could be enhanced by about a factor of 2. In summary, however, the single charged-Higgs 
production, e + e~ — > eu e H ± , could increase the potential of a LC for the detection of the 
heavy MSSM Higgs-boson spectrum only for parameters beyond the typical benchmark 
scenarios. 
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Chapter 5 

Higgs boson decays to SM fermions 



After the detection of a scalar particle it is mandatory as a next step to measure its couplings 
to gauge bosons and fermions and also its self-couplings very accurately, in order to establish 
the Higgs mechanism and the Yukawa interactions experimentally. The determination of the 
trilinear Higgs-boson self-couplings might be possible at a future linear e + e~ collider with 
high luminosity [41]. 

In this section we concentrate on the Feynman-diagrammatic corrections (i.e. beyond the 
effective coupling approximation) to the coupling of the lightest MSSM Higgs boson to SM 
fermions. Since the b-, the c-quark and the r-lepton are the heaviest particles for which the 
decay h — *■ ff is kinematically allowed, it is of particular interest to calculate the correspond- 
ing decay rates and branching ratios with high precision [11,82,99,106-110]. We analyze 
these decay rates and branching ratios, taking into account the Higgs-boson propagator cor- 
rections, where at the one- loop level the full momentum dependence is kept [99]. We also 
take into account the one-loop vertex corrections resulting from gluon, gluino, and photon 
exchange together with real gluon and photon emission as given in Ref. [11]. Only the purely 
weak 0(a) vertex corrections have been neglected (see e.g. Ref. [Ill] for a recent evalua- 
tion). We numerically investigate the effect of the two-loop propagator contributions and 
the one-loop gluino-exchange vertex correction. The latter is supplemented with a resumma- 
tion of the leading contributions of 0((a s tan/5) n ), see Sect. I2.3.H Some phenomenological 
consequences are investigated. 

We show analytically that the Higgs-boson propagator correction with neglected momen- 
tum dependence can be absorbed into the tree-level coupling using the effective mixing angle 
from the neutral CP-even Higgs boson sector [99]. The result in this approximation is then 
compared with the full result. 



5.1.1 Calculation of the decay amplitude 

Our main emphasis in this section is on the fermionic decays of the light Higgs boson, but 
for completeness we list the expressions for both h and H. The amplitudes for the decays 
h,H—*ff can be written as follows: 



5.1 Higher-order corrections 
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A(h^ff) = yJZ h {T h + Z hH T H j , (5.1) 
A(H^ff) = \[z~ H (r H + Z Hh T^) , (5.2) 

with Z h , Z H , Z hH and Z Hh given in eqs. Ij4.34|)-fj4.37|). 

5.1.2 The a e s approximation for decay amplitudes 

The dominant contributions for the Higgs boson self-energies can be obtained by setting 
a 2 = 0. see Sect. l2~TO 

t(q 2 ) -> S(0) = £ . (5.3) 
With the approximation (|5.3|) (see also eqs. (|2.48|) . ()2.49p ) one deduces 

£hH 

ZhH = = — tan Act, (5.4) 

Zuh = * — = + tanAa, (5.5) 

M H -m\ + E hh 



and Zh can be expressed as 



1 + 



Analogously one obtains 



— = cos 2 A a . (5.6) 

1 + tan 2 Aa 

Z H = cos 2 Aa . (5.7) 
At tree level, the vertex functions can be written as 

for d-type fermions, (5.8) 



t i iemr sin a s^(d) 

1 h — o — id w — L f sin a 

In— o — TT w = f cos a 

n 2s v ,M w cos/3 / 

m — ierrif cos a s~i(u) 

1 h = o — IT — = C)' cos a 
1// = ^ — tj- — = C ) sin a 



for u-type fermions . (5.9) 
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Incorporating them into the decay amplitude yields the corresponding a e g approximation: 
A eS (h^ff) = ^Z~ h (T h + Z hH T H ) 



= Cf cos Aa (sin a — tan Act (— cos a )) 

= Cf sin(ct + Aa) 

= sina e ff (for d- type fermions), (5.10) 

Af&(h—>ff) = Cf cos a e ff (for u-type fermions), (5.11) 

Acs(H —>■ ff) = — cos ct eff (for d- type fermions), (5-12) 

A e g(H —> ff) = sin a e fr (for u-type fermions) . (5.13) 

5.1.3 Evaluation of the Higgs-boson propagator corrections 

The Higgs boson self-energies employed in the field renormalization constants, see eqs. ()4.34j) - 
()4.37|) . have been evaluated in the Feynman-diagrammatic approach according to 

£ s (q 2 ) = ££\q 2 ) + £V (o) , s = h,H,hH, (5.14) 

where the momentum dependence has been neglected only at the two-loop level, while the 
full momentum dependence is kept in the one-loop contributions, see Sect. 12.3.41 

For the numerical evaluation in a first step of approximation for the calculation of the 
decay width T(h — > //) the momentum dependence is neglected everywhere in the Higgs 
boson self-energies, the q 2 = approximation (see eq. (|5.3Jl ): 

S s (g 2 )^S«(0) + Si 2) (0) , s = h,H,hH. (5.15) 

This corresponds to the a e ff-approximation, as described in Sect. 15.1.21 

In a second step of approximation we approximate the Higgs boson self-energies by the 
compact analytical formulas given in Ref. [17]: 

± s (q 2 ) = S« approx (0) + s( 2 )™ x (o) , s = h,H,hH . (5.16) 

Here the full result of the self-energy corrections is approximated by an expansion in terms 
of m t /MsusY and Xj/Msusy; yielding relatively short expressions which allow a very fast 
numerical evaluation. In the following, this approximation is labeled by a e fr(approx). 

5.1.4 Decay width of the lightest Higgs boson 

At the tree level, the decay width for h — > // is given by 

r„(^//) = iv c ^(i-i^y irj". (5.17) 

The electroweak propagator corrections are incorporated by using the higher-order decay 
amplitude (|5.1j) 

3 

M f 4 w? \ ~ z 

r 1 = r l (h^ff) = N C j±\l- 1 ^) \A(h^ff)\ 2 - (5.18) 
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The a eff -approximation is given by 




4m 2 f 



3 
2 



ri, eff = r lfiS (h -> //) = iv c -^ 1 



A eS (h^ff)\ 2 . 



(5.19) 



In this section we consider only those electroweak higher-order contributions that enter 
via the Higgs boson self-energies. These corrections contain the Yukawa contributions of 
O {G prnf I 'Myy) , which are the dominant electroweak one-loop corrections to the Higgs-boson 
decay width, and the corresponding dominant two- loop corrections, see Sect. 15.1.31 The pure 
weak 0{a) vertex corrections are neglected (they have been calculated in Refs. [11,111]. 
They were found to be at the level of only a few % for most parts of the MSSM parameter 
space, see also Ref. [99]). 

QED corrections 

Here we follow the results given in Refs. [11,106-108]. The IR-divergent virtual photon 
contribution is taken into account in combination with real-photon bremsstrahlung yielding 
the QED corrections. The contribution to the decay width induced by 7-exchange and 
final-state photon radiation can be cast into the very compact formula 



QCD corrections: gluino contributions 

In this section we will first focus on the decay <f> — > bb,(j) = h, H only, since for this decay 
channel the gluino contributions are especially relevant. Concerning the evaluation of the 
corresponding gluino-exchange Feynman diagrams for — > bb, we follow the calculation given 
in Ref. [11] (see also Ref. [99]), similar results can also be found in Ref. [109]. The additional 
contributions to the decay amplitude induced by gluino-exchange are incorporated, following 
the effective Lagrangian eq. (|2.54|) : 



Ar^ = rf • <p = h } H , 

where for m 2 f <C M? the factor 8Tt, has the simple form 



(5.20) 




(5.21) 




(5.22) 




(5.23) 
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for real Z hH) Z Hh (i.e. neglecting the imaginary part in eqs. (|5.22|) . (j5.23|l ): r~ and T? are 
given by 



9 



-,11 



AT> 



r 



H 



+ SUmg) + Am b - 2m 2 b + ^vA m l)) 



(5.24) 



. (5.25) 



AT- ' denote the gluino vertex-corrections, whereas S b represents the gluino contribution 
to the bottom self-energy corrections. Explicit expressions for these terms can be found 
in Ref. [11]. 

For large values of tan/3 in combination with large values of the gluino-exchange 
corrections to A{<p —>■ bb) can become very large. They are resummed to all orders of 
(a s tan/3) n via the inclusion of Am& into A^ H (eqs. (|5.22|) . (|5.23j) V The results for A^ ,H 
constitutes the currently best available evaluation of the gluino-exchange corrections to — * 
bb. It contains the resummation of potentially large corrections ~ (a s tan/3) n to all orders 
as well as the full evaluation of the g — b vertex corrections including momentum effects. 



In the case of 
resummed part: 



Ah 

A g 



H 



A 1 ' 

A 9 



"p/l 

a 



1 + 
1 + 



f f ) f 7^ b the expressions are similar, but do not contain the 
1 



r\ + z hH T H 
i 

r_ff + z Hh r h 



AT^ 



q *=Ml 



Re [r h g + z hH rf], 

Re [If + Z Hh r h ] , 
+ E^(m^) - 2ml (^U m l) + ^U m D) 



r 



H 



r 



ATi 



1 2 =M 2 h 



+ S* (mg) - 2m 2 b (4^K) + ^v,M)) 



(5.26) 
(5.27) 

(5.28) 
(5.29) 



QCD corrections: gluon contributions 

The corresponding results have been obtained in Refs. [11,82,106-110]. The additional 
contribution to the decay width induced by gluon exchange and final-state gluon radiation 
can be incorporated into (|5.22j) by 



g<j 



A g 



m q {Ml 



m n 



1 + 



2vr 



r 9 + 8 



(5.30) 
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The running quark mass m g (M?) is calculated via 

c(q 2 ) 



m q (q 2 ) 
<q 2 ) 



m 



q c(m q ) 



2tt 



2\ \ -7o/2/3o 



l CgiTo ~ A)7i) a ^(g 2 ) 



A 2 



; 7T 



(5.31) 
(5.32) 



where m q is the pole mass and m q {m 2 q ) = m q . The coefficients in eq. (|5.32j) are: 

33 - 2N f 
Po - § , 

Pi = m-^.N f , 



7o 
7i 



404 40 Ar 

1 N f 

3 9 i 



(5.33) 



where Nf — 5 for / = c, 6, which are considered here. The strong coupling constant a s is 
given up to two loops by: 



a s (q 2 ) 



Air 
PoL q 



(5.34) 



where L q = \og(q 2 / Aq CD ) . (For the numerical evaluation Aqcd = 220 MeV has been used.) 
Numerically, more than 80% of the gluon-exchange contribution is absorbed into the running 
quark mass. 



Decay width and branching ratio 

The result for the decay width of h — > ff, including strong corrections, is given by 



M h A Am 2 
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99 



8 7T 



N 



M. 



H 



c- 



i7T 



Ml 



4mf\ 2 



Ml 



Including also the QED corrections, the full decay width is given by 



r(A-&&) = r^ + Ar;, 

T(H^bb) = Tg + AT*. 



(5.35) 
(5.36) 



(5.37) 
(5.38) 



Summing over f = b,c,T and adding T(h — > gg) (which can be numerically relevant [112]), 
results in an approximation for the total decay width 



r tot = Yl r(h->ff) + r(h 

f=b,c,r 



(5.39) 
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We do not take into account the decay h — > AA (see e.g. Ref. [113] for a detailed study). 
Although it is dominant whenever it is kinematically allowed, it plays a role only for very 
small values of tan/3 (tan/3 < 1.5) which will not be considered here because of the limits 
obtained at LEP2 [6]. We also assume that all other SUSY particles are too heavy to 
allow further decay channels. In addition, we neglect the decay h — > WW* which can 
become substantial for > 120 GeV. The quantitative change in our results due to this 
approximation is small, the qualitative change is negligible. 
The fermionic branching ratio is defined by 

Rf = B R(h -> //) = r(/ \T — • (5-40) 

1 tot 

The results of this section, including all decay modes, is incorporated into the Fortran code 
FeynHiggs2.1. In a comparison with the code Hdecay [114] we found good agreement. 

5.2 Phenomenological implications 

Concerning the numerical evaluation of the Higgs-boson propagator corrections, we follow 
Sect. 15 .Ol For tan/3 we mostly concentrate on two representative values, a relatively low 
value, tan/3 = 3, and a high value, tan/3 = 40. For sake of comparison we also consider an 
intermediate value of tan /3 = 20 in some cases. If not indicated differently, the other MSSM 
parameters are chosen as follows: \i = —100 GeV, M 2 = M SUSY = rrig, rrig = 500 GeV, 
A^ = A t . With "no mixing" we denote the case X t = 0, whereas "maximal mixing" denotes 
X t = 2Msusy- The other SM parameters are given in the Appendix. The mass Ma of 
the CP-odd Higgs boson is treated as an input parameter and is varied in the interval 
50 GeV < Ma < 500 GeV. In some cases, also tan/3 is varied in the interval 2 < tan/? < 50. 

The corresponding values for follow from eq. (j2.41|) . M^, derived in this way, subse- 
quently enters the numerical evaluation of the formulas presented in section 15.11 Thus the 
variation of in the plots stems from the variation of Ma (or tan/?) in the above given 
range. 

5.2.1 Effects of the two- loop Higgs-propagator corrections 

We first focus on the effects of the two-loop Higgs-boson propagator corrections in comparison 
with the one-loop case. They have been evaluated at the one- and at the two-loop level as 
described in Sect. 15.1.31 Fi gure ISTl shows the results for T(h — > bb) for a common scalar quark 
mass rriq = Msusy = 1000 GeV and tan/3 = 3 and tan/3 = 40 in the no-mixing and the 
maximal-mixing scenario. The QED and the QCD gluon and gluino vertex contributions 1 are 
also included. 

In the small tan/? scenario, larger values for T(h — > bb) are obtained for maximal mixing. 
The two-loop corrections strongly reduce the decay width. In the large tan/3 scenario the 

lr The gluino vertex corrections are included without the resummation formula, since either the small 
value of \/i\ or of tan/3 results in very small values of Araj. In Sect. I5.2T^1 however, also the full correction 
including resummation is taken into account. 
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Figure 5.1: T(h — > bb) is shown as a function of M^. The Higgs-propagator corrections 
have been evaluated at the one- and at the two-loop level. The QED, gluon and gluino 
contributions are included. The other parameters are /i = —100 GeV, M 2 = M SU sy, m g = 
500 GeV, Ab = A t , tan/3 = 3,40. The result is given in the no-mixing and maximal-mixing 
scenario. 

variation is mainly a kinematical effect from the different values of Mh at the one- and two- 
loop level. The absolute values obtained for F(h — > bb) are three orders of magnitude higher 
in the tan /3 = 40 scenario, which is due to the fact that T(h — > bb) ~ 1/ cos 2 (3 . 

In Fig. 15.21 the three decay rates T[h — > bb), T(h — > t + t~) and T(h — > cc) are shown as a 
function of Mh- The results are given in the no-mixing scenario for Msusy = 500 GeV and 
tan (3 = 3,40. 

In the low tan/? scenario T(h —>■ bb) and T(h — > t + t~) are lowered at the two-loop level, 
while T(h — > cc) is increased. The decay rate for h — > bb is about one and two orders of 
magnitude larger compared to the ones of h — » t + t~ and h — > cc, respectively. In the large 
tan/? scenario the shifts are again dominated by the kinematical effect from the different 
values of at the one- and two-loop level. In the maximal-mixing case, which is not 
plotted here, we find qualitatively the same behavior. 

We now turn to the effects of the two-loop corrections to the branching ratios. For not too 
large values of Mh, T tot is strongly dominated by T(h — > bb). For large values of Mh the decay 
into gluons becomes more relevant. In Fig. 15.31 we show the branching ratio BR(/i — > bb) as 
a function of Mh and M4. For values of Ma ^ 250 GeV there is a non- negligible difference 
between one-loop and two-loop order, where at the two-loop level the branching ratio is 
slightly enhanced. Compared in terms of Mh there is nearly no change for small values of 
Mh in the low and in the high tan f3 case. Here BR(/i — > bb) is changed by less than about 
1%, see Fig. 15.31 BR(/i — > r + r~) is increased by less than about 2%. BR(/i — > cc) can be 
increased at the two-loop level by O(50%), but remains numerically relatively small. For 
tan/3 = 40 the main difference arises at the endpoints of the spectrum, again due to the 
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Figure 5.2: T(h — ► 66), r(/i — > t + t~) and r(/i — ► cc) are shown as a function of M^. The 
Higgs-propagator corrections have been evaluated at the one- and at the two-loop level. The 
QED, gluon and gluino contributions are included. The other parameters are \i = —100 GeV, 
M 2 = Msusy, fn-g = 500 GeV, Ab = A t , tan/5 = 3,40. The result is given in the no-mixing 
scenario. 

fact that different Higgs boson masses can be obtained at the one-loop and at the two-loop 
level. For tan/? = 3, however, also several GeV below the kinematical endpoints there is a 
sizable effect on BR(/i — * 66). Thus, in the experimentally allowed region of Mh, the two-loop 
corrections can have an important effect on BR(/i — > 66). 

5.2.2 Vanishing decay rate for h — > bb 

The search for the Higgs boson, especially at e + e~ colliders, often relies on b tagging, since on 
the one hand the lightest CP-even Higgs boson decays dominantly into 66, and on the other 
hand 6 tagging can be performed with high efficiency For some combinations of parameters, 
however, T(h —>■ bb) can become very small and thus BR(/i — > bb) can approach zero as a 
consequence of large Higgs-boson propagator corrections or large gluino vertex-corrections, 
making Higgs boson search possibly very difficult for these parameters. Higgs boson searches 
then have to rely on flavor-independent decay modes. In order to have reliable predictions 
for these regions of parameter space a full calculation of the one-loop hbb vertex corrections, 
including all 0(a) contributions, is necessary [111]. 

Effects of two-loop propagator corrections 

We first demonstrate the effect of the two-loop propagator corrections on the values of the 
parameters, especially of Ma, for which BR(/i — > 66) goes to zero. We also show the impact 
of the inclusion of the momentum dependence of the Higgs boson self-energies that is often 
neglected in phenomenological analyses of the decays of the lightest CP-even Higgs boson. 
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Figure 5.3: BR(/i — > bb) is shown as a function of Ma and for the same settings as in 
Fig. I5.ll The QED, gluon and gluino contributions are included. 



In Fig. I5.4I BR(/i — > bb) is shown as a function of Ma- The Higgs boson self-energies are 
evaluated at the one-loop and at the two-loop level with and without momentum dependence 
(see eq. ()5.3|) ). The other parameters are tan/3 = 25, Msusy = 500 GeV, m~ g = 400 GeV, 
M 2 = 400 GeV, X t = 400 GeV, A b = A u ft = -1000 GeV. The inclusion of the two- 
loop propagator corrections shifts the Ma value for which BR(/i — > bb) becomes very small 
by about —35 GeV. The inclusion of the momentum dependence of the Higgs boson self- 
energies induces another shift of about —6 GeV. In order to have reliable phenomenological 
predictions for the problematic Ma values the two-loop corrections as well as the inclusion 
of the momentum dependence is necessary. Note that the inclusion of the gluino vertex 
corrections as well as the purely weak vertex corrections can also have a large impact on the 
critical Ma values. 
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Figure 5.4: BR(/i — > 66) is shown as a function of Ma- The Higgs boson self-energies are 
evaluated at the one-loop and at the two-loop level with and without momentum dependence 
(see eq. (|5.3j) ). The other parameters are tan/3 = 25, M SU sy = 500 GeV, rrig = 400 GeV, 
M 2 = 400 GeV, X t = 400 GeV, A b = A t , fj, = -1000 GeV. 



Effects of 0(a s a b ) 

We now turn to the effects of the additional contributions induced by 0{a s a b ) corrections 
and the corresponding resummation as discussed in Sect. 12.3.41 We recall the tree- level 
couplings of the lightest CP-even Higgs boson to the up-type and down-type SM fermions: 

r „ _ iem u - cos a ^ d _ i e m d sin a 

h ~2s w M w sin/? ' h ~ 2s w M w cos/3 ' [ ' ' 

The main source (besides the SM QCD corrections and the Higgs boson propagator cor- 
rections) modifying the hbb and hr + r~ vertices are the corrections modifying the relation 
between the bottom quark or r lepton mass and the corresponding Yukawa couplings. 

The potentially large corrections to the couplings come from the tan /3-enhanced 
threshold effects in the relation between the down-type fermion mass and the corresponding 
Yukawa coupling [43], already mentioned in Sect. 15.1131 A simple way to take into account 
these effects is to employ the effective Lagrangian formalism of Refs. [44,55], see Sect. 12.3.41 
where the coupling of the lightest Higgs boson to down-type fermions (expressed through 
the fermion mass) can be written as 



d , iem d sinews 



cfr 2s w M w cos 



1 "T~ — (1 + c °t a eS COt (3) 

1 + Am d 



(5.42) 



where Am^ contains (as described in Sect. I2.3~i|) the tan /3-enhanced terms, and other sub- 
leading (i.e. non tan /3-enhanced) corrections have been omitted. In the case of the coupling 
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to the bottom quarks, the leading contributions to Am b are of 0(a s ) and 0(a t ), coming 
from diagrams with sbottom-gluino and stop-chargino loops, respectively. In the case of the 
r leptons, the leading contributions are of 0(a T ), coming from sneutrino-chargino loops. 
The terms containing Am^ in eq. ()5.42j) may be relevant for large tan /? and moderate values 
of Ma- When Ma is much bigger than Mz, the product cota e fj cot (3 tends to —1, and the 
SM limit is again recovered. 

The effects of the higher-order corrections to the couplings of the lightest Higgs boson 
to the down- type fermions appear very pronounced in the "small a e fj" scenario [53], see the 
Appendix. 

In the upper row of Fig. 15.51 we show the ratio sin 2 a e g/ cos 2 /? in the M^-tan/3 plane, 
evaluated in the "small ct e fr" scenario. As explained above, replacing a by a e g in the tree- 
level couplings of eq. (j5.41j) gives rise to regions in the M4-tan (3 plane where the effective 
coupling of the lightest Higgs boson to the down-type fermions is significantly suppressed 
with respect to the Standard Model. The region of significant suppression of sin 2 a e fr/ cos 2 (3 
as evaluated without the inclusion of the sbottom corrections beyond one-loop order, is shown 
in the upper left plot. The upper right plot shows the corresponding result as evaluated 
with FeynHiggs2.1, where the 0(ai,a s ) and 0(ab(a s tan/?)") corrections as well as the new 
O (a 2 ) ones are included. The new corrections are seen to have a drastic impact on the 
region where sin 2 a e ff/ cos 2 /? is small. While without the new corrections a suppression of 
70% or more occurs only in a small area of the M^-tan /3-plane for 20 < tan/3 < 40 and 
100 GeV < Ma 200 GeV, the region where sin 2 a e fj/ cos 2 /? is very small becomes much 
larger once these corrections are included. It now reaches from tan (3 > 15 to tan (3 > 50, and 
from Ma ^ 100 GeV to Ma ^ 350 GeV. The main reason for the change is that the one- 
loop 0(ab) corrections to the Higgs-boson mass matrix, which for large tan/3 would prevent 
a c ff from going to zero, are heavily suppressed by the resummation of the 0(ab(a s tan/3) n ) 
corrections in the bottom Yukawa coupling. This kind of suppression depends strongly on 
the chosen MSSM parameters, and especially on the sign of fi. 

In order to interpret the physical impact of the effective coupling shown in the upper row 
of Fig. 15.51 the Am^ terms in eq. (J5.42)) as well as further genuine loop corrections occurring 
in the h — > bb process have to be taken into account as described in Sect. 15.1.31 The effect 
of these contributions can be seen from the plot in the lower row of Fig. 15.51 where the ratio 
T(h — > 66)mssm / r(/t — > bb)sM is shown, which has been evaluated by including all terms of 
eq. (|5.42|) as well as all further corrections given in eq. (|5.26|) . The region where the partial 
width T(h — > bb) within the MSSM is suppressed compared to its SM value is seen to be 
somewhat reduced and shifted towards smaller values of Ma as compared to the region where 
sin 2 a e s / cos 2 (3 is small. 

5.2.3 Effects of the gluino vertex corrections 

In this subsection we focus on the effect of the gluino-exchange contribution to the hff 
vertex corrections in a more general way, i.e. for parameters where hff does not go to zero. 

Fig. 15.61 shows T(h — > bb) in three steps of accuracy: the dotted curves contain only the 
pure self-energy correction, the dashed curves contain in addition the QED and the gluon- 
exchange correction. The solid curves show the full results including also the gluino-exchange 
correction. The results are shown for X t = (no mixing), /z = —100 GeV, rrig = 500 GeV, 
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Figure 5.5: Regions of significant suppression of the coupling of h to down-type fermions 
in the M^— tan /3-plane within the "small a e fr" benchmark scenario. The upper left plot 
shows the ratio sin 2 a c fr/ cos 2 (3 as evaluated with FeynHiggsl.O (i.e. without 0(a.bOi s ) and 
0(ab(a s tan/3) n ) corrections), while the upper right plot shows the same quantity as evalu- 
ated with FeynHiggsl.3 (i.e. including these corrections). The plot in the lower row shows 
the ratio T(h &&)mssm / r(/i — > bb) SM , i.e. the partial width for h ^ bb normalized to its 
SM value, where the other term in eq. (|5.42j) and further genuine loop corrections are taken 
into account (see text). 
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Figure 5.6: T(h — ► bb) is shown as a function of for three values of tan/3. The Higgs- 
propagator corrections have been evaluated at the two-loop level in the no-mixing scenario. 
The dotted curves shows the results containing only the pure self-energy corrections. The 
results given in the dashed curves in addition contain the QED correction and the gluon- 
exchange contribution. The solid curves show the full result, including also the gluino cor- 
rection. The other parameters are \x = —100 GeV, M2 = trig, Ab = A t . 




Figure 5.7: T(h —>■ bb) is shown as a function of fi for tan/3 = 20 and two different 
values of m- g in the no-mixing and the maximal- mixing scenario. The other parameters are 
rriq = 500 GeV, M 2 = 500 GeV, A b = A u M A = 100 GeV. 
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rriq = 200, 1000 GeV, in the left and right part of Fig. 15.61 respectively. For tan (3, three 
values have been chosen: tan/3 = 3,20,40. 

The left plot of Fig. 15.61 corresponds to a small soft SUSY-breaking scale, rrig = 200 GeV. 
The effect of the gluon contribution is large and negative, the effect of the QED correction 
is small. For this combination of rrig, rriq and fi the effect of the gluino correction is large 
and positive, as can be seen from the transition from the dashed to the solid curves. For 
tan/3 = 40 it nearly compensates the gluon effect, for tan/3 = 20 it amounts up to 20% 
of the gluonic correction, while for tan f3 = 3 the gluino-exchange contribution is negligible. 
Note that we have chosen a relatively small value of /i, u, = —100 GeV. For larger values of 
\fj,\ even larger corrections can be obtained. Hence neglecting the gluino-exchange correction 
in the large tan/3 scenario can lead to results which deviate by 50% from the full 0(a s ) 
calculation (see also Sect. 15.1. "4^1 . The right plot of Fig. 15.61 corresponds to rriq = 1000 GeV. 
The gluino-exchange effects are still visible, but much smaller than for rriq = 200 GeV. The 
same observation has already been made in Ref. [11]. For X t = 2rriq (maximal mixing) we 
find qualitatively the same behavior for the gluino-exchange corrections as for X t = 0. 

In Fig. 15.71 the pure gluino-exchange effect is shown as a function of \x. This effect in- 
creases with rising 2 nig and where for negative (positive) /i there is an enhancement 
(a decrease) in T(h —>■ bb). The size of the gluino-exchange contribution also depends on 
Ma, where larger effects correspond to smaller values of Ma, see also Ref. [11]. The small 
difference between the curves where the decay rate has been calculated without gluino con- 
tribution is due to the variation of Mh induced by different values of m g which enters at 
0(aa s ). Fig. 15.71 demonstrates again that neglecting the gluino contribution in the fermion 
decay rates can yield (strongly) misleading results. 

As a final remark, one should note that the gluino exchange contribution has only a 
relatively small impact on BR(/i — > bb). However, it can have a large influence on BR(/i — > 
t + t~) or BR(/i — > WW*). Both branching ratios are expected to be measurable at the per- 
cent level, see Tab. 11.31 While the Higgs-propagator contributions are universal corrections 
that affect T(h —>■ bb) and T(h — > t + t~, h — > H^W 7 *) in the same way (i.e. the influence on 
the effective coupling is the same in both cases), the gluino corrections, which influence only 
T(h bb), can lead to a different behavior of the decay widths. This has been analyzed in 
the mSUGRA, GMSB and AMSB scenario in Refs. [115, 116]. 

5.2.4 The (^-approximation 

Finally, we investigate the quality of the a e ff-approximation. In Fig. 15.81 we display the 
relative difference between the full result ()5.14|) and the a e s result, where the external mo- 
mentum of the Higgs self-energies has been neglected, see eq. ()5.15|) . The relative difference 
AT(h -> bb) = (T lull (h -»■ bb) - T a ^(h -> bb))/T lull (h -> 66) is shown as a function of M A 
for rriq = 1000 GeV and for three values of tan/3 in the no mixing and the maximal mixing 
scenario. Large deviations occur only in the region 100 GeV < Ma < 150 GeV, especially 
for large tan/3. In this region of parameter space the values of Mh and Mh are very close 
to each other. This results in a high sensitivity to small deviations in the Higgs boson 

2 This is correct for all values of trig considered in this work. A maximal effect is reached around trig ~ 
1500 GeV. The decoupling of the gluino takes place only for very large values, rrig > 5000 GeV. 
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self-energies entering the Higgs-boson mass matrix eq. ([2.45)1 . Another source of differences 
between the full and the approximate calculation is the threshold Ma = 2m t = 350 GeV in 
the one-loop contribution, originating from the top-loop diagram in the A self-energy, see 
Ref. [99] for more details. Here the deviation can amount up to 6%. This threshold is only 
present in the pure on-shell renormalization (see also Sect. I2.3T4")) . which has been used for 
these plots. Making use of the hybrid DR/on-shell renormalization (as incorporate in the 
latest version of FeynHiggs), these thresholds are absent. 

In Fig. 15.91 we compare the a e ff result (|5.15j) with the a e fr(approx) result ()5.16|) . where 
the Higgs boson self-energies have been approximated by the compact analytical expression 
obtained in Ref. [17]. Figure IS~9~1 displays the relative difference in the effective mixing angles, 
(sin a e s— sin a eff (approx))/ sin a e fr- Via eq. ()5.10|) sin a e fj directly determines the decay width 
T(h — > bb). The result is shown for rriq = 1000 GeV, for three values of tan f3 in the minimal 
and the maximal mixing scenario. Apart from the region around Ma ~ 120 GeV (compare 
Fig. 15.8)1 both effective angles agree better than 3% with each other. 

Concerning the comparison of the a e fj-approximations in terms of M^ (which is not 
plotted here), due to the neglected external momentum or the neglected subdominant one- 
or two-loop terms, Mh receives a slight shift. Besides this kinematical effect, the decay rate 
is approximated rather well for most of the Mh values: independently of rriq, the differences 
stay mostly below 2-4%, for the no-mixing case as well as for the maximal-mixing case. 
Only at the endpoints of the spectrum, due to the different Higgs-boson mass values, the 
difference is not negligible. 
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Figure 5.8: AT{h -> bb) = (T iul \h -> bb) - T a ^(h -> bb))/T iul \h -> bb) is shown as a 
function of Ma for three values of tan/3. The QED, gluon- and gluino-contributions are 
neglected here. The other parameters are \i = —100 GeV, M 2 = rriq, m- g = 500 GeV, 
Af, = A t , tan/? = 3,20,40. X t has been set to X t = (no mixing) or X t = 2rriq (maximal 
mixing) . 
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Chapter 6 
Conclusions 



We have reviewed the current status of higher-order corrections in the MSSM Higgs sector. 
This comprises the corrections to the Higgs boson masses in the real and complex MSSM, the 
existing corrections to Higgs boson production cross sections at the LC as well as corrections 
to Higgs boson decays to SM fermions. 

Concerning Higgs boson masses and couplings in the MSSM with real parameters, the 
higher-order corrections are known to be huge. The corrections described in this report lower 
the intrinsic error due to unknown higher-order corrections down to about ~ 3 GeV for the 
lightest Higgs boson mass. Likewise a precision at the per-cent level for the couplings can 
be reached. In the case of the MSSM with complex parameters, uncertainties of more than 
~ 5 GeV for the lightest Higgs boson mass have been eliminated. The effects of the full one- 
loop corrections have been shown to be substantial for all Higgs boson masses and couplings. 
We have investigated the higher-order contributions to Higgs boson production channels at 
the LC and found corrections of up to ~ 20% for the channel e + e~ — > hZ. In the WW fusion 
channel especially the non-decoupling behavior of the heavy CP-even Higgs boson has been 
investigated. Taking loop corrections into account the reach of the LC could be extended by 
up to ~ 200 GeV. Concerning the Higgs boson decays to SM fermions, we have shown that 
the corrections arising in the Feynman-diagrammatic approach are subtantial and have to be 
taken into account in an accurate analysis. This holds especially for gluino-loop corrections 
for the decay h — > bb. 

To summarize, in all investigated cases the higher-order corrections are substantial and 
must not be neglected in phenomenological analyses. Once a MSSM-like Higgs boson will 
have been discovered, and the corresponding quantities like masses, couplings and production 
cross sections will be measured at the LHC and (hopefully) furtheron at the LC, it will be 
mandatory to have even more precise evaluations at hand. In view of the complexity of 
these calculations in the past, and the time needed to perform them, further higher-order 
corrections will have to be evaluated continously in the future. Only then the required 
precision will be available at the same time as the experimental data will be determined 
from the collider experiments. 

At the same time it will be necessary not only to improve on the higher-order corrections 
but also on the reduction of the parametric uncertainties due to the imperfect knowledge of 
the input parameters. This is especially true for the top quark mass. Here certainly the LC 
precision will be necessary to match even the anticipated LHC precision in M^. 
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The tasks for theorist described above are huge. However, it still seems to be possible 
to achieve the required calculations, provided that there is enough financial and man power 
support from the whole community [117]. Then (but only then) it can be possible that 
theoretical and experimental high-energy physics has a bright future. 
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Appendix 



For our numerical results, the following values of the SM parameters are used (all other 
quark and lepton masses are negligible): 

Gp = 1.16639 x 1(T 5 , m T = 1.777 GeV, 

M w = 80.450 GeV, m t = 174.3 GeV, 

M z = 91.1875 GeV, m b = 4.25 GeV, ^ ' 

Y z = 2.4952 GeV, m c = 1.5 GeV. 

For our numerical evaluation we mostly rely on four benchmark scenarios that have been 
defined in Ref. [53] for Higgs boson searches at hadron colliders and beyond. The four 
benchmark scenarios are (more details can be found in Ref. [53]) 

• the "m™ ax " scenario, which yields a maximum value of M h for given Ma and tan (3, 

m t = 174.3 GeV, M SU sy = 1 TeV, fi = 200 GeV, M 2 = 200 GeV, 
Xi = 2M SUS Y, A T = A b = A t , m g = 0.8M SVSY , 

• the "no-mixing" scenario, with no mixing in the t sector, 

m t = 174.3 GeV, M SUS y = 2 TeV, /i = 200 GeV, M 2 = 200 GeV, 

X t = 0, A T = A b = A t , m~ g = 0.8 M SI jsy , (3) 

• the "gluophobic-Higgs" scenario, with a suppressed ggh coupling, 

mt = 174.3 GeV, M SUS y = 350 GeV, fi = 300 GeV, M 2 = 300 GeV, 
X t = -750 GeV, A T = A b = A t , m~ g = 500 GeV , (4) 

• the "small-o; c fj" scenario, with possibly reduced decay rates for h — > bb and h — > r + r _ , 

mt = 174.3 GeV, M SUS y = 800 GeV, /i = 2.5 M SUSY , M 2 = 500 GeV, 
X t = -1100 GeV, A T = A b = A t , m g = 500 GeV . (5) 

As explained above, for the sake of simplicity, Msusy is chosen as a common soft SUSY- 
breaking parameter for all three generations. 
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